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Bestimmung der Linien gleicher Hauptspannungssumme 
mittels Interferenzen gleicher Dicke *. 


Von A. Dose und R. Landwehr. 


1. Einleitung. In der Spannungsoptik werden die Spannungsverhiltnisse in Bauteilen, soweit 
es sich um elastische Probleme handelt, an einem durchsichtigen Modell in polarisiertem Licht 
untersucht. Ebene Probleme werden dabei bevorzugt. Die iiblichen Verfahren liefern im allge- 
meinen nur die Differenz der Hauptspannungen, die der Schubspannung proportional ist. Man 
erhalt ein Bild der Linien gleicher Schubspannung, der Isochromaten. Mit besonderem Aufwand 
kann man auch die beiden Hauptspannungen von Punkt zu Punkt einzeln bestimmen. Eine 
solche haufiger durchgefithrte Messung beruht z. B. auf der Anwendung des Interferometers nach 
Zehnder-Mach (Favre1). Doch ist es einfacher, noch eine zweite Grée zur Bestimmung heranzu- 
ziehen. Dazu geht man im allgemeinen von der Summe der Hauptspannungen aus. Die Kombi- 
nation mit ihrer Differenz liefert dann die Werte der einzelnen Hauptspannungen. Zur Feststel- 
lung der Spannungssumme sind verschiedene Verfahren bekannt geworden, und zwar rein experi- 
mentelle, halbempirische und rein theoretische *. Fast alle Verfahren verlangen aber einen relativ 
groBen Aufwand an Zeit und Mitteln. Am einfachsten ist noch die Bestimmung der Hauptspan- 
nungssumme aus der Messung der Dickenanderung dd der Modelle der Dicke d (aus Material mit 
den elastischen Konstanten E, jw) nach der Gleichung 

o+o,=—ES. ( 
Diesen Zusammenhang benutzt schon Mesnager °, der iiberhaupt die Bestimmung der Haupt- 
spannungssumme in die Spannungsoptik eingefiithrt hat. Die Dickenanderung wird entweder 
mechanisch [mit mechanischer Ubersetzung (Huggenberger), mit optischer Ubertragung (Coker), 
mit interferenzoptischer Anzeige (Mesnager usw.) oder elektrischer Anzeige] oder rein optisch, und 
zwar mittels Interferenzen, gemessen. Da die auftretenden Dickenanderungen im elastischen Be- 
reich maximal etwa von der GréBenordnung 1/100 mm sind, ergibt sich eine erforderliche MeB- 
genauigkeit von einigen 1/10 wu. Das ist aber gerade die Gréfienordnung der Lichtwellenlange. 
Deshalb sind die Interferenzmethoden fiir diesen Zweck besonders geeignet. 


2. Die bisher bekannt gewordenen Interferenzverfabren. Iu der Tabelle 1 sind die verschiedenen 
fiir die Bestimmung der Hauptspannungen, insbesondere der Spannungssumme, benutzten Inter- 
ferenzverfahren zusammengestellt. 

Unter ,,Beziehung“‘ ist die Gleichung angefiihrt, nach der die Bestimmung von 0, 0, bzw. 
o, +o, erfolgt oder zu erfolgen hatte. Schon Mesnager* hat Interferenzen zur direkten punkt- 
weisen Bestimmung der Dickenanderung benutzt, indem er die zwischen zwei Glasplatten, die an 
den Armen eines Fithlhebels angebracht sind, vorbeiwandernden Fizeauschen Streifen auszahlt. 
Handliche Interferenzfiihlhebel beschreiben spater Vose* und Schaid®. Schaid hat durch Be- 


* Vorlaufige Mitteilung: R. Landwehr und A. Dose, Naturwiss. 36 (1949), S. 342. 

1 M. H. Favre, Schweiz. Bauztg. 90 (1927), S. 291, 307; Rev. d’Opt. 8 (1929), S. 193, 241, 289; auch 
Diss. Ziirich 1929. ‘ 

2 E.G. Coker u. L. N. G. Filon, A Treatise on Photoelasticity. 5. 170 ff., Cambridge 1931; L. Féppi u. 
H. Neuber, Festigkeitslehre mittels Spannungsoptik, S. 31—67. Miinchen 1935; G. Mesmer, Spannungs- 
optik, S. 23, 73—85, 149—158. Berlin 1939, M. M. Frocht, Photoelasticity, Bd. II, 5S. 202 ff. New York 
1948; L. Féppl u. E. Ménch, Praktische Spannungsoptik, S. 39 ff. Berlin-Géttingen-Heidelberg 1950; 
M. Hetényi, Handbook of Experimental Stress Analysis, S. 874ff., 906 ff. New York u. London 1950; 
A, Kuske, Verfahren der Spannungsoptik, 5. 25 ff. Diisseldorf 1951; L. Villena, Fotoelasticidad, Madrid 
1943; A. Pirard, La Photoelasticité, Paris und Liittich 1947; H. T. Jessop u. F. B. Harris, Photo- 
elasticity, Principles and Methods, London 1949. aes t 

3 4, Mesnager, Congrés international des méthodes d’essai des matériaux de construction, Paris 1 
(1900), S. 149; Buda-Pesth-Congress of the International Association for Testing Materials 1901, Ann. des 
Ponts et Chaussées 8. Sér., 1, 4. Trim. (1901), S. 128. 

4 R. W. Vose, J. Appl. Mech. 2 (1935), A 99. 

5 R. J. Schaid, Master of Science Thesis, Northwestern Univ. 1947. 
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pee male 1 die Scharfe der Interferenzstreifen und damit die Genauigkeit 

Tank * hat das Verfahren von Favre? durch Hinzufiigung eines A/4-Plattchens derart modifi- 
ziert, da} er — allerdings auch nur punktweise — unmittelbar die Spannungssumme arcs 
konnte *. Sinclair® und Bubb ® haben das ganze Feld der Linien gleicher Hauptspannungssum . 
also der Isopachen direkt aufnehmen kénnen; Sinclair wieder mit dem Vedado Maal thisvides 
meter und Bubb mit einem Michelson-Interferometer. Durch Benutzung eines besonderen Mei: 
Methacryl-Harzes, des Lucites, mit extrem geringer Doppelbrechung konnte Sinclair dieses Ver. 
fahren weiter vereinfachen’. Maris® hat zuerst mit Hilfe einer einfachen Interferenz-Anordnun 
— und zwar auch schon iiber das ganze Gesichtsfeld hinweg — die Spannungssumme Honea 
indem er die Interferenzen gleicher Dicke zwischen dem Modell und einer Vergleichsflache oe 
Messung benutzt. Unabhangig von ihm hat spiter Tesar® einen gleichen Vorschlag gemacht 
- Frocht!® hat das Verfahren weiter durchgefiihrt und auf verschiedene praktische Probleme ee 
wandt. Auf die experimentellen Schwierigkeiten weist schon Maris hin. Frocht mu8 besolidors 
MaBnahmen treffen, um das Modell auch bei Belastung zur Vergleichsebene parallel zu erhalten 
Das ganze Bild der Isopachen wird direkt nur bei einer ganz ebenen Oberflache des Modells Balin 
bar. Doch kénnte man auch hier bei geringen Abweichungen das weiter unten beschriebene Diffe- 
renzverfahren anwenden. 


Fabry" hat ein Verfahren unter Benutzung der in dem Modell selbst entstehenden Interferenzen 
gleicher Dicke angegeben, um beide Spannungen getrennt zu bestimmen. Dieses mufte dazu 
allerdings so parallel sein, da ein deutliches Interferenz-Streifensystem entsteht. Da er ferner 
ein Material mit einer im Verhdltnis zur elastischen Verformbarkeit hohen spannungsoptischen 
Empfindlichkeit benutzt, namlich Glas, erhalt er infolge der Doppelbrechung ein doppeltes Strei- 
fensystem, das sich nur unter Benutzung eines Analysators trennen laB8t. Fabry hat dieses Ver- 
fahren im wesentlichen nur vorgeschlagen und die Ausfiihrbarkeit an einem einzigen Versuch be- 
statigt, aber noch nicht die Bestimmung der Spannungen an einem Beispiel wirklich durchge- 


fiihrt. 


3. Das neue Verfahren !?. Mit dem in der vorliegenden Arbeit angegebenen Verfahren wird die 
Bestimmung der Spannungssumme weiter vereinfacht, so daB sie beziiglich des zeitlichen und 
apparativen Aufwandes den iiblichen spannungsoptischen Methoden nicht mehr nachzustehen 
braucht. Es benutzt wie Fabry die in dem Modell selbst entstehenden Interferenzen gleicher 
Dicke, und zwar zunachst im reflektierten Licht. In derselben Weise wie in einer Luftschicht! 
interferieren jeweils zwei entsprechende an der Vorder- bzw. Riickseite der Platte reflektierende 
Lichtstrahlen. Dadurch, da8 man — ahnlich wie Sinclair — ein Material verwendet, das eine 
groBe elastische Verformbarkeit (kleinen Elastizitatsmodul) mit einer kleinen spannungsoptischen 
Konstante verbindet, z. B. Plexiglas, erhalt man nun bei nicht zu hohen Belastungen nur ein ein- 
ziges Streifensystem. Wenigstens konnte in allen beobachteten Fallen die Zahl der Ordnungen 
auf einen ausreichenden Wert gebracht werden, ohne da’ Doppelbrechung auftrat. Berechnet 
man die Zahl der Streifen, die bei der auBerst zulassigen Belastung an der Elastizitatsgrenze!4 
0 = 100 kg/cm? entstehen, so erhalt man fiir eine Plattendicke von 1 cm und 4 = 589 my rund 
eine Isochromate und 50 Isopachen. Bei der iitblichen Zahl der Isopachen (10 bis héchstens 20) — 


Tied nel aor) Seely 
F. Tank, Helv. phys. Acta 9 (1936), 5. 611. 
Siehe FuBnote 1 von Seite 73. 
Bei Favre und Tank fehlt das von der Dickenanderung abhingige Glied. 
D. Sinclair, Proceedings Tenth Semi-Annual Eastern Photoelasticity Conference December 9, 1939, 
S: D. Sinclair, Phys. Rev. (2) 57 (1940),S. 564; D. Sinclair and P. B. Bucky, Trans. Am, Inst. Mining 


Ge 
and Metallurgical Eng. 139 (1940), S. 224. 
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F. W. Bubb, Journ. Opt. Soc. Am. 30 (1940), 5. 297. 
D. Sinclair, Journ. Opt. Soc. Am. 30 (1940), S. 511. 
H. B. Maris, Journ. Opt. Soc. Am. 15 (1927), 5. 207. 
9 V. Tesar, Rev. d’Opt. 11 (1932), 5. 97. 
10 M. M. Frocht, Journ. Frankl. Inst. 216 (1933), S. 73; — Proc. 5th Int. Congress Appl. Mech. New 
York 1939, S. 221; Journ. Appl. Phys. 10 (1939), S. 248; vgl. auch FuBnote 2 von Seite 73. 
11 Ch. Fabry, C. R. Acad. Sci. 190 (1930), S. 457. 
12 In dem Buch von A. Kuske, (FuBnote 2 von §. 73), $. 35, Z.4 und 5 von oben, sind die hinweisenden 
Zahlen zu vertauschen. 
18 Vel. R. Landwehr u. G. Grabert, Ing.-Arch. 18 (1950), S. 1. 
14 4. Kuske gibt inseinem Buch S. 38 fiir eine Belastung von 5 bis 10 min. Dauer 150 kg/cm? als Auferste 
Grenze an, bevor merkliches Kriechen eintritt. 
6* 
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man arbeitet meist noch am Anfang des elastischen Bereichs — kann die Doppelbrechung also 
héchstens nur eine ganz geringe Verbreiterung der Streifen bewirken. Da die Interferenzerschei- 
nungen an das Modell gebunden sind, fallt ferner die durch die Unsicherheit der Lage, z. B. bei 
Frocht, sich ergebende Unbestimmtheit fort. Dadurch, da das Lichtbiindel direkt an der Riick- - 
seite des Modells reflektiert wird, ist die mégliche Versetzung zwischen dem hin- und zuriick- 
gehenden Strahl, die bei einigen anderen Verfahren auftritt!, auf ein Minimum zuriickgefihrt. 

Zur Erzeugung und Auswertung der Interferenzen ist unter Umstanden schon eine einfachere 
Anordnung ausreichend, wenn man die nétigen Bedingungen einhalten kann. Fir das Auftreten 
eines geniigend scharfen Interferenzstreifensystems ist erforderlich, da das monochromatische 
Licht méglichst senkrecht und ausreichend parallel auf das Modell fallt. Um in einer Platte von 
1 cm Dicke noch ein deutliches Streifensystem zu erhalten, verlangt die Koharenzbedingung bei 
einem senkrechten Auffall des Achsenstrahls, daB die Abweichungen von der Senkrechten nicht 
mehr als 14 Bogenminuten betragen. Bei einem Auffallswinkel von 1° ist die Toleranz sogar nur 
1,5 Bogenminuten®. Die Beobachtung erfolgte im reflektierten Licht. Natiirlich kann das Ver- 
fahren mit Hilfe einer geeigneten Anordnung auch in der schon von Fabry angegebenen Weise in 
durchfallendem Licht angewandt werden. Dazu ist das Modell durchlassig zu verspiegeln. Bei 
Benutzung von Mehrfachinterferenzen ist noch eine gréBere Genauigkeit zu erwarten. Es emp- 
fiehlt sich, die von Tolansky, allerdings fiir die Untersuchung kleinerer Objekte, angegebenen 
Hinweise und Forderungen zu beachten®. 


Abb. 1. Versuchsanordnung. 


4, Benutztes Geraét. Als besonders geeignet erwies sich das groBe Interferenzgerat der As- 
kania- Werke, Berlin, bei dessen Konstruktion auf die Méglichkeit der Beachtung der geforderten 
Bedingungen besonderer Wert gelegt wurde*. Als Lichtquelle dienten Spektrallampen der Firma 
Osram. Fiir die vorliegende Untersuchung wurde das Licht der roten Kadmiumlinie (A= 644 m 1) 
benutzt (Cd.-Lampe -+ Filter Schottglas RG 1, 2 mm), da man hiermit in dicken Schichten die 
scharfsten Interferenzen erhielt. Abb. 1 zeigt den Strahlengang der Versuchsanordnung. Die 
Bogensaule L der Lampe wird bei dem Gerat auf eine Irisblende Bin der Brennebene der groBen 
Beleuchtungslinse 0, abgebildet, die das Lichtbiindel mehr oder weniger parallel machen. Das 
Biindel fallt dann nach Reflexion an einer halbdurchlassig verspiegelten Flache S, von oben auf 
das Modell M. Nach Reflexion an dem Modell durchsetzt das Biindel den durchlassigen Spiegel 
und gelangt in den Beobachtungskopf. Mit Hilfe besonderer Ansatze kann das Modell zusammen 
mit dem Interferenzstreifensystem nun entweder — je nach Wahl des Okulars, mit 1,7- oder 
2,8facher VergréBerung — beobachtet oder im MaBstab 2 : 1 (auf Plattenformat 6 x 9 cm) auf 
der photographischen Schicht Ph aufgenommen werden. (Bei dem neuen Gerat erfolgt der Uber- 


1 7Z. B. F. W. Bubb, (FuBnote 6 von §.75), vgl. dazu D. Sinclair, Journ. Opt. S A 
5.567; F. W. Bubb, Journ. Opt. Soc. Am. 31 (1941), S. 506. yy cai oe 
2 R. Landwehr, Optik 5 (1949), S. 355. 
3S. Tolansky, Phys. Bl. 4 (1948), 5. 472; Multiple Beam-Interferometry, Oxford 1948, dort weitere 
Literaturangaben. 
4 R. Landwehr, Z. Instrumentenkde. 62 (1942), S. 73. Uber ei Gera i i 
4 cm freier Offnung vgl: Druckschrift Phyo. 131. ey EGlok viele Gate ee 
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ang von der visuellen Beobachtung zur photographischen Aufnahme einfach durch Umlegen 
ines Hebels.) Als zweckmaBig etwies sich auch die verzeichnungsfreie Abbildung des Gbjecies 
sesonders vorteilhaft zur Aufnahme der Kinspannvorrichtung ist der groBe freie Raum unterhalb 
les Beobachtungskopfes und der geraumige feinverstellbare Tisch. Zur sehr empfindlichen punkt- 
yeisen Ausmessung der Dickendnderung bei parallelen Modellen ist die Kinrichtung zur Unter- 
uchung mit Interferenzen gleicher Neigung (Haidingerschen Ringen) besonders gut brauchbar. 


5. Modell-Werkstoff. Ein Vergleich der spannungsoptischen Konstanten fiir verschiedene 
ebrauchliche Werkstoffe zeigt die nachstehende Tabelle 2. S, und S, geben die Differenz bzw. 
ie Summe der Hauptspannungen an (unabhingig vom Vorzeichen), die fiir eine Streifenordnung 
n = 1 bei einer Modelldicke von lem gilt. Die jeweils vorliegende wirkliche Spannung ergibt 
ich also durch Multiplikation mit der Zahl der Ordnungen und Division durch die Dicke (in em). 
p Wurde aus der Beziehung (1) unter Beriicksichtigung, daB 


5 mA 

od = oe (2) 
ilt (m Ordnungszahl der Interferenzstreifen), berechnet, wobei m= 1 und d =1 em gesetzt wur- 
en. Die elastischen Konstenten! und der Brechungsindex n (fiir 2 = 589 myc) sind ebenfalls in 
ie Tabelle aufgenommen. In der 6. Spalte sind auch noch die oft? angefiihrten Konstanten C 
=A/S., absol.) angegeben (10-7 cm2/kg ~ 1 brewster). S, gilt entsprechend der im allgemeinen tib- 
chen Anordnung fiir durchfallendes Licht, S, —wie im vorliegenden Falle fiir reflektiertes Licht 
fiir durchfallendes Licht sind die Zahlen zu halbieren). Man erkennt die geringere spannungs- 
ptische Wirksamkeit bei Glas und Plexiglas und die noch wesentlich geringere bei Lucite (Spalte 
und 7), dagegen die grofe interferenzoptische Empfindlichkeit bei allen Kunstharzen (Spalte 8). 
Besonders geeignet fiir die Feststellung der Spannungssummen ist also Lucite. Aber auch 
texiglas laBt sich, wie die vorliegenden Versuche zeigen, noch gut verwenden. Die Methyl- 
lethakrylharze haben allgemein vor den Kunststoffen der Phenol-Harzgruppe voraus, daB sie 
eringere Kigenspannungen haben und sich ihr Spannungszustand iiber langere Zeit hinweg un- 
erandert erhalt. Eine elastische Nachwirkung machte sich bei den Interferenzversuchen nicht 
emerkbar. Wegen der sich ausschlieBenden Forderungen ist es fiir genauere Messungen wohl 
nvermeidlich, fiir die Beobachtung der Isochromaten und der Isopachen zwei verschiedene Mo- 
elle zu benutzen. Schon zur Bestimmung der Isoklinen ist ein Modell aus einem optisch weniger 
mpfindlichen Material angebracht. Wegen seiner guten Eigenschaften hat sich gerade Plexiglas 
ierfiir als besonders geeignet erwiesen. Im vorliegenden Fall wurden fir die spannungsoptischen 
ersuche Modelle aus Dekorit oder Trolon und fiir die interferenzoptischen solche aus Plexiglas 
erwandt. Héchstens fiir angenaherte Untersuchungen unter Verzicht auf héhere Ordnung lieBe 
ch wahrscheinlich mit einem einzigen Modell aus einem Werkstoff mittlerer spannungs- und 
\terferenzoptischer Empfindlichkeit, wie z. B. Reilit, auskommen. Doch konnte dies bisher nicht 
aher untersucht werden. Reilit erwies sich titbrigens im spannungsoptischen Versuch nahezu 

nabhangig von der Dauer der Krafteinwirkung. 


6. Modell-Gestalt. Da die Modelle so planparallel wie moglich sein miissen, wurden die han- 
elsitblich gelieferten Plexiglasplatten mit Hilfe eines Priiffernrohrs der Firma Askania, das die 
eilwinkel mit einer Genauigkeit von mehreren Sekunden zu messen gestattet, aussortiert. Auf 
en ausgesuchten Platten wurden dann die Stellen mit dem geringsten Keilwinkel angezeichnet. 
ie weitere Uberpriifung der Platten exfolgte dann mit dem Interferenzger at. Abb. 2 zeigt einen 


1 Der Elastizitatsmodul ist abhingig von der Art und Dauer der Belastung. Z. B. gilt der in Zahlen- 
fel 1 fiir Plexiglas genannte Wert von 32000 kg/cm? (nach freundlicher Mitteilung der Fa, R6hm & 
aas) fiir eine Belastungsdauer von etwa 1 min. Bei steigender Belastungsdauer tritt eine merklich 
Bere Dehnung und damit eine Verminderung des Elastizitatsmoduls auf. In einem dynamischen Ver- 
ich ergeben sich véllig andere Werte, z. B. fiir Plexiglas 50000 kg/cm’. 

AuBerdem ist die Temperaturabhangigkeit der mechanischen und optischen Eigenschaften zu he- 
hten; z. B. ergeben sich nach G. H. Lee und C. W. Armstrong, Journ. Appl. Mech. 5 (1938), A. 11—12, 
ich Umrechnung folgende Anderungen pro 1° C 


BT 61—893 Marblette, ungehartet Marblette, gehartet 
E in kg/cm? 15 579 318 
S. in kg/em*?-cm 0,03 0,20 Ont 


2 7. B. L. Féppl, Kap.: Spannungsmessung an Werkstiicken (Ergebn. de techn. Rontgenkde., Bd. VI), 
eipzig 1938, S. 85; G. Mesmer, (Fr Bnote 2 von Seite 73), 5. 117; A. Peterlin und H. A. Stuart, Doppel- 
echung, insbes. kiinstliche Doppelbrechung (Hd.- u. Jahrbuch d. chem. Phys. Bd. 8), Leipzig 1943. 
Ib 15; A. Kuske, (FuBnote 2 von Seite 73), S. 37. 
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Ausschnitt aus einem Interferenzstreifensystem, das in einer brauchbaren Platte entsteht. Bei 
einem Keilwinkel von 1 Bogenminute betrigt die Zahl der auf 1 cm entfallenden Interferenz- 
streifen 13,5. Bei ebenen Begrenzungsflachen ist dann der Abstand zweier Streifen rund 0,74 mm. 
Als obere Grenze fiir den Keilwinkel kann etwa ein solcher von 3 Bogenminuten gelten. In diesem 
Fall 1a8t sich dann das unter Abschn. 8 angegebene Verfahren noch eben anwenden!. In Wirklich- 
keit sind die Platten natiirlich unregelmafig. Sie haben Sattelflachen und dergl. (vgl. Abb. 3). 
Die Plattenstiicke mit besonders geringem Keilwinkel wurden zur Anfertigung der Modelle benutzt. 

Bei der Ubertragung der wirklichen Verhaltnisse auf das Modell sind auch hier natiirlich die 
Abnlichkeitsgesetze zu beachten (vel. Ménch?). Fiir die Bestimmung eines FormanderungsmaB- 
stabes muB zumindest der Elasti- 
zitatsmodul bekannt sein. Er 1abt 
sich ahnlich wie im Fall der Iso- 
chromaten mit Hilfe eines Eich- 
versuches feststellen (s. u.). Bei der 
Bearbeitung eines Problems verfahrt 
man zweckmabigerweise so, da man 
mit Hilfe der bereits bekannten Ver- 
fahren und dieses neuen die Span- 
nungsverteilung zunachst einmal an 
dem Modell ausfiihrlich untersucht 
und erst dann die Ubertragung auf 
den wirklichen Zustand sinngemaB 
vornimmt. Die Bearbeitung der Mo- 
delle erfolgt unter Beachtung der 
bekannten VorsichtsmaBregeln unter == == 


méglichster Vermei une. You Vor Abb. 2. Interferenzstreifensystem in einer guten Kunststoffplatte (Reilit), 
spannungen, Dicke 5 mm, Ausschnitt. Vergr.: 1,7. 


<< Se 


7. Einspannvorrichtung. Zur Aufnahme der Modelle wurde ein besonderer rechteckiger 
Doppelrahmen gebaut (etwa 25 x 15cm), der an den Ecken unter Verwendung von Distanz- 
stiicken zusammengehalten wurde. In diesem konnte dann auf der einen Seite die Haltevorrich- 
tung fiir das Modell in geeigneter Weise mit Schrauben befestigt werden, wahrend auf der gegen- 
iiberliegenden Seite und seitwarts Fiihrungen fiir eine Belastungsschraube angebracht waren. 
Die GréBe der Belastung wurde durch Eichung ermittelt und konnte an einer Trommel abgelesen 
werden. Dieser Rahmen wurde zusammen mit dem eingespannten Modell auf den Tisch des Inter- 
ferenzgerates gelegt, der zur Vermeidung jeglichen Reflexlichtes mit einem schwarzen Velours- 
papierstiick belegt war. Die Justierung geschah dann unter Betatigung der beiden vorne an dem 
Tisch befindlichen Feinstellschrauben. Das Modell stand senkrecht zu dem auffallenden Strahlen- 
biindel, wenn das Reflexbild der genannten Irisblende sich in der Mitte der Aufnahmedffnung be- 
fand. Bei dem benutzten Gerat konnte eine genaue Kontrolle mit Hilfe des Ansatzes zur Beob- 
achtung Haidingerscher Ringe erfolgen. Das Blendenbild erschien direkt in dessen Okularge- 
sichtsfeld. — Da der Durchmesser des Gesichtsfeldes des Interferenzgerates nur 120 mm betrug, 
muBten bei gréBeren Modellen mehrere aneinander anschlieBende Aufnahmen gemacht werden. 
Das Modell lieB sich dazu mit dem Tisch, der mit Kugellagern auf Zylindern lauft, verschieben 
und in jeder Stellung festklemmen, ohne daf} eine neue Justierung erforderlich gewesen ware. 


8. Punktweise Bestimmung der Dickenanderung. DieBestimmung der Dickenanderung 6d 
des Modells 1aGt sich an jeder belicbigen Stelle durch visuelle Beobachtung vornehmen. Dazu zahlt 
man nur die Zah] m der an dem betreffenden markierten Punkt unter dem Einflu8 der Belastung 
vorbeiwandernden Interferenzstreifen. Die Form des Ausgangsinterferenzsystems ist dabei voll- 
kommen gleichgiiltig. Die Dickenanderung bestimmt sich dann nach (2) fiir Plexiglas und 
, = 644my zu: 2,16 x 10->.mem. Die Spannungssumme ist 


0+ 02 = Spy (3) 


1 Die Einbeziehung von Platten mit noch gréSeren Formfehlern lat sich wabrscheinlich durch die 
Einbettung der Modelle in eine geeignete Immersionsflissigkeit erreichen. E. Zehender hat ein solches 
Verfahren mit groBem Erfolg zur Oberflachenpriifung im mikroskopischen Bereich verwandt (Z. VDI 94 


1952) S. 456). ‘ 
: 2 i Arce Z. angew. Phys. 1(1949), S$. 306; L. Féppl und E, Ménch (FuBnote 2 von Seite 73), 5. 74ff. 
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(S, s.Tabelle 1, din cm). Mit Hilfe eines Okularschraubenmikrometers lassen sich natiirlich auch 
Zwischenwerte als Bruchteile eines Streifenabstandes bestimmen. 

Die Bestimmung des Vorzeichens der Dickenanderung, also die Bestimmung, ob an der be- 
treffenden Stelle Zug (+) oder Druck (—) vorliegt, ist leicht, wenn die Richtung der Dickenzu- 


nahme an dieser Stelle bekannt ist. Wenn es sich — wie in den vorliegenden Fallen — um einen | 
gréBeren Keilwinkel zwischen der oberen und unteren Begrenzungsflache (von einigen Minuten) | 


handelt, kann dieser u. U. mit einer Mikrometerschraube vorher an dem Modell bestimmt werden. 
Eine optische Messung des Keilwinkels 1aBt sich mit dem Priiffernrohr vornehmen (vgl. Abschn.6). 


Bei einer geeigneten optischen Interferenzanordnung, wie Z. B. bei dem genannten Interferenz- | 
gerit, bei der die von den beiden Begrenzungsflachen reflektierten Strahlenbiindel ein beobacht- | 


bares Bild einer definierten Blende oder Marke liefern, 148t sich auch aus deren Lage zueinander 
die Lage des Keilwinkels bestimmen. Nur muf vorher festgestellt werden, welche Flache an dem 
Zustandekommen des betreffenden Reflexbildes beteiligt ist. Zu diesem Zweck kann man die 
Reflexion der unteren Flache herabsetzen bzw. ganz ausléschen, indem man die Flache anhaucht 
oder ein Stiick schwarzen Seidenpapiers mit einem Oltropfen andriickt. 

Bewegen sich nun bei Belastung die Interferenzstreifen in Richtung auf das diinnere Ende der 
Platte zu, so bedeutet das eine Dickenzunahme. An der betreffenden Stelle herrscht also eine 
Druckspannung. Wandern die Interferenzstreifen dagegen zu den dickeren Stellen, so zeigt das 


eine Abnahme der Dicke, also das Vorhandensein einer Zugspannung an. Da sich das Streifen- | 
system als Folge der Belastung héchstens um etwa 20 Streifen, im allgemeinen aber um weniger, | 
verschiebt, laBt sich die Frage, ob die Streifenordnung zu oder abnimmt, meistens leicht aus dem || 
Vergleich mit dem Streifenverlauf im Ausgangszustand entnehmen. Ist das Vorzeichen der | 
Dickenadnderung an einer Stelle des Modells von voinherein aus anderen Uberlegungen bekannt, | 


so kann man von dort ausgehend die Vorzeichen an den iibrigen Stellen bestimmen. Beim Uber- 
schreiten einer Null-Isopache (Sd =0) kehrt sich das Vorzeichen um. Sind die Modelle sehr gut 
parallel, d. h. ist ihr Keilwinkel kleiner als 30 Bogensekunden, so aft sich eine sehr genaue punkt- 
weise Messung der Gré®e und des Vorzeichens der Dickenanderung bei dem benutzten Gerat mit 
Hilfe des Ansatzes zur Beobachtung der Haidingerschen Ringe vornehmen. Die oben erwahnte 
Blende dient dann zur Einengung der MeBstelle. Bei einer Zunahme der Dicke um A/2n quillt 
gerade ein Ring aus dem Zentrum hervor. Mit dem Okularschraubenmikrometer sind in diesem 
Falle noch ganz geringe Dickenunterschiede (1 x 1077 cm) zu messen. 


9. Ermitilung des ganzen Feldes der Isopachen. Das ganze Feld der Isopachen wirde man in 


dem Modell direkt erhalten, wenn es so vollkommen planparallel hergestellt werden kénnte, dafi | 
es ohne Belastung fast gar keine Streifen mehr zeigen, also im wesentlichen nur hell oder dunkel | 


erscheinen wiirde. Bei Belastung wiirden dann die Interferenzen gleicher Dicke direkt die Linien 
gleicher Dickenanderung = Isopachen darstellen. Die Herstellung eines derart optisch plan- 
parallelen Modells verlangt aber einen auBerordentlichen Aufwand. Die zu diesem Zweck vorge- 
nommenen Versuche iiber eine weitere Bearbeitbarkeit der Flachen konnten jedoch nicht fortgesetzt 
werden. Sie zeigten zunachst nur, daf es (unter Benutzung des Poliermittels Plexipol) gelang, 
gut ebene Flachen zu polieren, wobei allerdings die Parallelitat der Platte verloren ging. 


Da aber nur die Dickenanderung von Interesse ist, kann man auch fiir den Fall, daB die | 


Modellplatte einen kleinen Keilwinkel hat oder die Begrenzungsflachen titberhaupt unregelmaBiger 


sind, einen Kunstgriff anwenden. Voraussetzung ist nur ein einigermafen deutliches Interferenz- 
streifensystem im Ausgangszustand. Man kombiniert eine photographische Aufnahme des Mo- 


dellsim belasteten Zustand (verschobenes Interferenzsystem) mit einer solchen des unbelasteten 
Modells (unverschobenes Interferenzsystem). Allerdings diirfen die Gestaltanderungen nicht zu 


grof sein, sonst tritt bei der Uberlagerung eine Verschiebung des Isopachensystems auft. Bei der | 


1 Bei den bisherigen Belastungsversuchen, die weiter ausgewertet wurden, traten seitliche Ver- 
schiebungen von max. etwa 0,13 mmein, Da die niedrigen Streifenordnungen meist einen grofen 


Abstand haben und sich nur die héheren Ordnungen enger zusammendringen, ist der prozentuale — 


Fehler fast immer von derselben GroBenordnung und betragt nicht mehr als 5%, meist sogar wesent- 
lich weniger. Nur in den wenigen ganz ungiinstigen Fallen, in denen auch die Streifen kleiner Ord- 
nung sehr geringen Abstand haben und zudem senkrecht zur Verschiebungsrichtung liegen, kann ein 
Fehler von einer halben Ordnung auftreten. — Die schnelle Wanderung des Systems der Uber- 
lagerungsstreifen bei Verschiebung der beiden zur Deckung gebrachten Bilder ist nur eine Eigenschaft 
der Moirée-Streifen, die auch anderweitig benutzt werden kann (z. B. H. Voigt ,,Optik‘s 9 (1952) 
5: 438), und bedeutet nur einen Wechsel von hellen zu dunklen Streifen, hat aber — bis auf den 
Mad Seegte a allenfalls auftretenden Einflu8 — keine grundsatzliche Anderung des Streifensystems 
zur Folge. 


| 
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starken Deformierbarkeit der Kunstharze gentigen im allgemeinen schon geringere Belastungen, 
um eine ausreichende Dickenanderung zu erreichen. Im vorliegenden Fall wurde meistens so ver- 
fahren, daB man von einer der Aufnahmen — vornehmlich von der des unbelasteten Zustandes — 
eine Kontaktkopie herstellte und diese dann (Schicht auf Schicht) mit der Aufnahme des be- 
lasteten Zustandes zur Deckung brachte. Als Kriterium kann das Aufeinandertreffen der Abbil- 
dungen eines Kratzers oder dergleichen dienen. Auf der VergréBerung in einem VergréBerungs- 
apparat erscheint dann deutlich ein System von Uberlagerungsstreifen (Moirée-Effekt), das un- 
abhangig von der Form des zugrundeliegenden Interferenzstreifensystems ist’, Die Uberlage- 
rungsstreifen sind besonders deutlich bei einem relativ engen, aber noch scharfen Interferenz- 
streifensystem. Sie sind ferner um so scharfer, je kleiner ihr Abstand ist und je mehr sie zur 
Richtung der Interferenzstreifen senkrecht verlaufen. Abb. 3 zeigt eine solche VergréBerung, 


Abb. 3. Isopachensystem in einem Balken auf zwei Stiitzen mit Mitteleinzellast, Ausschnitt. WVergr.:) 17x. 


und zwar das eine Ende eines Balkens auf zwei Stiitzen mit einer Einzellast in der Mitte als Be- 
lastung. Nach einer solchen VergréBerung wurde die in der vorlaufigen Mitteilung? gebrachte 
Zeichnung angefertigt, nur wurde dort der Ubersichtlichkeit halber jede zweite Isopache fortge- 
lassen. Die Modellplatte hatte die Dimensionen 40 x 130 X 7mm. Man erkennt deutlich in dem 
Ausgangsinterferenzsystem zwei Sattelflachen, die aber das symmetrische Uberlagerungssystem 
mm keiner Weise stéren. Das auftretende Isopachennetz ist also unabhangig von dem Abstand der 
Interferenzstreifen. 

Daf das durch Uberlagerung entstehende Streifensystem wirklich ein System gleicher Dicken- 
anderung, also ein System von Isopachen ist, erkennt man direkt aus der folgenden Abb. 4. Die 
Ordnung eines mittleren Streifens des Ausgangsinterferenzsystems sei k (die Dicke also k-A/2n), 
die Ordnung eines mittleren Streifens des durch Belastung verschobenen Systems sei k +m 
(Dickendnderung also m-A/2n). Die Ordnungszahl der Kreuzungsstellen beider Systeme ist 
jeweils die Differenz der Ordnungen der beiden sich schneidenden Interferenzstreifen. Daraus 
ergibt sich eine Schar von Linien gleicher Dickendifferenz (= Isopachen) mit den Ordnungen 
...m—l1,m,m-+1... (unabhangig von der Ordnung k des Ausgangssystems). Die Dicken selbst 
bzw. deren Unterschiede erhalt man daraus durch Multiplikation mit 4/2n (n Brechungsindex des 
Modellmaterials). Bei der Auswertung mu man darauf achten, in welcher Weise die beiden Auf- 
nahmen zur Deckung gebracht werden. Man muBf die folgenden Méglichkeiten unterscheiden: 
a) bei direkter Betrachtung von 2 Originalen sind die Isopachen die hellen Streifen, b) bei Ver. 


1 Nach freundlicher Mitteilung von Herrn Prof. Mesmer (St. Luis/ USA) und Dr. Drouven (Frankfurt), 
kann man vereinfachender Weise beide Zustande auf derselben photographischen Platte aufnehmen: 
2 R, Landwehr und A. Dose (FuBnote * von Seite 73). 
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gréBerung von zwei itbereinanderliegenden Originalen sind die Isopachen die schwarzweifen 
(hellen) Streifen zwischen dem dunklen Ubergangsgebiet, c) bei direkter Betrachtung der 
Kombination eines Originals mit einer Kopie liegen die Isopachen zwischen den hellen Uberlage- 
rungsstreifen, d) bei VergréSerung dieser Kombination liegen die Isopachen im Ubergangsgebiet, 
das sind aber hier die hellen Streifen (zwischen den schwarzweiBen Gebieten). 


10. Ermittlung des Streifenwertes. Sind die 

> elastischen Konstanten geniigend genau bekannt 
oder nach anderen Verfahren zu bestimmen, so 
1a4Bt sich der interferenzoptische Spannungswert S, 
direkt mit Hilfe der Gleichungen (1) und (2) fiir | 
m=1 und d=1lem berechnen (vgl. Tabelle 1). 
Fir A = 644 my ist S,= 1,96 kg/em?-cm. Dabei 

ist aber auf die erwahnte Abhangigkeit der Kon- 
stanten von Temperatur und Belastungsdauer zu 
achten. Deshalb ist es zweckmaBiger, einen Kich- 
versuch anzustellen, indem man eine einfache, 
definierte Belastungsprobe an einem Modell aus 
demselben Material und von derselben Dicke unter 
denselben Bedingungen wie das Hauptmodell vor- 
Abb. 4. Nachweis der Uberlagerungsstreifen als Isopachen. nimmt und fiir geniigende Temperaturkonstanz 
wihrend des ganzen Versuches sorgt. Im Rahmen | 

der vorliegenden Untersuchung wurde ein Versuch mit einem Zugstab gemacht, der in dem- 
selben Einspannrahmen belastet werden konnte. Die Dimensionen des Stabes waren: Lange 
116 mm, Breite an den beiden Einspannenden 30 mm, Lange des geraden Mittelteiles ungefahr — 
A 70mm und Breite in der Mitte 

| | 14,8mm. Der Stab war aus der- 

selben Platte hergestellt wie das 
Modell desHauptversuches(Platte 
801, Dicke 8mm). Eine visuelle 
Messung ergab fiir den Streifen- 
wert bezogen auf die Dicke 1 cm 
bei A = 644. my als Mittel aus zwei 
Einzelmessungen 1,99 kg/cm? (fiir 
d = 0,8 cm: 2,49 kg/cm?). Die 
Streuung (von + 15%) der ein- 
zelnen an verschiedenen Stellen 
aufgenommenen Mefwerte deutet 
darauf hin, daB es auch bei der | 
Kichung nach Isopachen mit | 
einem Zugstab sehr darauf an- 
kommt, da Belastung und Mes- 
sung genau auf der Mittellinie | 
vorgenommen werden. 


Abb. 5 zeigt Ausschnitte aus | 
dem in einem Zugstab entstehen- 
den Isopachensystem bei drei ver- 
schiedenen Belastungsstufen. Er- 
ganzt man das Isopachensystem 
bis zu den Nullpunkten an den 
scharfen Ecken,wo die Verjiingung 
beginnt, so erhalt man fiir einen 

Abb. 5. OGRA ee en Zugstabes Punkt in der Mitte auf dem Stab 

einen Streifenwert von 2,32 kg/cm? 

bei 0,8 cm Dicke bzw. von 1,86 kg/cm? bei lem Dicke. Tragt man die Zahl m der Ordnungen 
in Abhangigkeit von der Belastung auf, so erhalt man eine Gerade durch den Nullpunkt, also 
eine sehr gute Proportionalitat (Abb. 6). Bei bekanntem Spannungswert und bei bekannter 
Poissonscher Zah] kann man natiirlich mit Hilfe der genannten Beziehung aus der gemessenen 
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Spannungssumme den Elastizitatsmodul bestimmen. Aus dem gemessenen Streifenwert findet 
man im vorliegenden Fall einen Elastizitatsmodul von rund 30 200 kg/em?.1 


11. Einige allgemeine Eigenschaften der Isopachen 2. Entsprechend der Gréfe der spannungs- 
optischen Konstanten (vgl. Tabelle 2) treten bei derselben Randspannung im Fall der Isopachen 
weniger Streifen auf als im Fall der Isochromaten. AuBerdem e1 folgt die Zu- oder Abnahme auf 
Grund ihrer Bestimmungsgréfen im allgemeinen auch weniger steil, als das bei den Isochromaten 

der Fallist®, Kine besonders wichtige Rolle spielen die Null-Isopachen, also die Linien, langs deren 
die Summe der Hauptspannungen verschwindet, d.h. also, wo die beiden Hauptspannungen 
entgegengesetzt gleich sind. Es sind die Linien, die die Gebiete, in denen Zug und Druck vor- 
herrschen, voneinander trennen. 

In dem Momenten-Nullpunkt eines rechteckigen auf Biegung beanspruchten Balkens ist 
bei verschwindender Normalkraft stets die Spannungssumme Null, wie sich leicht aus der 
Bildung der Spannungssumme nach den entsprechenden Ansiitzen ergibt ¢. 

A Auch im Isopachennetz sind singulire Punkte aus- 

gezeichnete Punkte. Inihnen stehen Isochromaten und 
Isoklinen (Linien gleicher Hauptspannungsrichtung) 
aufeinander senkrecht, falls nicht beide Hauptspan- 


20 


4 
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Abb. 6. Abhangigkeit der Ordnungszahl m der Abb. 7. Isopachensystem in einem Balken zwischen zwei ebenen 
Isopachen von der Belastung (Nennlast). Belastungsstempeln. 


nungen einzeln gleich Null sind. Vorspringende Ecken sind sowohl singulare Punkte der Haupt- 
spannungslinien als auch der Isopachen. Im allgemeinen fallen aber singulare Punkte der Iso- 
pachen nicht mit solchen der Hauptspannungslinien zusammen. Man unterscheidet auch bei 
den Isopachen singuldre Punkte verschiedener Ordnungen. Die singularen Punkte 1. Ordnung 
haben vier aufeinander senkrechte ‘Asymptoten an die Isopachen. Singulare Punkte zweiter 
Ordnung haben 6 verschiedene unter 60° versetzte Asymptoten. 

Ist ein freier Rand gradlinig, so haben Isochromaten und Isopachen die gleiche Richtung (die 
Isokline fallt mit dem Rand zusammen). Steht die Isokline aber auf dem Rand senkrecht, so 
verlaufen die Isochromaten und die Isopachen parallel zum Rand. An diesen Stellen hat die 
Randspannung o, ihren gréften (positiven oder negativen) Wert, wahrend die Spannung o, an 
den lastfreien Randern verschwindet. 


12. Beispiele. Zur Orientierung wurde das System der Uberlagerungsstreifen ‘n einigen Fallen 
erzeugt, bei denen der Verlauf des Isopachensystems bereits weitgehend bekannt ist. Abb. 7 zeigt 
das Isopachensystem, das in einem Balken (Modell aus Plexiglas) zwischen zwei ausgedehnten 
Lasten auftritt. Da ein harter Stahlstempel mit ebenen Endflachen (von 16 mm Breite) be- 
nutzt wurde, greift die Last nicht gleichmafig ither den Querschnitt der Angriffsflache an, son- 
dern nimmt nach den Enden der Laststrecke zu’. 


1 Siehe FuBnote 1 von Seite 77; H. Peukert gibt fiir ungerecktes und unbelastetes Plexiglas vom Typ 
M 33 an: E = 30400 kg/cm? [Z. VDI. 93 (1951), S. 382]. 

2 Vel. M. M. Frocht, (FuBnote 2 von Seite 73), insbes. aber L. Féppl und H. Neuber, (FuBnote 2 von 
Seite 73). 

3 M. M. Frocht, Journ. Appl. Phys. 10 (1939), 5. 255. 

4 G. Mesmer, (Fufnote 2 von Seite 73), S. 178 ff; L. Féppl und E. Ménch, (FuBnote 2 von Seite 73), 
5. 46 ff. 

5 E.G. Coker und L. N.G. Filon, (FuBnote 2 von Seite 73), 5.354; G. Mesmer, (FuBnote 2 von Seite te) 
S.197; M. M. Frocht, Journ. Frankl. Inst. 216 (1933), S. 82ff. Bei gleichmafigem Lastangriff wiirden 
die zwischen den Stempeln liegenden Isopachen zu den Enden der Laststrecke hin gehen, wie eine bereits 
vorliegende theoretische Durchrechnung zeigt. 
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Die Aufnahme der in einer Halbebene unter der Einwirkung einer Einzellast entstehenden 
Isopachen (Reilitplatte, Belastung itber eine Schneide) ergab in Ubereinstimmung mit der 
Theorie! in der Umgebung des Lastangriffspunktes Kreise durch diesen Punkt. Allerdings lieferte 
ein Randstreifen wegen zu groer Dickenunterschiede keine Interferenzstreifen. Von dem Fall 
des Balkens auf zwei Stiitzen mit Mitteleinzellast wurde schon in der vorlaufigen Mitteilung ein 
Bild des auftretenden Isopachensystems gebracht. In diesem Bild befinden sich singulare Punkte 
an der freien Ecke, wo 0, +0 
—Qist, und im Mittelpunkt 
der Sattelfliche in der Nahe 
dieser Ecke. 

13. Quantitatives Beispiel: 
gelochter Triger. Die Brauch- 
barkeit des Verfahrens wurde 
ferner noch an einem Beispiel 
nachgewiesen, namlich an dem 
Fall eines gelochten Tragers, 
der auf der einen Seite ein- 
gespannt war und durch eine 
Querkraft am freien Ende auf 
Bicgung beansprucht wurde. 
In diesem Falle wurde eine 
Abb. 8. Isopachensystem in einem durch eine Querkraft beanspruchten, einseitig rechnerische Kontrolle durch- 
eingespannten gelochten Riese ocen fake A ee an cinem gefihrt. Das Modell wurde aus 

einer der besseren Plexiglas- 
platten von 8 mm Dicke, aus der auch der Kichstab entnommen war, hergestellt (Dimensionen: 
Lange — von der Einspannstelle aus — 107 mm, Breite 41 mm, drei Lécher von 25 mm Durch- 
messer in 33mm Abstand). Interferenzaufnahmen wurden von dem Modell auBer im un- 
belasteten Zustand bei drei verschiedenen Belastungsstufen (bei den Trommelstellungen 2, 10 
und 24 der Belastungsschraube) hergestellt. Leider waren nur in der einen Halfte des Modells 
die Interferenzstreifen einigermaBen 
deutlich und daher die Isopachen 
auch nur dort sichtbar. Abb. 8 zeigt 
eine Zeichnung des durch das be- 
schriebene Uberlagerungsverfahren 
bei der mittleren Belastungsstufe 
entstandenen Isopachensystems. 
Hine Hilfe bei der Festlegung des 
Isopachenverlaufes im einzelnen ist 
x hier die bekannte Tatsache, dai am 
~~ freien Rande der Lécher Zug und 
Druck abwechseln, also vier Null- 

Abb. 9a. Zur Berechnung der Werte Abb. 9b, Zur Berechnung der Rand- . ss 
eines Potentialfeldes nach Gl.(4).  werte eines Potentialfeldes nach Gl. (5) punkte vorhanden sein miissen. Auch 
lassen sich an dem Isopachenbild 
schon von vorneherein Stellen angeben, an denen Zug, und solche, an denen Druck herrschen 
muff. Aus den Aufnahmen bei verschiedenen Belastungsstufen sieht man aufSerdem, an 
welchen Stellen sich das Isopachensystem am wenigsten und an welchen Stellen es sich am 
meisten andert. Auch die Festlegung der Ordnungen wird durch Aufnahmen bei mehreren Be- 
lastungsstufen erleichtert (Mesmer betont diesen Vorteil bei der Bestimmung der Isochromaten?). 
In der Mitte der Sattelflachen zwischen den Léchern liegen s'‘ngulaére Punkte 1. Ordnung. Die 
GiéBe der Hauptspannungssumme 1a4B8t sich nun an jeder Stelle aus der angegebenen Ordnung 
durch Multiplikation mit dem vorher (Abschn. 10) gefundenen Streifenwert bestimmen. (Dieser 

betragt fiir eine Dicke des Modells von 0,8 cm : 2,32 kg/em?.) 

Um nun feststellen zu kénnen, wieweit das experimentell gefundene Isopachenbild mit dem 
wirklich vorhandenen iibereinstimmt, wurde das Isopachenbild herangezogen, das sich auf Grund 
der experimentell bestimmten Randwerte der Isochromaten berechnen laBt. Diese Berechnung be- 


1G: Mesmer, (FuBnote 2 von Seite 73), S.195; D. Sinclair and P. B. Bucky, (FuBnote 5 von Seite 75), 
Se); ant), 
2 G, Mesmer, (FufBnote 2 von Seite 73), S. 139. 
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itzte die Tatsache, daB die Linien gleicher Hauptspannungssumme Aquipotentiallinien dar- 
ellen, daB also die Gesetze der Potential-Theorie angewandt werden kénnen. Vor anderen gra- 
uischen Methoden! hat die Rechnung den Vorzug, sich ohne besondere Schwierigkeiten durch- 
hren zu lassen. Auch ein neuerdings von A. Kuske? angegebener zeichnerischer Weg fiihrt erst 
ich der Konstruktion verschiedener Kurvenscharen zum eigentlichen Ziel. 


Bei einem quadratischen Netz von Feldpunkten (Abb. 9a) gilt niherungsweise 
UW (x,y) = 1/4 (u, + uy + uz, +uy,). (4) 


iegt nun als Rand eine gezeichnete oder analytisch gegebene Kurve vor, so muB folgende ab- 
-wandelte Differenzengleichung verwendet werden (Abb. 9b) 


gee elena a : 
POS SESS SUEY eee 


gl. hierzu G. H. Shortley und R. Weller*). Da die Laplacesche Gleichung auch in mehrfach 
isammenhangenden Bereichen Giiltigkeit besitzt, so ist die Anwendung der obigen Gleichung in 
*m vorliegenden Falle gestattet ¢. 


& 1 
Mo Speers lili (lease) U4 (5) 


Abb. 10. Isopachensystem in einem durch eine Querkraft beanspruchten einseitig 
eingespannten gelochten Trager (errechnet auf Grund der Randwerte nach dem 
Isochromatenbild eines Trolonmodells). 


: Durch wiederholte Anwendung der Gleichungen (4) und (5) wurde dann, ausgehend von den 
oh aus den Isochromatenmessungen ergebenden Randbedingungen, unter Zugrundelegung 
nes geeigneten Netzes von Punkten, das ganze Feld der Hauptspannungssumme aufgebaut. 
ie GréBe des Netzes ist von Bedeutung, sowohl im Hinblick auf die Konvergenz des Verfahrens, 
s selbstverstandlich auch auf den Umfang der auszufiithrenden Rechenoperationen. Priifung 
:s Momentengleichgewichts und der Querkrafte in verschiedenen Schnitten diente zur 
ontrolle. 


_ Abb. 10 zeigt das Ergebnis dieser etwas langwierigen Rechnungen. Ihnen lagen Versuche an 
nem Trolon- und e‘nem Dekorit-Modell zu Grunde. Die Ubereinstimmung des gerechneten 
reifensystems mit dem interferometrisch an einem Plexiglasmodell gewonnenen System (Abb. 8) 
: durchaus befriedigend. Die Ubereinstimmung in der Zah] der Streifen ergab sich zufiallig fiir 
m gerechneten Fall und den empirisch bestimmten Fall der mittleren Belastungsstufe. Die 
rdnungszahlen stimmen fast genau iiberein. Die Abweichungen in den Ordnungszahlen liegen 
1 Mittel unter 4°/. Wesentliche Abweichungen kommen nur am oberen Rand in der Nahe der 
slastungsstelle und an der rechten AuBenkante vor, was wahrscheinlich z. T. auf unterschied- 
she Eigenspannungen in den beiden Modellen (Plexiglas bzw. Trolon) zuriickgeht. Bei spannungs- 
‘tischen Versuchen wird man im allgemeinen mit Fehlern von 5 bis héchstens 10% rechnen 


1 Vgl. FuBnote 2 von Seite 73, insbesondere L. Féppl und H. Neuber. 

2 A, Kuske (FuBnote 2 von Seite 73), S. 28 ff. 

3 G. H. Shortley and R. Weller, Journ. Appl. Phys. 9 (1938), 5S. 334, 348. 

4G. H. Shortley and R. Weller, Journ. Appl. Mech, Trans. ASME 6 (1939), S, A-71. 
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kénnen, so daB die vorgeschlagene Methode hinsichtlich ihrer Genauigkeit den berechtigten Fol; 
derungen entsprechen wird. Einige bereits durchgefiihrte diesbeziigliche Versuche bestatigte, 
diese Genauigkeit des Verfahrens. Auf die Wiedergabe muf infolge des Verlustes einiger Unte 
lagen im Augenblick verzichtet werden. Es sei erwahnt, daB die Héhe der nach dem bb 
schriebenen Verfahren erhaltenen Randspannungen in der GréSenordnung mit den Werte 
iibereinstimmte, die von anderen nach einem mechanischen Verfahren gewonnen wurden 1, | 

Die Felder der einzelnen Spannungen o, und o, erhalt man endlich in der itblichen Weis 
durch Kombination des Isochromatennetzes mit dem Isopachennetz ?. 


(Eingegangen am 7. August 1952.) 


Anschrift der Verfasser: Dr. A. Dose, Aachen, EifelstraBe 20; 
Dr. R. Landwehr, Menden/Sauerland, Schattweg 36. 


1 EK. Siebel und E. Kopf, Forsch. Ing. Wes. 5 (1934), Forsch. Heft 369. 
2 E.G. Coker und L, N. G. Filon, (FuBnote 2 von Seite 73), S. 178. 
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| Die Kreisplatte auf elastischer Unterlage 
unter einem an starrem Mittelstiick angreifenden Moment. 


Von H. S. Gedizli!. 


_ Es sei nach Abb. 1 eine elastische, isotrope, ebene Kreisplatte vorgegeben, deren konstante 
icke h klein sei gegen ihren Halbmesser ro. In der Mitte habe sie ein starres Stiick vom Halb- 
esser r;, an dem ein Drehmoment M angreife. Sie ruhe auf einer elastischen Unterlage, deren 
eaktionskonstante (Baugrundziffer) c sei ®. In der Liésung dieser Aufgabe, wie sie Fliigge 3, 
tmoshenko * und spater Beyer® angaben, verblieb eine unbestimmte Konstante, auBerdem 
eicht diese Lésung fiir kleine Dickenverhiltnisse h/r) und groBe Werte des Parameters 


Cc 
wt = TE (1) 
theblich von der Wirklichkeit ab. Dabei ist die GréBe D die allgemein gebrauchliche Platten- 
eifigkeit. Die im folgenden gegebene Lésung fiihrt rascher und genauer zum Ziel. 

Die Differentialgleichung der Durchbiegung w einer unbelasteten Kreisringplatte auf elasti- 
ther Unterlage lautet® in Polarkoordinaten r, 0 


0 | 1 @ i @ 
or? r2 Op? | r Or 


4 
AAw+Sw=0 mit A= Se ee (2) 
0 


le 
gant (3) 
e dimensionslose Veranderliche @ ein, geht sodann in die entstehende Gleichung mit dem Ansatz 


w (o, 0) = V(e) cos 0 (4) 


n und formt die entstehende gewohnliche Differentialgleichung fiir V (@) in der von Kirchhoff 7 
agegebenen Weise um , dann erhalt man als allgemeine Lésung 


V(o) = Ay Jive) + Ay Siz (0) + As ii (Q) + Aa Cie (Q) (5) 
0 J¥ (0). ...+ Gite (@) die von Télke® eingefithrten Kombinationen von Zylinderfunktionen sind, 
ad die Integrationskonstanten A; mit Hilfe dreier Randbedingungen und einer Gleichgewichts- 
edingung bestimmt werden miissen. 
- Unter Beriicksichtigung der Ersatzscherkraft nach Kirchhoff hat man am freien Rand der 
latte die beiden Randbedingungen ° 


‘ihrt man in sie mittels 
: 


| m, (tT) =9, | (%) = 9 (6) 

e fiir w auf 

2 2 

(Gert RE) mo 9 

ad 

| & Lew. low -3—yp Cw , 2—pn Ow oy, (8) 
or® ip Ge r? Or TO 7 CP OOP}, = 


1 Ubersetzt von R. Eppler. 
) 2 F. Kégler u. A. Scheidig, Baugrund und Bauwerk, 8. 59 und 185, Berlin 1948. 
| 8 W. Fliigge, Bauingenieur 10 (1929), S. 221. 
4 §. Timoshenko, Theory of Plates and Shells, S. 263, New York 1940. 
5 K. Beyer, Die Statikim Stahlbetonbau, S. 671—672, Berlin 1948. 
6 F, Schleicher, Kreisplatten auf elastischer Unterlage, S. 13, Berlin 1926. 
7 Siehe S. Timoshenko, a. a. O. S. 414—415 und insbesondere K. Margugerre, Neuere Festigkeitspro- 
eme des Ingenieurs, S. 157—158 (Artikel K. Klotter), Berlin 1950. 
8 F. Télke, Besselsche und Hankelsche Zylinderfunktionen nullter bis dritter Ordnung vom Argument 


/i, S. 89 —90, Stuttgart 1936. fe: 
9 K. Girkmann, Flachentragwerke, S. 154—160, Wien 1948, und Proceedings of Symposia in Applied 
athematics, Bd. 3 S. 117—124 (Artikel K. O. Friedrichs), New York 1950. 
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fihren, wobei js die Poisson’sche Querkontraktionsziffer ist. An der Ubergangsstelle zum starre 


Mittelstiick muB ferner 
Ow eae ( 
Or r= r, wm i r= rT 


gelten. Endlich gibt noch die Gleichgewichtsbedingung zwischen dem aufgebrachten Moment 
und den Reaktionskraften der Unterlage 


To 2% 
Me | i cw(reos #)rd ddr 
0 6 


. r 
oder mit 0, = %— 
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Cormac 


Ms C [Ot 2 au Sy * 
Ma" vle) + [ V@ er dol = "he M (10 
Qi 
starres Mitelstiick 
Platte 70 
WE 
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X 
a 
0 
Abb, 2. Ergebnisse der Beispielrechnung ; 
w = 2,621-1077 7, cos 6 M, 
Abb. 1. Skizze des behandelten elastischen i yet hee nals Cee 
Problems, dr = — 1,68+ 10-3, cos  M. 


Die vier Randbedingungen (7) bis (10) fiihren auf ein System von vier linearen Gleichungen fii 
die A;. Zur Auflésung desselben mit Hilfe der Cramerschen Regel fiihren wir die folgenden 


Kiezancen ein: 


ae Ji Ned dJf, Tx 
ims Se 0 do o BO 
fJt , 1 Se 38—p a sy 
Ee ad 0 det 0? 2 BS 08 Ji (2) me ? 
= dJ¥, Shi (@) 
fi dg @ ist f 


fP= 4 Tule) + | Ther do, 


Analog werden fiir 1 = 2, 3,4 die GréBen F;, FF, f; und f* mit den in (5) bei den A; stehende 
Funktionen gebildet. Alle diese GréS8en kénnen mit Hilfe der Télkescnen Tafeln berechnet 
werden. Mit den Determinanten 


Her ny | OME EAE, Beis | 
Spe Fi FS FS Fi ae 0 FE FF Fi nee F¥ FS FF 0 a2 
ae k ss OU Seelam lobe Eva xe pated) 


fat han | 
Feet ft ft 


| uf Joehe } 
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wird dann die Lésung 


Oi 
A; = Oo . (13) 


Als Beispiel werde eine Platte mit w = as % = 2,2, rp =9m und r, = 3,6 m berechnet. 
Man findet 
A, = — 9,9226 M*, A, = — 3,8309 M*, A, = — 0,8918 M*, A, = 5,3999 M* 


woraus die in Abb. 2 dargestellten Ergebnisse folgen. 


Es sei noch bemerkt, daB die Herleitung der Gleichgewichtsbedingungen (10) nur dann in der 
obigen Form erfolgen darf, wenn die Platte so stark vorbelastet ist, daB sie sich nicht abheLen 
kann. Die entsprechende Lésung ist bei Schleicher! angegeben. 


(Eingegangen am 3], Juli 1952.) 
Anschrift des Verfassers: H. S. Gedizli, Cebeci Guvey Sk. 16 Ankara, Tirkei. 


| 1 F, Schleicher, a. a. O. S. 52. 
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Niherungsweise Berechnung und Stabilitatsverhalten 
von laminaren Grenzschichten mit Absaugung durch Einzelschlitze. 


Von W. Wuest. - 


1. Einfithrung. Die Grenzschichtabsaugung ist urspriinglich als wirksames Hilfsmittel zur 
Verhinderung der turbulenten Strémungsablésung und damit zur Steigerung des Hichstauftrie- 
bes von Tragfliigeln oder zur Erhéhung des Wirkungsgrades von Diffusoren angewandt worden. 
Spiter hat man erkannt, daB durch Absaugung der wandnahen Grenzschicht auch die Stabilitat 
der laminaren Grenzschicht wesentlich gesteigert werden kann. Eine Verlangerung der laminaren 
Laufstrecke ergibt aber eine Verringerung der Oberflachenreibung und damit eine Verkleinerung 
des Widerstandes von umstrémten Kérpern. Dabei wird das Stabilitatsverhalten in mehrfacher 
Hinsicht durch die Absaugung beeinfluBt. Einmal wird die Grenzschichtdicke iiberhaupt verkleinert, 
so daB jede mit einer Grenzschichtgrofe (z. B. Verdrangungsdicke) gebildete Reynoldssche Zahl 
des Umschlags zur Turbulenz erst weiter stromabwarts erreicht wird. Andererseits wird aber 
gerade durch die Absaugung der wandnachsten Schichten das Geschwindigkeitsprofil vélliger und 
damit die kritische Reynoldssche Zahl nach den vorliegenden Ergebnissen der Stabilitatstheorie 
fiir laminaren Grenzschichten wesentlich vergréBert. Dariiber hinaus ist bei Absaugegrenzschich- 
ten jenseits der Stabilitatsgrenze die Anfachung der Stérungen wesentlich geringer, so daB der 
tatsichliche Umschlag zur Turbulenz noch weiter hinausgeschoben wird. Bei sehr starker Ab- 
saugung durch Einzelschlitze kann auch der durch die Potentialstrémung vorgegebene Druck- 
verlauf beeinfluBt werden, indem ein zusatzliches periodisches Senkenfeld iiberlagert wird. Auck 
dieser EinfluB wirkt stabilisierend. Jedoch sind so starke Absaugungen nicht mehr lohnend. 


Bei der Grenzschichtabsaugung sind zwei Hauptarten zu unterscheiden, namlich die kon- 
tinuierliche oder flachenhafte Absaugung an halbdurchlassigen Wanden und die diskon- 
tinuierliche Absaugung an Einzelschlitzen, die auch in groBer Zahl aufeinander folgen 
kénnen. 


a) Kontinuierliche Absaugung. Die bisher vorliegenden Berechnungen beziehen sich 
fast durchweg auf kontinuierliche Absaugung langs der Wand. Die langsangestrémte ebene Platte 
mit homogener Absaugung, die nach einer geniigend grofen Anlaufstrecke zum sogenannter 
asymptotischen Absaugeprofil mit unverdnderlicher Grenzschichtdicke fiihrt, ist von 
H. Schlichting 1 untersucht worden. Der Anlaufvorgang selbst ist von R. Iglisch * exakt berechnet 
worden. Exakte Lésungen sind ferner von H. Schlichting und K. Bufmann ° fiir die ebene Stau- 
punktsstrémung mit homogener Absaugung sowie fir die langsangestrémte Platte mit einer 
Absaugegeschwindigkeit v, (x) ~1/|/* untersucht worden. Fir beliebige Druckverteilung und 
beliebige Verteilung der Absaugegeschwindigkeit hat H. Schlichting* ein Naherungsverfahren 
nach Art des Pohlhausenverfahrens bei undurchlassiger Wand entwickelt. Die Stabilitat der 
laminaren Grenzschichten mit kontinuierlicher homogener Absaugung ist verschiedentlich * ° mit 
dem bereits erwahnten Ergebnis untersucht worden. 


Versuche mit kontinuierlicher Absaugung sind in Langley Field an einem Tragfliigelmodell 
durchgefiihrt worden ?. Das Modell war mit einer 2,4 mm starken Schicht aus gesinterter Bronze 
itberzogen. Obwohl also das Modell nicht in aerodynamischem Sinne glatt war, konnte doch die 
laminare Grenzschicht iiber den ganzen Bereich aufrecht erhalten werden. Bei sehr hohen 
Reynoldsschen Zahlen (iiber 6 - 10%) war jedoch die notwendige Absaugemenge so gro8, da ins- 


gesamt kein Gewinn erzielt wurde. 


1 H. Schlichting, Luftfahrtforschung 19 (1942), 5.179 u. 293. 

2 R. Iglisch, Exakte Berechnung der laminaren Grenzschicht an der langsangestrémten ebenen Platte | 
mit homogener Absaugung, Schriften Dtsch. Ak. Luftfahrtforschung 8 B (1944), S. 1—51. 

3 H. Schlichting und K. BufSmann, Exakte Lésungen fiir die laminare Grenzschicht mit Absaugung 
und Ausblasen, Dtsch. Akad. Luftfabrtforschg. 7 B (1943), 5. 25. 

4 H. Schlichting, Ing.-Arch. 17 (1949), S. 199. 

5 K. BuBmann und H. Miinz, Die Stabilitat der laminaren Reibungsschicht mit Absaugung, Jahrb. 
Dtsch. Luftfahrtforschung 1942 I 36. 

6 J. Pretsch, Umschlagbeginn und Absaugung, Jahrb. Dtsch. Luftfahrtforschung 1942 I 1—17. 

7 A. E. v. Doenhoff and L. K. Loftin, Journ. aero. Sci. 16 (1949), 5. 729. 
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b) Diskontinuierliche Absaugung durch Einzelschlitze. Wesentlich eingehender 
ist das Verhalten von Grenzschichten mit einzelnen oder in groBer Zah] aufeinanderfolgenden 
Absaugeschlitzen experimentell untersucht worden. H. Holstein! untersuchte zwei Tragfligel- 
profile mit einer gré®eren Anzahl von hintereinander in bestimmten Abstanden angeordneten 
Absaugeschlitzen. Weitere Versuche wurden mit einer perforierten ebenen Platte durchgefihrt ?, 
die aber wegen technischer Mangel keinen Erfolg brachten. Ferner sind zwei Arbeiten von 
W. Pfenninger® zu erwabnen, in denen Tragfliigel mit mehreren hintereinander angeordneten 
Absaugeschlitzen untersucht wurden. Die zweite Arbeit von W. Pfenninger wurde an einem 
Tragfliigel mit 17% Dicke durchgefithrt, auf dessen Oberseite 24 und auf dessen Unterseite 16 
Schlitze angeordnet waren, die nach der Hinterkante des Fliigels dichter aufeinanderfolgten. 
Pfenninger schlieBt aus seinen Untersuchungen, daB die mit der Grenzschichtdicke gebildete 
kritische Reynoldssche Zah] durch die Absaugung nicht verandert wurde, so da also die Wirkung 
der Absaugung nur auf einer Verringerung der Grenzschichtdicke und auf der Senkenwirkung der 
Absaugeschlitze beruhte. Weitere Versuche mit einer groBen Anzahl von Schlitzen * haben ge- 
zeigt, da die Schlitzbreite nicht gréBer als die Grenzschichtdicke sein sollte. Auch spielt die 
sorgfaltige Bearbeitung und Glattung der Schlitzéffnungen eine wesentliche Rolle. 

Rechnerisch ist nur der Fall eines einzelnen Absaugeschlitzes ausgehend von der asympto- 
tischen Lisung naherungsweise behandelt worden 5, 

Im folgenden wird die Grenzschichtentwicklung an einer ebenen Wand naherungsweise unter- 
sucht. Dabei sollen an der Wand in bestimmten Abstanden Absaugeschlitze angeordnet sein. 
Die Absaugemenge wird als so klein vorausgesetzt, daB das Senkenfeld, das durch die Absauge- 
schlitze hervorgerufen wird, keinen merklichen EinfluB auf die Grenzschichtentwicklung hat. 
Der Druck wird also langs der Wand als konstant angenommen. Zur Berechnung wird ein Nahe- 
rungsverfahren angewandt, bei dem sowohl Impulssatz als auch Energiesatz angewandt werden. 
Die auf diese Weise gewonnenen Grenzschichtprofile sollen dann auf ihre Stabilitat untersucht 
werden. 


| 


2. Ableitung des Naherungsverfahrens. Bei den bekannten Naherungsverfahren zur Berech- 
nung laminarer Grenzschichten werden die Geschwindigkeitsprofile durch eine passend gewahlte 
inparametrige Kurvenschar angendhert. Die Kurvenparameter und die Grenzschichtdicke wer- 
en dann mit Hilfe der v. Kdérmdénschen Impulsgleichung und der aus der Prandilschen Grenz- 
schichtgleichung folgenden Randbedingung: 
| E du dU ou 
to(55), = Uae + * (a5), am 
erechnet. L. S. Leibenson * und K. Wieghard’ haben aufer der Impulsbedingung noch eine Ener- 
iebedingung aufgestellt und durch Wahl von zweiparametrigen Kurvenscharen die Rechnung 
erbessert, allerdings auch den Rechenaufwand vergréfert. A. Walz® hat gezeigt, dafi die 
eae auch dann gesteigert werden kann, wenn man beim einparametrigen Ver- 
fahren Impuls- und Energiesatz erfiillt, jedoch auf die Randbedingung (1) verzichtet. Im folgen- 
en wird ebenfalls ein einparametriger Ansatz fiir die Geschwindigkeitsverteilung benutzt. So- 
at die Geschwindigkeit der AuBenstrémung als konstant vorausgesetzt wird, ist es aber trotz- 
dem méglich, im schlitzfreien Bereich die Randbedingung (1) zusammen mit Impuls- und Energie- 
satz zu erfiillen und damit eine verhaltnismaBig gute Annaherung zu erzieclen. 
Die Geschwindigkeitsprofile werden also in folgender Form angesetzt: 


uj/U = F(n, K), 


1 FT, Holstein, Messungen zur Laminarhaltung der Grenzschicht an einem Tragfliigel, Ber. Lilienthal- 
Ges. Luftfahrtforsch. — H. Holstein, Messungen zur Laminarhaltung der Reibungsschicht durch Ab- 
saugung an einem Tragfliigel mit Profil NACA 0012/64, Forsch. Ber. Zentr. wiss. Berichtswesen Luft- 
fahrtforsch, Nr. 1654. i ; 
_ 2 H. Holstein und B. Pekarek, Messungen zur Laminarhaltung der Reibungsschicht durch Absaugung 
an einer perforierten ebenen Platte, Zentr. wiss. Ber. Wes. Luftfahrtforsch. UM 3222 (1945). 
_ 8 W. Pfenninger, Untersuchungen iiber Reibungsverminderungen an Tragfliigeln, insbesondere mit 
Hilfe von Grenzschichtabsaugung, Mitt. Nr. 13 Inst. Aerodyn. Ziirich 1946; W. Pfenninger, Journ. Aeron. 
Sci. 16 (1949), S. 227. 
4 Siehe FuBnote 7 von S. 90. as 

5 W. Wuest, Ing.-Arch. 17 (1949), 5. : : 
i 8 AN ae eae Form der Integralbedingung in der Theorie der Grenzschicht (russ.). 
Arbeiten ZAHI, H. 240 (1935). 

7 K. Wieghardt, Ing.-Arch. 16 (1948), 5. 231. 

8 A, Walz, Ing.-Arch. 16 (1948), S. 243. 
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wobei K (x) der Profilparameter und 7 = y/6 (x) ist (0 (x) ist ein spater noch definiertes Ma fiir 
die Grenzschichtdicke). Ferner fithren wir folgende GréBen ein: 


Verdrangungsdicke : ‘ 
6 =3 | (1G) dn = 9 (), (2) 
Impulsverlustdicke : : 
=8[ 5 (1 vr) a = 5 fa(K) (3) 
Energieverlustdicke : 


peat [e( (7) ) di, = of, (K), o 


Dimensionslose Dissipationsfunktion: 


co 


D(K) =2 i(° an dn, (5) 
on 
0 
Dimensionslose Wandschubspannung: 
du/U 
EK) ie ‘ 6) | 
Ali linear tal (6 
Mit diesen Abkiirzungen nimmt der auf Absaugegrenzschichten erweiterte Impuls- und Energie- 
satz folgende Form an: 
d (U?6 dU vU 
(A) 4, (0 —u) = — U6 RE, (7) 
d(U*6 Uu Of 
oAt oe — 09 (U8 — 0) = 2 7 E+v5D; ) (8) 


Uy, und vy sind dabei die Komponenten der Absaugegeschwindigkeit an der Wand. 

Wir wenden diese Bezichungen auf die Strémung an einer mit Absaugeschlitzen versehenen 
ebenen Platte an. Dabei wird eine die Potentialstrémung beeinflussende ,,Senkenwirkung“ der 
Absaugeschlitze vernachlassigt und demnach die AuBenstrémung als U = const angenommen. 
Ferner wollen wir voraussetzen, daB die Absaugung nur normal zur Wand stattfindet, so daB 
Ug = 0 (v, hat einen negativon Wert!). Damit geht die Impuls- und Energiegleichung iiber in 


dosage y ee 


PPMRD pe aE Re? (9) | 
dz vy 
te oS OSL (10) 


aa ABO U x 
oder mit am saa Z und > =& sowie (3) und (4) 


dZ “Peele 

A ety cla ae ee (11) 
dZ y Pic Samia ee 

ZESe tags Se — BZ =D. (12) 


Hieraus findet man durch Auflésung nach dZ/dé und dK/é 
df ah Ef, — Dis ely Z0 


dé Z (fe fi —fef}) : ae 
dK _ Ef, —Df, + (vo/U) Z(fs —f) (14) 
dé 2 (ffs —fef 5) ; 


Die Striche bezeichnen dabei Ableitungen nach dem Parameter K. Durch Division von (13) und 
(14) und Integration findet man schlieBlich 


K 
ie (EFL — DF) + (vf U) ZU —F) 
pp CS | (Ef, = DA) eda Co 


Ky 
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3. Lésungsansatz fiir Geschwindigkeitsprofil. Jeder einparametrige Ansatz fiir die Geschwin- 
digkeitsprofile in der Grenzschicht ist mehr oder weniger willkirlich. Man kann jedoch versuchen, 
sich dem erwarteten Verlauf der Geschwindigkeitsverteilung miglichst gut anzupassen. Im fol- 


_ genden beschaftigen wir uns vor allem mit der Entwicklung einer laminaren Grenzschicht an einer 
 langsangestrémten ebenen Platte, die an einzelnen Stellen durch Absaugeschlitze unterbrochen 
ist. Bei sehr groBem Abstand der Absaugeschlitze wird sich das Grenzschichtprofil jeweils weit- 


‘weitgehend dem Blasiusschen Plattenprofil annahern, andererseits wird sich bei einem Absauge- 
schlitz sehr grofer Breite der Grenzfall des ,,asymptotischen Absaugeprofils ergeben. Beide 
Grenzfalle sollen im einparametrigen Ansatz enthalten sein. Dabei wollen wir dem Vorgang von 
Schlichting} folgend, das Blasiusprofil in folgender Weise 


annahern: 


K=0 
reals fix 0S yn <3 | 
F (16) a x 
=7 
Das asymptotische Absaugeprofil hat andererseits die Form: a ty 
| Gol-e. (17) 
Wenn eine bereits ausgebildete Grenzschichtsstrémung in 
den Bereich eines Absaugeschlitzes kommt, mu8 man sich c x 
den Vorgang so vorstellen, da®B zunachst die wandnachsten ‘ 2 


Abb. 1. Einige Sonderfalle fiir die Verteilung 


Schichten beeinfluBt werden. Erst weiter stromabwartsgreift —4.; Andee cecschyind: Feit dues de Wee 


die Umgestaltung auch auf die a4uBeren Bereiche der Grenz- 
schicht iiber. Das gleiche wird auch der Fall sein, wenn der Absaugeschlitz wieder aufhdrt. 


' Aus diesem Grunde werden wir fiir den Schlitzbereich und den schlitzfreien Bereich nicht 
die gleichen einparametrigen Ansitze verwenden kénnen. Denn bei Verwendung des gleichen 


Ansatzes wiirden ja im einen Fall zuerst die 4uBeren Bereiche der Grenzschicht umgestaltet wer- 
den, wenn der Profilparameter riickwarts durchlaufen wird. 

a) Lésung im schlitzfreien Bereich. Wir wollen zunichst annehmen, daB® im Aus- 
gangszustand das asymptotische Absaugeprofil vorhanden ist und fiir «> 0 (vgl. Abb. la) die 
Absaugegeschwindigkeit Null wird. Die Geschwindigkeitsverteilung in der Grenzschicht wird 
sich dann in Richtung auf das Blasiussche Plattenprofil umgestalten, wobei offenbar zunachst die 
wandnachsten Schichten beeinfluBt werden. Dies legt folgenden einparametrigen Ansatz nahe: 


u 


uc 1-— e —F(K,y), (18) 
wobei 
Repl oil Ke oe) fa 0 = = Ol K, 
F (K,n) ( e sin ns | 8 
F (K,n) = — K e-0!VK fir 14> 3/K. 


Dieser Ansatz hat offenbar die Eigenschaft, daB die zweite Ableitung der Geschwindigkeit an der 
Wand unabhangig von K verschwindet. Fiir kleine Werte von K verschwindet F'(K, ») fiir 


n> 3 /K sehr schnell, so daB sich also Anderungen zunachst nur im wandnahen Bereich auswirken, 
Die Grenzschichtprofile, die sich mit diesen Ansitzen ergeben, sind in Abb. 2 fiir verschiedene 
Werte des Parameters K aufgetragen. In Tabelle 1 und 2 sind die GréBen fi, fy, f;, D, E zusammen 
mit den in (13) bis (15) gebrauchten Verbindungen fiir verschiedene Werte von K zusammen- 
gestellt. Es ist bemerkenswert, daB Ef,— Df, fir K = 0,968 verschwindet. Wenn man im An- 


- satz (19) das exakte Blasiusprofil statt der Naherung durch die Sinusfunktion einsetzen wiirde, 


verschwande Ef,— Df, erst fiir K = 1. Der Parameter K darf also bei unserer Naherung streng 


- genommen nur Werte zwischen 0 und 0,968 durchlaufen. 


Wegen v, = 0 nimmt das Integral (15) folgende numerisch leicht auswertbare Form an 
K 


Zz Ef; —Df; 20 
z= 0| { —Bpape ax}. (20) 
0 


1 H,. Schlichting, vgl. FuBnote 4 von 5. 90. 
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} 

| 
Die Lésung Z/Z, = 6/6) ist in Abb. 3 in Abhangigkeit von K aufgetragen. Der Wert Z) = 69 U/y | 
ist dabei der Grenzschichtdicke 6, am Anfang der schlitzfreien Strecke zugeordnet. Durch nume- | 
rische Integration der Differentialgleichung (14) kann sofort der Zusammenhang zwischen K und | 
—£ — Ux/y und damit auch zwischen Z und & ermittelt werden. Das Ergebnis dieser Rechnung | 
ist in Abb. 4 aufgetragen. Wenn K sich dem Wert K = 0,968 nahert, geht der Integrand gegen | 
| 


Unendlich. Wenn man die verschiedenen Funktionen in der Umgebung dieser singularen Stelle | 


j ki 
—— See eee Sa 
0 7 a 8 ¥ 0 G2 OF Go 08 10 
Abb. 2. Einparametrige Grenzschichtprofile nach Gl. (19) Abb. 3. Abhangigkeit der Grenzschichtdicke vom | 
in Abhangigkeit vom Profilparameter K. Profilparameter K im absaugefreien Bereich. 
linear annahert, erhalt man fiir groBe Tabelle 1. 


Werte von K (0,8 < K < 0,968) folgende Grenzschichtgréfen in Abhdngigkeit vom Profilparameter K. 


Naherungsformel : x f. f h a E 
Z\2 Zo 
(Gate O28) 22> ~(7,) wmnlsesk ne (ll) 0 1,00000 | 0,50000 | 0,83333 | 1,00000 | 1,00000 
0,1 1,00285 | 0,4965 0,82679 | 0,99093 | 0,84935 
b) Lésung im Schlitzbereich. On? 1,00806 | 0,4897 0,81289 | 0,97534 | 0,78695 
Wir wollen annehmen, da® die Grenz- _—-0,3'—*||:1,01481 | 0,4816 | 0.79527 | 0.95718 | 0.73907 
: B . 0,4 1,02280 | 0,4726 0.77381 | 0,93776 | 0,69870 
schicht zu Beginn des Absaugeschlitzes 0:5 1.03187 | 0.4630 0.75453 | 0.91781 | 0.66313 
4d 0.6 1,04189 | 0,4529 0,73251 | 0,89792 | 0,63098 
20 5 0.7 1,05279 | 0,4425 0,70847 | 0,87817 | 0,60142 | 
% 0,8 1,06450 | 0,4319 0,68605 | 0,85904 | 0.57389 | 
0,9 1,07696 | 0,4210 0,66179 | 0,84040 | 0,54805 — 
1,0 1,09014 | 0,40986 | 0,63662 | 0,82245 0,52360 
18 
Tabelle 2. 
16 K | ffs | fof | BDF, [DSi— EF 
SS 
0 03333339)/00 OZ3c0D0n nO 0 
44 0,1 .0,33029 | 0,02400 | 0,21024 | 0,0202 0,0035 
0,2 0,32319 | 0,0640 0,16208 | 0,0362 0,0095 
0,3 0,31367 | 0,0890 0,12679 | 0.0484 0,0155 
48 0,4 0,30118 | 0,1015 0,09747 | 0,0512 0,0195 
0,5 0,29147 | 0,1110 0,07499 | 0.0496 0,0224 
0,6 0,27407 | 0,1175 0,05553 | 0,0472 0,0248 
ry, 0,7 0,26078 | 0,1234 0,03749 | 0,0451 0,0266 
ANG G2 OY 06 48 40 0,8 0,25385 | 0,1280 0,02270 | 0,0426 0,0280 
: 5 4 0,9 0,24033 | 0,1331 0,00888 | 0,0409 0.0296 
Abb. 4. Verlauf der Grenzschichtdicke fiir den Fall 9 2 % ? ? ? 
Sick bh ik oe ee 1,0 | 0,22609 | 0,1370 | 0,00376 | 0,0406 | 0,0312 


entsprechend (19) folgende Form erreicht hat: 

= l—e1—F(K,»). 
Im Bereich des Absaugeschlitzes soll die konstante Absaugegeschwindigkeit — v, vorgesehen 
sein. Die einfachste Annahme fiir die Entwicklung des Geschwindigkeitsprofils im Bereich des 
Absaugeschlitzes wiirde darin bestehen, daB auch hier der Ansatz (19) gilt und der Parameter K 
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wieder von K, nach Null lauft. Dieser Ansatz wiirde aber nicht beriicksichtigen, da zunichst die 
wandniachsten Schichten von der Umgestaltung erfaBt werden. Besser ist daher folgender Ansatz: 


u 


Vas 1! en F (K,, 7) + F(K, n) (22) 


i 
Hierbei ist die Funktion F (K, 7) wie in (19) definiert. Der Parameter K durchlauft die Werte von 
(0 bis K,. Die Randbedingung (1) wird durch diesen Ansatz nicht erfiillt. Im Schlitzbereich ist 
also das hier angewandte Verfahren identisch mit dem einparametrigen Naherungsverfahren von 
A, Walz1, Fiir die weitere Rechnung ist dies aber nicht schwerwiegend, da es vor allem darauf 
-ankommt, im schlitzfreien Bereich eine méglichst gute Annaherung zu erzielen. 


O48 940 


aye 915 


guy O00 


O42 95 G65 


0 G2 GF 46 GE 40 0 G2 OY 46 G8 40 
Abb. 5. f, und df,/dK im Schlitzbereich. Abb. 6. f,; und df,/dK im Schlitzbereich. 


Die GréBen f,, df,/dK, f, und df,/dK sind in Abhangigkeit der beiden Parameter K und K, 
 graphisch bestimmt worden und in Abb. 5 und 6 aufgetragen. Fiir D und E erhalt man 


1 ee aS 
WK al Kav) 


D = Deto:0,13587-K? a 4( 


2a ‘K a K K,3? _ 2% n/6 VK, oe 3VKIK, 4+/K 
mage Pk one eee 8 a ere RI alae (23) 
2m gp Wo VRePVRIK 4 /K, 20K ajo o°!® +/K 
3 J K, + Kx?/36 re a 
E = 1— 0,47640 (/K, —/K) . (24) 


Wir wollen als Beispiel den in Abb. 1b dargestellten Fall betrachten, bei dem im Ausgangszustand 
das Blasiusprofil (K, = 1 bzw. genauer K, = 0,968) vorhanden ist und dann ein Bereich homo- 
gener Absaugung beginnt. Wir wollen als Berechnungsbeispiel Z)/Z. = 5 annehmen, d.h. die 
Ausgangsgrenzschichtdicke betrage das Fiinffache der asymptotisch bei unendlich langem Schlitz 
erreichten Grenzschichtdicke. Abb. 7 zeigt den errechneten Verlauf von Z/Z,. Diese Funktion 
sinkt vom Anfangswert 5 mit zunehmendem K schlieBlich fiir K = K, auf den Wert 3,39 ab. 
Damit ist aber noch nicht der Endzustand erreicht, da ja Z/Z.. asymptotisch gegen den Wert | 
gehen mu8. Man muB also fiir den weiteren Verlauf der Grenzschichtentwicklung passende An- 
nahmen machen. Dabei kann man insofern etwas groBziigig vorgehen, als man eine genaue Dar- 
stellung der Strémungsvorgange in unmittelbarer Umgebung des Schlitzes von einer Grenz- 
schichttheorie ohnehin nicht erwarten kann, weil hier schon die Grundvoraussetzungen nur un- 
geniigend erfiillt sind. Rechnerisch am einfachsten ist die Annahme, da8 nach Erreichen von 
K = K, die Profilform unverandert bleibt und nur die Grenzschichtdicke standig weiter ab- 


1 Siehe FuBnote 8 von S. 91. 
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nimmt. Wir denken uns also die Grenzschichtentwicklung im Schlitzbereich in zwei Schritte zer- 
legt. Der erste Schritt vollzieht sich nach Gl. (15), wobei K von 0 bis K, ansteigt. Beim zweiten 
Schritt spielt sich dann die Grenzschichtdicke bei festgehaltener Profilform (K = K,) auf den 
geforderten Wert ein. Bei alleiniger Anwendung des Impulssatzes (11) ergibt sich dann mit K, f, 


und F£ konstant 
1 d(U*/Zu%) — fo\*/, 9 Z 
4 dé i (3) ( Za) 


pend (O) GF pm Bae, 09) 


Vo 


mit dem Integral 


In Abb. 7 kommt der zweite Schritt nur als senkrechte Linie CA mit K = K, zum Ausdruck. 


4, Periodische Absaugegrenzschichten. Bei sehr vielen, in gleichem Abstand aufeinander- 
folgenden Absaugeschlitzen (vgl. Abb. 1c) wird sich nach einer gewissen Anlaufstrecke eine perio- 
disch von Schlitz zu Schlitz wiederkehrende Geschwindigkeitsverteilung einstellen. Der Verlauf 
der Grenzschichtdicke 6 (bzw. des Verhaltnisses 6/6, = Z/Z,.) istin Abb. 8 in Abhangigkeit vom 
Profilparameter aufgetragen, und zwar im 
schlitzfreien Bereich (Kurvenstiick c) itber K, 
im Bereich des Schlitzes dagegen itber K,—K, 


so da® sich eine geschlossene Kurve ergibt. 


a Abb. 8 bezieht sich auf den Fall K, = 0,8 und } 
Z3/Z ie 
3 
Z 
52 B 
a 
4 
if 
k 
Ky-k A Ky-K 
0 FPL QT hy Po GPa e 49 a2 a4 06 08 
Abb. 7. Verlauf der Grenzschichtdicke fiir den Fall Abb. 8. Verlauf der Grenzschichtdicke fiir periodisch 
nach Abb. 1b. aufeinanderfolgende Absaugeschlitze (Fall nach Abb. 1c). 


Es bereitet nach dem vorhergehenden keine Schwierigkeiten, die den drei Teilabschnitten 
a, b, c zugeordneten Langen x,, x), x, zu berechnen und damit auch den Schlitzabstand a — us 
-++ x, -+- x, und die Schlitzbreite s = x, + x, zu bestimmen, die einer geschlossenen Kurve im | 
Diagramm nach Abb. 8 zugeordnet sind. Im vorliegenden Beispiel erhalt man 
2 
U4 _ 4,88 es und 2 = 0,71 eae 
v Vo v Vo 
Damit betragt in diesem Fall das Schlitzverhaltnis A = s/a = 0,145. | 
Es erweist sich als zweckmaBig, statt der Absaugegeschwindigkeit Vy am Schlitz die mittlere 
Absaugegeschwindigkeit v, = A v, einzufiihren. Wir definieren ferner den »Absaugeparameter‘ — 


© durch 
Ua/v,\ v a \2 
ames (a) = walso) (26) 


Fur das Rechenbeispiel ergibt sich 
O = 4,88 - 0,1452 = 0,103 . 


Dabei ist 6.., die Grenzschichtdicke, die sich asymptotisch bei homogener Absaugung mit der 
mittleren Absaugegeschwindigkeit v, einstellen wiirde. 

Durch Vorgabe von © und 4 ist eine geschlossene Kurve im Diagramm nach Abb. 8 im all- 
gemeinen bestimmt. Damit ist auch der GréBtwert des Formparameters K gegeben, der im_ 
Punkt B erreicht wird. Dieser Wert Ky. = , ist aber entscheidend fiir die Stabilitat des Ge- 
schwindigkeitsprofils. Es zeigt sich nun, daB das Schlitzverhaltnis, solange es itberhaupt klein ist, 


fast gar keinen Kinflu8 hat. In Abb. 9 ist fiir K, = 0,8 der Absaugeparameter @ iiber dem Schlitz- 


= 
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verhaltnis A aufgetragen. Fitr Werte J < 0,15 ist O praktisch konstant. Diese Aussage bezieht 
sich natiirlich nur auf die Grenzschichtstrémung. Auf die Strémung im Schlitz selbst wird 
natiirlich die Schlitzbreite einen wesentlichen Einflu8 haben. Doch liegt dies auBerhalb der vor- 
liegenden Betrachtung. In Abb. 10 ist der sich aus der Rechnung ergebende Verlauf des Ab- 
saugeparameters @ fi kleine Schlitzverhiltnisse J in Abhingigkeit von K, aufgetragen. Nach- 
dem dieser Zusammenhang hergestellt ist, besteht die niichste Aufgabe darin, die kritische 
Reynoldssche Zahl der Grenzschichtprofile zu berechnen. Erwartungsgema8 nimmt die kritische 
Reynoldssche Zahl mit zunehmendem Profilparameter K monoton ab, so daB also am Ende der 
schlitzfreien Strecke das instabilste Profil vorhanden ist. Es besteht also die Aufgabe, die Funk- 
tion Rez,;: (Kmax) zu bestimmen. 


___ 5. Die Stabilitat der ebenen Laminarstrémung. Mit Hilfe der von Tollmien! entwickelten 
Theorie der kleinen Schwingungen soll die Stabilitat der Laminarstrémung hinter einem Ab- 
saugeschlitz untersucht werden. Diese Methode ist schon mehrfach ® ° sehr eingehend behandelt 
worden, so dafs es erlaubt ist, auf eine ein- K 

gehende Darstellung und Begriindung an i 


a ; 98 + : ; 
dieser Stelle zu verzichten. Bemerkens- | 
werterweise ist fiir die Stabilitatsunter- 
ssuchung nur die Form der értlichen Ge- 1} 
schwindigkeitsverteilung maBgebend, wah- 
rend der EinfluB der x-Abhangigkeit der 
{ 08 l 
@ 
65 ib 
ge OF i 
G3 
b2 a2 ial 
Of | | i a 
| 
A 0- Fle) 
eo 
0 G2 U4 06 G007 2 rr a aa 2 UMA CMEENOT 
Abb. 9. Abhangigkeit des Absaugeparameters @ Abb. 10. Abhangigkeit des Absaugeparameters © vom Maximalwert K, 
vom Schlitzverhiltnis (fir K, = 0,8). des Profilparameters fiir kleine Schlitzverhaltnisse (A < 0,15). 


‘Stérungsamplitude vernachlissigbar gering ist und die x-Abhangigkeit der Grundstrémung nur 
mittelbar iiber einen etwaigen Druckgradienten EinfluB gewinnt. In unserem Falle kénnen aller- 
dings in unmittelbarer Umgebung der Absaugeschlitze Schwierigkeiten auftreten, weil sich dort 
‘die laminare Grundstrémung sehr rasch andert. Jedoch kénnen wir die unmittelbare Umgebung 
der Absaugeschlitze ausschlieBen, weil sich erwartungsgem48 die geringste Stabilitat am Ende 
ider schlitzfreien Strecke ergibt. 

| Wie die Versuche von Schubauer und Skramstad* gezeigt haben, ist Voraussetzung fiir die 
Giiltigkeit der Stabilitatstheorie, daB die Anfangsstérungen sehr klein sind (Turbulenzgrad 
< 0,2%). Bei vielen Versuchen mit hintereinander angeordneten Absaugeschlitzen ist diese Be- 
dingung wahrscheinlich nicht erfiillt gewesen. Die Stabilititstheorie ist bisher — von Sonder- 
fallen abgesehen — nur auf Geschwindigkeitsverteilungen angewandt worden, die sich ganz oder 
tiickweise durch Binome ausdriicken lassen. Die im folgenden mitgeteilten Formeln beziehen 
sich dagegen auf Geschwindigkeitsverteilungen, die durch beliebige Potenzreihen dargestellt 
werden kénnen. Man hat dabei aber im allgemeinen die Schwierigkeit, daB diese Potenzreihen 


1 W. Tollmien, Nachrichten Ges. Wiss. Géttingen, Math.-Phys. KI]. 1929, S, 21. 
2 C. C. Lin, Quart. Appl. Math. 3 (1945), S. 117. 

3 H. Holstein, Z. angew. Math. Mech. 30 (1950), S. 25. 

4G. B. Schubauer and H. K. Skramstad, Journ. Aeron. Sci. 14 (1947), 5, 69. 
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wegen des eigentiimlichen asymptotischen Verhaltens der Grenzschichtgleichung im wandfernen 
Bereich immer schlechter konvergieren. Um praktischen Bediirfnissen zu dienen, wird im folgen- 
den ein Weg gezeigt, bei dem es geniigt, wenn die Geschwindigkeitsverteilung in Wandniahe analy- 
tisch bekannt und durch eine Potenzreihe ausdriickbar ist, wahrend im wandfernen Bereich die 
Geschwindigkeitsverteilung nur irgendwie graphisch oder numerisch vorgegeben sein mu, ohne} 
durch eine analytische Beziehung ausdriickbar zu sein. Zur Darstellung dieser Erganzungen der 
Stabilitatstheorie mégen im folgenden die wichtigsten Beziehungen der Tollmienschen Theorie} 
in aller Kiirze zusammengestellt werden. 

Der laminaren Grundstrémung U(y) seien zweidimensionale Stérungen tiberlagert, wobei die } 
Stromfunktion einer Partialschwingung lautet: 


F(x, 9,1) =—igpgen res: 
Hierbei ist 7=y/d und a=ay 6 (a9 Kreisfrequenz der Stérung) reell, wahrend c = c, +1 ¢; im all- 
gemeinen komplex ist, wobei ¢, die Phasengeschwindigkeit der Stérung und c; ihre Anfachung | 
bezeichnet. Da wir uns auf die Berechnung der Stabilitatsgrenze, also auf neutrale Schwingungen 
beschranken, ist c; = 0 und c= c,. Die Amplitude p(n) der Stérung gehorcht dabei der Tollmien-; 
schen Stérungsgleichung | 


u—c te 2 


oP — 17 =a ieee eeeaangae (27)) 


mit der Reynoldsschen Zahl Re = Ud/v. Die allgemeine Liésung der Differentialgleichung (27)) 
kann in folgender Form angesetzt werden: 


P =O, Py + ly Po +3 V3 + €4 Pa; (28) 


hierbei sind g, und ~, mit dem Wandabstand ,langsamveranderliche“ Lésungen, 3 und pq ,.schnel! | 
veranderliche‘‘ Lisungen. Dabei wachst aber nach Tollmien :, iiber alle Grenzen, so daB man mit 
Riicksicht auf die Randbedingungen c, = 0 setzen mu. Die iibrigen drei Integrationskonstanten } 
ergeben sich aus den Randbedingungen p= yp’ = 0 an der Wand (ny = 0) und yg’ =a fiir groBen h 
Wandabstand. Die letzte Randbedingung folgt daraus, daf fiir groBen Wandabstand u” gegen) 
Null geht und auch die Reibungswirkung vernachlassigbar wird, so daB dort nach (28) 


/ 


O10" 0 ee pe ne=te aa) 
wird. 

Die langsamveranderlichen Lésungen gy, und @, entsprechen im wesentlichen den Lésungen |} 
der reibungslosen Stérungsgleichung, die man erhalt, wenn a Re so groB wird, daB die rechte Seite} 
von Gleichung (27) vernachlassigt werden kann. Die reibungslose Differentialgleichung hat einen } 
Pol erster Ordnung an der Stelle u = c, der wir den ,,kritischen‘‘ Punkt 7 = y* zuordnen. Physi-| 
kalisch ist diese Stelle dadurch ausgezeichnet, da die Fliissigkeit stets in der gleichen Stérungs- 
phase schwingt. 

Wir fithren als neue Veranderliche 7, = (7 — 7*)/m ein,wobei m ein noch nicht naher bestimm- 
ter Mafstabsfaktor ist. Das laminare Grenzschichtprofil sei als Potenzreihe gegeben und in 
folgender Weise entwickelt : 


u—c ~ hee oe 
i=l 
Die Lésungen ¢y, und g, sind dann gegeben durch | 
Pr = My 2s On Ns (30) 
P= D>) bam EEO (In 4, — 1) . (31) 
- | 


Das imaginare Glied iz ist nur fiir negative Werte von 7, vorhanden, wahrend es fiir positive 7; 
verschwindet. Bei Beriicksichtigung der Fliissigkeitsreibung vollzieht sich dieser Phasensprung 
beim Durchgang durch die ,,kritische*‘ Schicht auf einer endlichen Wegstrecke. Ebenso ver- 
schwindet auch die logarithmische Singularitat von y, bei Beriicksichtigung der Fliissigkeits- 
reibung. Die Koeffizienten a, und b, lauten fiir eine beliebige, entsprechend (29) als Reihenent- 
wicklung gegebene Geschwindigkeitsverteilung : 
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le, = | Geis Ares ole. oy al Be Ot By 
4 4 apie © 24 6, | 18 B,’ 


ee | Bs 1 /Bs\? _, 1 Bs 1 
et 120 8, ane Y 60 Bot pat 


ae WR Gee ot, Baty a he By abeaks 
ae 2400(2,) — 5400 B? ' 340 B,|°* + foa00 7,” (30a) 


1 1 (Be 11 BE Bs 1 /Bs 1 Bs Ba 13 Bs] 12 
4 ~ 5040 a Bl FIN He r 27600 Ba seal cae 2520 p3 Ieee 
Bs 


_ 31 (Ba Bs] _1 
| | 0) F apie at + aap a8 
eee ee 


ie Be) 1 Be Bs 1 Bp, 1 By 
w= [a(n] oo tin] et 


1 


oh 


a — 1 /B2\* I BE Bs nS Bs\? 1 Bs Bs 1 B; j 61 / Be\? 1 /Bs ms 

i ma 2 py 7 eae a4 pe | 72 | | male | rz (f°) a tare (31a) 
Zz Bi 1 63h, 1 Pa62 1 BEB. 1 BeBa BB, | 1 Be 

Ps ee cota + a pi 72 ~«O#BSSCOB:COBR + 96 pi 1440 6? + 480 a 


109 /Bs 7 Be Bs 1 By 2 

te li are pe) das B2 + oo Ho race 

\Hierbei ist 1; = am. In den Reihenkoeffizienten (30a) und (31a) sind auch die von Pretsch fiir 
binomische Geschwindigkeitsverteilungen angegebenen Koeffizienten als Sonderfalle enthalten. 
~ Die Reihenentwicklungen (30) und (31) sind nur im wandnahen Gebiet brauchbar. Wegen der 
dritten Randbedingung braucht man aber auch eine Darstellung der Lésungen y, und g, am 
AuBenrand der Grenzschicht. Statt nun die Funktionen gy, und q, nach der Koordinate 7, zu 
sntwickeln, wobei die Entwicklungskoeffizienten selbst wieder von «0, = a) 6m (do Kreisfrequenz 
der Stérung) abhangig sind, kann man die Funktionen auch nach «, entwickeln, wobei die Koeffi- 
ienten umgekehrt von 7, abhangen!: 


Pa Vito + 0? gy, + of Pip + 08 Gig +°°°, 


(32 
Pz + Poo + 0? Yor + OF Yoo + 08 Mog +°->. | ) 
Die v,, (v = 1,2) gewinnt man iterativ aus 
Uy an ne pe 
7 Mn — P= Pans (33) 


Pir n = 0 verschwindet die rechte Seite und man erhAlt als Anfangsglieder der Entwicklungen (32) 


1 
Pi ayy Vibes (34) 


| 
we "G 
| 


Un 


u—c [ dn 
Po = — Bi U wi 3yue (39) 


3] 


Die Konstante f, ergibt sich durch Vergleich mit den Reihenentwicklungen (30) und (31), in denen 

man a, = 0 gesetzt hat. Die Integrationsgrenze 7 in (35) ist noch nicht naher bestimmt. Wir 

c<énnen uns aber in folgender Weise davon befreien: Wir berechnen q 9 aus der Reihenentwick- 

ung (31) durch Nullsetzen von «, und zwar fiir einen Wert 7,, der noch im Konvergenzhbereich der 
eihenentwicklung (31) liegt. Aus rechnerischen Griinden ist es vielfach zweckmaBig, den Punkt 

4 = 1 zu wahlen, so daB 

P20 (1) = 1 + beg + O39 + O49 + °° (36) 


1 Diese Méglichkeit ist bereits in der Tollmienschen Originalarbeit diskutiert worden, ohne daf aber 
tin rechentechnischer Nutzen daraus gezogen wurde. 
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wird. Andererseits ist dann aber 
u—c dy 
P20 (1) ar — fi U inh (wo ou ? 


so da (35) umgeformt werden kann in 


1 
0 = bi |G (Ga), | are): (35a 
7 
Fiir groBe Werte von 7, verhalten sich 99 und Yoo asymptotisch wie 
Pi0~- =e. (37 
By ~ 38 
Poo T= 4 1) (38) 


Die Konstante a berechnet sich dabei zu 


a= go()(—) g(a) +4— i (=e) Yen (39) 


Das Integral kann graphisch ermittelt werden. Andererseits geht aber die reibungslose Stérungs-) 
gleichung fiir grofen Wandabstand iiber in y’’—a?p=0 mit den Lésungen y =Coja7 bzw. 
Gina, und durch Vergleich mit (37) und (38) erhalt man schlieBlich als Lésungen fiir groBen 
Wandabstand 


=o | 
Vice ae Co} om » (40)) 
al 

Pz ™~ bs c Co] a1 — a Sin a, m| . (41) 

1—4/U Oy 

Fiir die spatere Rechnung ist folgende Beziehung wichtig, die jetzt leicht berechnet werden kann: 
ate Pot %1 Pe ae By 2 1 \ 

i) Pit Py ( aie! 7 (« a : (42) 


Bei den langsam veranderlichen Lésungen ¢y, und @, hat die Reibung nur auf @, einen KinfluB, une 
zwar nur in der Umgebung der ,,kritischen‘‘ Schicht, wo die Phasengeschwindigkeit der Stérung 
mit der Strémungsgeschwindigkeit in der Grenzschicht iibereinstimmt. Die schnell veranderliche 
Lésung y, ist dagegen wesentlich reibungsbedingt und auf die Wandnahe beschrankt, wahrend 
sie mit zunehmendem Abstand von der Wand schnell abklingt. Fiihrt man die neue Verdnderliche’ 


Va : 
Uae (43) 
ein mit | 

—1/3 

€ =| Re Cae / : 

a) Vi 
wobei wieder 
Re ALO ao=o,6=—, 
v dy m m 


so erhalt man aus der Stérungsgleichung (32) unter Vernachlassigung der Glieder mit dem Fak- 
tor ¢ die Differentialgleichung 


sti 1 
ey Ok (44) 
2 


Die gesuchte Lisung @, ist dabei diejenige schnell veranderliche Lésung dieser Differential- 
gleichung, die fiir groBen Wandabstand verschwindet. Die Lésung q, ist bemerkenswerterweise 
von der speziellen Form des Geschwindigkeitsprofils ganz unabhangig und kann daher universell 
berechnet werden. Die bereits genannten drei Randbedingungen fihren auf folgende Determi- 
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nantenbedingung: 
Piw Paw P3w 
Y u dys w dig 
Piw P20 “dy dy, | =9- (45) 
(Pi +O Pr) (2 +41 P2)a 0 


Die Striche bezeichnen dabei Ableitungen nach 7,, der Index w die Wand und a den AuSenrand 
ler Grenzschicht. y,ist am AuBenrand der Grenzschicht schon abgeklungen. Die Auflésung der 
Determinante fiihrt auf folgende Gleichung, wenn man noch L(q,) aus (42) einfithrt : 


dps wv \K Pow —L Pw 
w 7 Tae 46 
Ps / djy E PZw — L Pw ( ) 

Nun ist aber 
dy w 
Paw | ai = —D (new) » 
IEF) J(EF) 
OF S ag 
Nee - |A=G04 
Nanas , 


awe = SS 
GI Ga G3 OY GS G6 G7 0 G7 G2 G3 oF GS GO 
Abb. 11. Graphische Bestimmung der Liésungen von Gl. (54) Abb. 12. Graphische Bestimmung der Lésungen von Gl. (54) 
fiir K — 0;2- fir K = 0,04. 


wo Dy ~) die in tabellierter Form vorliegende komplexe ,, Tietjenssche Funktion“ ist. Dividiert 


nan beide Seiten von Gleichung (46) noch durch 7,,, = —7*/e, so erhalt man schlieBlich mit den 
Abkiirzungen 

F’(taw) = — D (Naw)/Naw » (47) 

E(o, 7*) = VK g.w—L ow (48) 


| n* G20 —L pw 
die Bedingungsgleichung 
F (tw) = E (a, 9*) ° (49) 
Zur Ermittlung der Indifferenzkurve tragt man nun den imaginaren Teil von F iiber dem Realteil 
nuf (ausgezogene Kurve in Abb. 1] und 12). Hierzu zeichnet man zu gegebenem 7* und damit 
su gegebenem c/U die Funktion E(a,). Die Schnittpunkte beider Kurven erfiillen die Bedingung 
49). Man findet demnach bei gegebenem 7* die zu den Durchgangspunkten gehérenden Werte a, 
ind 7,.. Als Endergebnis sind schlieSlich die auf die Verdrangungsdicke 6, bezogenen Werte 
ye = a, 6,in Abhangigkeit von der Re-Zah] Re* = U6,/e gesucht. 

) 6. Anwendung auf Absaugegrenzschichten. Wir beniitzen die im vorhergehenden Abschnitt 
ntwickelte Rechenmethode zur Berechnung der Stabilitat der frither aufgestellten Grenzschicht- 


alae —F(K, 9"), (50) 
| 
and die Koeffizienten B; der Reihenentwicklung (34) lauten bei cinem Mafistabfaktor m = 1 
| cos — dibs (« ungerade) , 
; ats Seis rd 6 \K 
(— 1)! By = en" — KI? et VE : 
sik (i gerade) . 

6 VK 


| 
| 
| 
) 
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In Abb. 11 und 12 sind fiir die Werte K = 0,2 und K = 0,04 die komplexen Funktionen EF (taw} 
und E(a,,7*) aufgetragen. Im einen Fall ist also 72 =— n*/e, im anderen Fall ist auf jeder ein4 
zelnen Kurve a, der laufende Parameter. Die Schnittpunkte der F-Kurve mit der E-Kurvenschani 
jiefern zusammengehbrige Werte 1, Now: n*. Man bestimmt daraus 


a* = ot 0, = 0, fi (K), (52) 
_ U6; /_ tew\? 1 53) 
Reet a OO — Oe) afi (X) (53) 


Hierbei ist f,(K) aus Tabelle 1 zu entnehmen. Auf diese Weise sind die Kurven in Abb. 13 er- 
halten worden, Zusiitzlich wurde noch die Kurve fiir das Blasiusprofil (K = 0,968) nach einer 
fritheren Rechnung von Tollmien! und die davon etwas abweichende Kurve nach C. C. Lin? ein- 
getragen. Ferner ist auch die Kurve fiir das asymptotische Absaugeprofil (K = 0) nacht 
J.Pretsch® eingetragen worden. Das Ergebnis von Pretsch mit einer Stabilitatsgrenze vont 
Re =5,52 104 ist neuerdings angezweifelt worden’. Eine vom Verfasser durchgefiihrte Neu- 
berechnung ergab innerhalb der Rechen-} 
aio genauigkeit (die etwa --5 % betragen diirfte) 
ff “K R fast wieder die friheren Werte (genauert 
ee 5,8 104). Es sei ferner noch darauf hin-) 
, gewiesen, daB C. C. Lin® eine etwas proble- 
matische Naherungsformel zur Berechnung 
der Stabilitatsgrenze angegeben hat, die| 
wegen ihrer Hinfachheit in amerikanischen| 
Arbeiten gern benutzt wird. Bei der An-| 
wendung auf das vorliegende Problem liefert 


4 


Blasiusprotil 
\ nach Tolimien 
nach CC Lin 


02 die Naherungsformel folgende Werte im Ver- 
gleich zu der hier durchgefiithrten Stabilitats- 
rechnung : 

Q7 Naberungs- 

Stabilitatsrechnung formel 
nach C, C. Lin 
yg &=0 5,52 104( Pretsch) 3,9) 108 
v 4,1 104(Freeman) 
4-10 10" 20" 10° 2:10" Kee 0,04 0,68 104 Nees 104 
Abb. 13. Tollmiensches Stabilitaitsdiagramm fiir verschiedene K=0,2 0,178 104 A 0,406 104 
Werte von K. Ik SS 0,968 420 ( Tollmien) 450 


Nach der Linschen Formel geht also der Abfall der kritischen Reynoldsschen Zahl mit wachsen- 
dem Parameter K viel langsamer vor sich als nach der hier durchgefiihrten Stabilitatsrechnung. 
Wir werden aber fiir die folgenden Betrachtungen die Ergebnisse der Stabilitatsrechnung zu- 
erunde legen. Die Auftragung der obigen Werte im Diagramm (Abb. 14) 1a4Bt den raschen Abfali- 
der kritischen Reynoldsschen Zahl mit wachsendem Parameter K erkennen. Dies liegt offenbar 
daran, da fiir die Stabilitat der wandnahe Verlauf des Grenzschichtprofils ma®gebend ist. Dieser 
Bereich wird aber nach dem Ende des Absaugeschlitzes sehr rasch wieder umgeformt. 
Verbinden wir die in Abb. 10 dargestellte Abhangigkeit des Absaugeparameters 


aha? 2 ets (iF 6 

0=—*(8) = 0a, ( 3) 
von K,,4x = K, mit der in Abb. 14 dargestellten Abhangigkeit der Stabilitaétsgrenze von K, so 
gewinnen wir daraus die in Abb. 15 dargestellte Abhangigkeit der kritischen Re-Zahl, Re = U6,/» 
von (6,,/a) (6,,/6;). Dabei ist 6/0, nicht wesentlich von 1 verschieden. Gleichzeitig ist in diesem 


Kurvenblatt auch der Mengenbeiwert cg = — v,/U aufgetragen. Bei der Berechnung des Mengen- 
beiwertes wurde naherungsweise 6, /6, = 1 gesetzt. 


1 Siehe FuBnote 1 von S. 97. 

2 Siehe FuBnote 2 von S. 97. 

3 Siehe FuBnote 6 von S. 90. 

* P. Chiarulli and J. C. Freeman, Stability of the Boundary Layer. Hdqtrs Air Nat. Comm, Dayton. 
Techn. Rep. No. F-TR 1197-1A Aug. 1948, S. 1—104. 

5 Siehe FuBnote 2 von S. 97. 
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Wenn man bedenkt, daB beispielsweise bei einer Strémungsgeschwindigkeit von 100 m/s die 
aulassige Grenzschichtdicke etwa von der GréBenordnung 1/10 mm ist, erkennt man, da8 nur 
verhaltnismafig kleine Werte des Verhaltnisses technische Bedeutung haben. Aus Abb. 15 er- 
kennt man dann aber, daB eine wesentliche Erhéhung der kritischen Reynoldsschen Zahl nicht zu 
erreichen ist. 

Man muB also aus den vorstehenden Uberlegungen die SchluBfolgerung ziehen, daB® es durch 
Anwendung aufeinanderfolgender Absaugeschlitze nicht méglich ist, die kritische Reynoldssche 


Rez wesentlich zu steigern. Fir verhaltnismafig kleine 
08 Geschwindigkeiten mit groBer Grenzschichtdicke 
8 ie se wird man dagegen eher groBe Werte von 6,,/a 
6 | erreichen kénnen und daher auch die kritische Re- 
: | Zahl auf den mehrfachen Betrag steigern kénnen. 
a Die vorstehenden Uberlegungen sind in mehr- 
facher Hinsicht unvollkommen. Die Grenzschicht- 
v4 ie 
1 Re pnit y=-ufv 
7500 G003 
pa leeeties. [ie a | 
& td 
=! = at 
: | 7000 002 iN 
N 
8 
Ss 
S 
2 a & 
SS 
500 4007 
70° - | 
8 = | 
6 =i 
doo 
Y af * 
0 02 G4 06 Ge 10 0 GOS 470 GIS 
Abb, 15. Abhangigkeit der kritischen Reynoldsschen Zahl (Stabilitits- 
Abb. 14. Abhangigkeit der Stabilitatsgrenze vom grenze) von dem Verhiltnis von mittlerer Grenzschichtdicke zum 
Profilparameter K. Schlitzabstand. Die gestrichelte Kurve gibt den Mengenbeiwert 


y= v,/U an. 


iprofile sind nicht exakt berechnet, sondern mehr oder weniger willkirlich als einparametrige 
‘Kurvenscharen angenommen. Das periodische Senkenfeld der Absaugeschlitze ist nicht be- 
‘riicksichtigt worden. Verfasser beabsichtigt, in einer weiteren Mitteilung eine genauere Theorie 
‘der periodischen Absaugegrenzschichten zu geben und damit die genannten Mingel zu beseitigen. 
|Aber auch beziiglich der Anwendung der Tollmienschen Stabilitatstheorie auf periodische Ab- 
jsaugegrenzschichten kénnen Einwande erhoben werden. Man kann sich beispielsweise vorstellen, 
daB Stérschwingungen auf einem Teil der Wegstrecke angefacht werden, auf einem weiteren Teil 
'dagegen gedampft, ohne daB im ganzen gesehen ein Umschlag zur Turbulenz erfolgt. Die Stabili- 
tatsgrenze diirfte also bei periodischen Absaugegrenzschichten héher liegen, als sich durch formale 
Anwendung der Tollmienschen Theorie ergibt. Um einen tieferen Einblick zu gewinnen, wird es 
notwendig sein, die Tollmiensche Theorie fiir dieses spezielle Problem in der angedeuteten Rich- 


tung zu erweitern. 


(Eingegangen am 9, August 1952.) 
Anschrift des Verfassers: Dr. W. Wuest, Gottingen, Bittingerstrafe 6/8. 
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Strukturmodelle und Systematik der FlieBvorgiange. 
Von E. Mewes. 


1, Einleitung. Unter Last zeigen die verschiedenen Stoffe die verschiedenartigsten FlieB 
vorginge. Viele Arbeiten aus der Rheologie (= FlieBkunde) befassen sich damit, wobei abe 
noch vielfach eine Kluft zwischen der theoretischen und der praktischen Erforschung festzu 
stellen ist. Versuche, systematische Angaben zu machen !, sollten die beiden Zweige einander 
niher bringen. Die Methode, die FlieBgesetzmafigkeiten an Hand von mechanischen Modellen| 
zu erkliren, wird nicht alle Erscheinungen erfassen, doch lassen sich danach verschiedene Eigen 
schaften einordnen. 

Besonders iibersichtlich ist der Aufbau, solange nur parallel gerichtete Federn und Dampferp 
mit linearen GesetzmaBigkeiten als Elemente angesetzt werden. Mit der Belastung t und deri 
VerformungsgréBe y (Abb. 1) gilt dann fiir Feder allein 7 = cy und fiir Dampfer allein r= Dy. 
Ein D vor einer Variablen bedeutet Ableitung nach der Zeit t, so daB Dy die Verformungsge- 
schwindigkeit ist. Tragheitseinfliisse zu beriicksichtigen, war bei den rheologischen Vorgangen| 


i. a. nicht erforderlich. 
C T 30 
4& 
C7 
20 
N DY, 
T 
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rae Abb, 2a, Zusammenhang 
“i yy Abbe Ou bintaches zwischen Beanspruchung, Ver- 
Aprile dell 5 formung und Verformungs- 
Abb. 1. Mit der Zeit sich verformender Block S " KeliinsVaaee geschwindigkeit fiir neben- 


und seine Beanspruchung. stehendes Modell. 

2. Parallelschaltung von Feder und Dampfer. Die Elemente kénnen nun verschieden zu- 
einander geschaltet sein. Bei Parallelschaltung nach Abb. 2 — die parallelen Begrenzungs- 
striche mégen kennzeichnen, da® die dazwischen liegenden Elemente, Feder und Dampfer, 
gleich groBe Verformungen aufweisen sollen — tritt fiir eine zur Zeit t = 0 plétzlich einsetzende ; 
gleichbleibende Belastung 1, ein Verformungsverlauf nach der GesetzmaBigkeit 


y=yilt en) (1) 
mit 
Zz 
y= = (2) 
und | 
ge (3) 
auf. 


Je nach den Giéfen fiir die Steifigkeit c, und die Viskositat uw, kann auch die Einschwing- |} 
zeit T', ganz verschiedene GréBenordnungen aufweisen. Ist cS j,, so ist die Einschwingzeit 
sehr kurz, z. B. 1/;9) sec. Dann ist schon nach 1 sec nichts von Verformungsadnderungen mehr | 
zu spiiren. Ist uu, > c,, dann ist die Kinschwingzeit sehr lang, z. B. 100 Std. Dann ist erst bei 
einer langeren Versuchsdauer das Vorhandensein von Verformungsanderungen festzustellen. 
Sowohl fiir c; —> oo als auch fiir 4, — 00 wiirden Verformungen an diesem Elementenpaar ver- 
schwindend werden. Fiir c, = 0 entartet der Verlauf in 


i phar rect JD (4) 


* Z. B. Systematik der elastischen und plastischen Deformationen nach Aufstellung des British 
Rheologists’ Club, siehe Nature 149 (1942) S. 702, auch Z. Kunststoffe 34 (1944) S. 26. 
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far uv, = Oin 


ig ais (5) 


Wird ein System nach Abb. 2, das zunichst eine Verformung y, gegeniiber dem unbe- 
asteten Zustand aufweisen mége, von t = 0 an sich selbst iiberlassen (also ohne auBere Be- 
astung), dann tritt ein Verformungsverlauf auf nach der Beziehung 


; Pie ae (6) 
Es ist dann, ohne dafs eine Belastung vorhanden ist, solange nicht t> T, ist, fiir einige Zeit 
sine zeitlich sich andernde Ve1formung festzustellen. Das Modell steht also unter der Nach- 
irkung der vorher auf- 

? 

Das zeigt auch die fol- 
ende Betrachtung: Fiir 
elie bigen Verfo1 mungs- 
erlauf y(t) kann der 


CG, 


Pt—cqy+uDy. (7) t 


Jiese Gleichung weist 0 


such darauf hin, daB y, , e 
tir dieses Modell der in 
Abb. 2a skizzierte Zu- 
sammenhang Z 
0 /0 \ 
e «=fly, Dy) =\\WH % 
silt. Demgegeniiber gilt : g 


ai. kine 0 0 

vei beliebigem _—Be- ey ee: 

. . 3. Verlaufe fiir zwei Modelle bei einfach vorgegebenem Belastungs- 
o stungsverlauf an die- bzw. Verformungsverlauf, - : 
sem Modell fiir den Ver- 

jormungsverlauf die aus obiger Differentialgleichung ableitbare Integrallésung 


t 


ie = e—4/T, Je e/T; dt. (8) 
My 
ty 
: as Integral, das sich von einem vorangegangenen Zeitpunkt vollstandiger Entlastung ty bis t 
Irstreckt, zeigt an, daB fiir den Verformungsverlauf die Belastungen der vor dem Zeitpunkt t 
tegenden Zeit maBgebend sind. Wir haben also mit einem Nachwirken der Belastung zu tun. 
Viir verschiedene analytische Gesetzmafigkeiten t(t) kann nach (8) y(t) ermittelt werden. Je 
cach Einsatz der Belastungsgesetzmafbigkeit ist dabei z. B. ts > —co oder ty = 0 anzusetzen. 


_ 3. Vergleich von ParalleJschaltung und Hintereinanderschaltung. In Abb. 3 ist fiir eine genii- 
‘iigende Zeitlang gleichbleibende und dann aufhérende Belastung zunichst fiir das Modell mit 
sarallel geschaltetem Feder und Dampfer der Verformungsveilauf aufgezeichnet, darunter 
asselbe fiir ein Modell mit hintereinander geschalteten Flementen. Fiir dieses zweite gilt 
vahrend der Belastungszeit 


1 1 
Si ee aa 9 
taleg Ding!) (9) 
Nach der Belastung ist eine bleibende Verformung festzustellen von der GréBe 
eee dad (10) 
Ho 
nit t, als Belastungsdauer. Als Differentialgleichung gilt fiir dieses Modell allgemein 
IDs £ 
DI a 
ee Te (11) 


Danach kann fiir belicbigen Belastungsverlauf der Verformungsverlauf ermittelt werden, Fitr 
segebenen Verformungsverlauf erhalt man dagegen den Belastungsverlauf aus der aus obiger 
8 
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Differentialgleichung abgeleiteten Integralbeziehung 


Oe fo 
ee oy fs! dt. (12) 
Inshesondere ist fiir Dy=konst.  infolge 
ens ay 
eM” dt= Pe 4+C 
A) 
ae 
T=fyDy+Ccoe ™ . (13) 
Fur t= 0 ist To=bMoDytCeo,. also 
ay 
t=pyDy+(to—MoDy)e ™ ° (14) 
Ist fiir t >0, y=konst., soist 
r= yet (15) 
mit 
i (16) 
£9 


(Abb. 3). Man spricht von Nachlassen der Spannung oder von Relaxation nach Maxwell. Die: 


Zeitkonstante 7’, hei®t Relaxationszeit. Sie ist 4ahnlich aufgebaut wie die damit nicht zu ver- 
wechselnde Einschwingzeit des Modells mit parallelgeschalteten Feder und Dimpfer. Nach- 


ve 
? 
C; « Cy iI 
iy Co /o S 
Ay 
0 he 0 < a We: 


Abb. 4. Vorgegebener Abb. 5. Zusammengesetztes Abb. 5a. Verformungsverlauf Abb. 5b. System 


Verformungsyverlauf. Strukturmodell. nebenstehender Systeme bei aus zwei parallel- 
i plétzlich einsetzender, gleich- geschalteten 
bleibender Belastung. Maxwell-Gliedern. 


wirken und Nachlassen treten bei diesen einfachen Modellen jeweils bei verschiedenartigen 
Systemen auf (Abb. 3). Fiir das Modell mit hintereinander geschalteten Elementen ist Nach- 
lassen bei Dy = konst. nach (14) nur festzustellen, wenn 


T o> fo Dy (17) 
mit 

To = Co A Vo ; (18) 
im gegebenen Verformungsverlauf nach Abb. 4 ist. 


4. Zusammengesetzte Modelle. Fiir das zusammengesetzte Modell nach Abb. 5 ist 


dt A . 
y=2 =f fo eee) roe dt. (19) 
Daraus lat sich ableiten 
it 1 Li 1 (all Pe} 
D2 rae ey Ss 2 | LL ret) (es ere Ot 
aa IP Dy Co D T | (r, = Te: { Cy 7) Dz if i Ts 7 (20) 
Fir + = konst. wird 
] Poe 
D? ee Be eet lt i| 
ea ape? Dare paige (21) 
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woraus 
amy 0 +t (Dye —7)(1 et) (22) 
Mo Ho 


srhalten wird, Fitr ein vorher in Ruhe befindliches System ist 


1 1 
aaa Dies en ee 
Yo Hi Yo Mac eel 
30 daf 
i 1] 1 
eel oe ane ae (1 rae 3] (23) 


erhalten wird, was in Abb. 5a dargestellt ist. Das Glied mit e~/ weist auf das Auftreten von 
rlastischer Nachwirkung hin. Uber die Relaxation bei diesem Modell wird spater gesprochen. 


Als nachstes Modell wird ein aus zwei parallelgeschalteten Maxwell-Systemen aufgebautes 
ystem besprochen ! (Abb. 5b). Es ist y; die Durchfederung der einen Feder. Die Verfor- 
mung des zugehérigen Dampfers ist dann y —yy. Jetzt ist 


eryr= pa(Dy — Dy). (24) 


ral 


Fir das Gesamtmodell gilt 
t=cpe "1 Dye! dt + ene {Dy eM de. (25) 


Fiir die Ableitung der allgemeinen Differentialgleichung wird nun die Heavisidesche Ope- 
atorenmethode benutzt. Der Operator D wird wie ein algebraischer Einzelausdruck behandelt. 
n den Zwischenrechnungen wird D getrennt von den Veranderlichen geschrieben. Aus (24) 
ird erhalten 


yi(er +o D) =n Dy, (26) 
eer Dy bt 27 
fies aegey ee TaD (27) 
Die Gesamtbelastung vom System (Abb. 5b) ist 
T= r+ enn » (28) 
My qr 29 
r=Dy( 5 | ep) 2) 
(1+ Tr D) (1+ Tu D)t = [par (1+Tu D) +n (1 + Tr D) Dy, ey) 
t+(TM+ Ty) Drt Tr Ty Dt = (ur tun) Dy + (41 Tn + en 11) D’y. (31) 
D iese Gleichung hat dieselbe Form wie (20): 
a, Dy +a, D?y=1t+),Dr+b, D*t (32) 
= fr + peo = o> 
dg = a Tr +o Tr = fo 41 eA (33) 
k= ™)+Tr St et Monte 
b, = 11 Tr =T)T,. 


Die beiden Modelle sind also dynamisch nicht verschiedenartig. Fiir dynamische Gleichwertig- 
keit gilt 


Cy=a +en , Ho = fa + La» ; | 
ey Cy (Cy + Cy) (My + Hq)” et My (My + 1) Su Cyq)” | (34) | 
ef ~ (cy Myr — Py Cqq)” : (My CTT. sas cy Ly) 


1 Siehe Holzmiiller und Jenckel, Z. phys. Chem. 1940, Abt. A, S. 359. Dort ist allerdings cy = ¢y; 


gesetzt. 
ang om 
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Man kann auch verwenden 


a, dy 
Lo= => Cc=7> 
fo bay bs 
a, by __ & 
0= ab si 3a 
2% 1 
ee ay Ay Te (Be ie 
1 @y Gy b, — a? by — a2 by’ 1° a, a b, — a? by — a2 by 


bzw. 


b 
calee ame 
on 4 bé 


Fiir y = konst. wird 


1 eb 


th == Cy et + Cir eel 5 (38) 


Da fiir t->0o offenbar t —> 0 geht, steckt hierin Relaxation. Relaxation ist also auch bei Mo- | 


dell nach Abb. 5 festzustellen, wie Nachwirken bei diesem Modell. 


5. Allgemeine Modelle. Das allgemeinste Modell aus parallelgeschalteten Maxwell-Systemen 
istin Abb. 6 angedeutet. Hierfiir gilt 


0 ei (Dy ei dt (39) 
oder 2 
a, Dy +a, Dy +--++ay DYNy=7+), Dr+b, D?r+---+by DXt (40) 
mit 
Sy it = tn (ee ( a ee 
ag = Pa (Ti Tin) far er (41) 
b= T, 6=Mmlo+thilin+:--, 
boy Pn Pa 5 OR 


Eine andere Form des allgemeinsten Modells aus parallelgerichteten Federn und Dampfern ist 
in Abb. 6a angedeutet. Hierfiir gilt 


cy cy 


ze ates ue 
- =E +E" Dre i [ve dt (42) 


i 
v=1 ie 


frde=D47 (43) 
und 
USW. C; Vi + D Vita (44) 
1 n 


Hieraus ergibt sich mit 
Abb. 6. Abb. 6a. 6 
Abb. 6 u. 6a. Allgemeinstes Modell. N=n-+1 (46) 


die Differentialgleichung (40) mit 
=o, A=Mg(T,+T,+T3+--), 4; =Mo(T1 T,+7,T;+---), -*-, 
1 1 1 
b= Dosey sh Tach My ei tette| eh eee 
1 2 3 


—— ee 
— 
fs 
~ 
_— 


1 ] 
y=1(7 tot are 
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ee aC Ren Re Tees hen) ag Se Tay ie, 
eh 
| &=7,(T 7,477, +--)+%T,T. + ees 
i at ee ear ee oC | 


Diese Differentialgleichung (40) ist giiltig fiir beliebig vorgegebenen Belastungsverlauf oder 
Verformungsverlauf. Die Giiltigkeit fiir beide Modelle nach Abb. 6 und 6a kennzeichnet, daB 
lie Systeme dynamisch gleichwertig sind 1, wenn bei dem System von Abb. 6 ebenso viele 
Federn und Dampfer verwendet werden wie bei dem System von Abb. 6a und die aus den ver- 
schiedenen Formen fiir die Beiwerte a und 6 sich ergebenden Beziehungen zwischen den Gréfen 
:und uw beachtet werden. 


6. Verwickeltere Modelle. Es sind nun auch kompliziert erscheinende Aufbauten von Modellen 
»estehend aus parallel gerichteten Federn und Dampfern, die in verwickelter Weise zueinander 
seschaltet wurden, angegeben worden, um zu zeigen, daB es auch kompliziertere Falle von 
Strukturschaltungen gibt®. Es 1a6t sich nun aber zeigen, dafs auch fiir solche Struktur- 
nodelle sich gleichartige Differentialgleichungen (40) ergeben. Dies wird fiir das in Abb. 7 
largestellte Modell nachgewiesen. Dafiir ist, wenn die Verformung der Feder 1 mit y,, der 
feder 2 mit y,, der Federn 1 und 3 zusammen mit 
7, und die Gesamtverformung des Modelles wie 
onst mit y bezeichnet werden: 


| T = [lg(Dy — Dys), 
| T = ¢3(Y3 — V1) + Me(D ys — Dye), 
T= Yi +272, 
11 = ¢3(73 — 1) +a (Dy, — Dy). 


Jas sind vier Gleichungen mit den vier Unbekannten y bei gegebenem 7. Unter Verwendung 
er Operatorenmethode lassen sich die Unbekannten bestimmen, z. B. y aus 


Abb. 7. Verwickeltes System. 


(48) 


és D 0 0 — fs D T 0 0 — lz D 
0 — Cs gD f,D+e, v —¢s —-lgD yD +e, 
ie , (49) 
0 a eg 0 u ea 2 
0 —(mD+e+¢) wD C3 0 —(4D+eq+¢) (4D C3 
joraus erhalten wird 
a, Dy +a, D?y +4, DP y = byt +6, Dr +b, D?r +b; Dt (50) 
nit 


Ay = Cj Co Cz Mg » 

Gy = [(¢y + Cy) €3 fy + (C1 €y + eo + Cy 3) Ma] Ms » 

a = (C, + C2) [My Me Ms » 

| by = Cy €3, | (51) 
by = (€, + ep) €3 fly + (C4 Cg + €4 3 + 2 C5) fla + (C1 + 3) C2 Ma » 

by = (€y + C9) Mi Me +(e, + cy + 3) ly Mg + (¢ + €3) Me Ms» 

by = fy [2 Ms « 


Jie Glieder lassen sich auch umformen auf b,=1, a= fl3,°°* 

Das Modell liefert also auch nichts Neues. Es ist, wie alle anderen derartigen Modelle, nur 
in Sonderfall des allgemeinsten Falles fiir parallelgerichtete lineare Federn und Dampfer, der 
lurch das System von Abb. 6 oder von Abb. 6a verkérpert wird. 


1 T, Alfrey und P. Doty, J. appl. Phys. 16 (1945) S. 700. e 
2 J. M. Burgers, Mechanical considerations concerning viscosity and plasticity. Aus: First Report on 
fiscosity and Plasticity. Verhandl. der Kon. Akad, v. Wetensch. Amsterdam, Afd. Naturk. (I. Sect.) 


1 XV Nr. 3 (1935). 
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Mit den geschilderten Modellen 1aBt sich ein weiter Bereich von FlieBvorgangen erfassen 
Wir haben es in den einzelnen Fallen mit Sonderausfiihrungen des allgemeinsten Falles zu tun, 
wobei die Stoffkonstanten verschiedenste: Werte und GréSenordnungen annehmen kénnen, 
z. B. auch fiir einzelne Glieder unendlich gro®B werden kénnen, womit dann diese EinfluBteile 
verschwinden. Z. B. kann j4g—>cO gehen, wobei dann keine bleibenden Verformungen mehr 
auftreten wiirden, Es ist nun die Aufgabe, bei praktischen Behandlungen von FlieSvorgangen 
die Auswahl des anzusetzenden Modells so zu treffen, da méglichst mit geringem Berech- 
nungsaufwand die tatsichlichen FlieBeigenschaften recht gut wiedergegeben werden. 


x 


T 
Tt 
Cy 
T 
Va 
Abb. 8. System mit elasti- Abb. 8a, Belastungs-Ver- eee Modell Ha aoe ~ 
schen und plastischen Ver- formungs-Diagramm fiir wirkung (rein p eee e ea 5 q 
formungen ohne ZeiteinfluB, nebenstehendes Modell. Verformungen WO - 10. System mit 
ZeiteinfluB). ZeiteinfluB. 


7. Medelle mit Keilen. Die bisher nur aus Federn und Dampfern aufgebaute Modellreihe 
1aBt sich nun doch durch Hinzunahme von weiteren Gliedern erweitern, z. B. unter Beriick- 
sichtigung von fester Reibung. Hierbei wird zunachst das Modell von Abb. 8 erwahnt. Bei 
einer Druckbelastung tritt eine Kinfederung der oberen Feder von der GréBe y, = t/c, ein. Zu- 


t satzlich verschiebt sich der Konus, und 

zwar auch um einen der Belastung pro- 

of portionalen Betrag, da die Normalkrafte 

am Reibungskonus durch die untere Feder 

e proportional mit dem Eindriickungsweg 

* des Konus steigen sollen. Der Konus soll 
RMON ez 


selbsthemmend sein, so daB bei Entlasten 
Abb. 11. Modell fiir den ideal Abb. lla. Belastungs-Ver- nur die obere Feder zuriickgeht. Das Sy- 
plastischen Kérper. Lopmiaiee Diagram wees stem weist also einen elastischen undeinen — 
nebenstehendes Modell, | 
plastischen Anteil auf. Bei Wiederbelasten | 

wird die Feder wieder einfedern, der Konus sich dagegen erst nach Erreichen der Vorbelastung 
bewegen. Entsprechend gilt das Diagramm in Abb. 8a. Ein KinfluB der Vorbelastung liegt — 

also bei diesem System vor, ein ZeiteinfluB nicht. 

Ein Modell fiir den einfachen Fall rein plastischer Verformung ohne elastischen Verfor- 
mungsanteil ist in Abb. 9 gezeich- 
net. Weiter ist in Abb. 10 ein 
System mit Keil angegeben, bei 
dem die Belastungszeit von Ein- 
fluB auf die VerformungsgréBe 
ist. Damit lassen sich die ver- 
schiedensten Uberlagerungen 
vornehmen. 


8, Elastisch-plastischer Korper. 
Abb. 12. Modell fiir elastisch- Abb. 12a, Belastungs-Verformungs- Als machetes syeiers wind dal 
plastischen Kérper mit Verfestigung. Diagramm fiir nebenstehendes Modell. SO: de Saint- Venantsche wieder- 
gegeben. Dabei tritt nach einem 

elastischen Bereich ein plastischer hinzu, bei dem kein Anwachsen der Spannung mehr 
auftreten soll (Abb. 11a). Dieser Idealfall gibt eine Naherung fiir das Verhalten be- 
summter fester Kérper in beschrankten Bereichen wieder. Das dazugehérige Modell ist 
in Abb. 11 skizziert. Der Reibungsstein darin soll andeuten, da®B vor Erreichen einer 
Grenzreibung ty der Stein sich nicht bewegt, bei Bewegung er dagegen stets die Reibungs- 


)\ 


}XXI. Band 1953. Mewes: Strukturmodelle und Systematik der FlieBvorgange. 111 


pgr6Be ty tibertragt. In Abb. 12 ist noch ein System mit elastischem Anteil und_plasti- 
Feobem Anteil mit sog. Verfestigung wiedergegeben. Das Spannungs-Dehnungs-Schaubild da- 
ifiir ist in Abb. 12a wiedergegeben. Zu bemerken ist dazu, daB bei Zuriickgehen der Belastung 
fim Bereich 1 < Ty vollstandige elastische Riickfederung, im Bereich Te <t <2 ty nur Riick- 
pfederung der ersten Feder und im Bereich t > 2 ty auch Riickfederung der zweiten Feder bis 
jauf Yptmax = Tf/Cg auftritt. In Abb. 13 ist weiterhin ein System ,,mit Spiel‘‘ aufgezeichnet, bei 
dem sich die in Abb. 13a dargestellte GesetzmaBigkeit herausstellt, die mit gewissen Ideali- 
fsierungen praktisch bei Stahlzugversuchen auftritt. Man sieht, wie man durch Hinzunehmen 
verschiedener Besonderheiten immer mehr Varianten erfassen kann. 


Dy 
SS 
eas 
~ 
. id 
Abb. 14. Modell fiir Abb.14a. Zusammenhang 
* ; Bingham-Korper. zwischen Belastung und 
Abb. 13. Modell mit FlieBgrenze Abb. 13a. Belastungs-Ver- ; Verformungsgeschwin- 
und Verfestigung mit Spiel. formungs-Diagramm fiir digkeit fiir nebenstehen~ 
nebenstehendes Modell. des Modell. 


| 9. Systeme mit FlieBgrenzen. Durch Hinzunahme eines Gliedes mit Couwlombscher Reibung 
jkann das Auftreten einer FlieBfestigkeitsgrenze beriicksichtigt werden. Das Gesetz von 
| Bingham 
Tit Veit Dy =0, | 
(r= 1% +H Dy) (52) 
t>uy: Dy=te—y[ >” 


\(Abb. 14a) wird durch das Modell von Abb. 14 erfiillt, bei dem dem Dampfer ein Reibungsglied 


}parallel zugeschaltet ist. 
t 
qt 
_ t 
0 a al 


ft 
a 


ig 
c + i 
@ M, 7) 
Abb. 15. Kelvin- Abb, 15a. Verformungsverlauf an Abb. 16. Maxwell-Bingham- Abb. 16a. Spannungsverlauf an 
Bingham- Modell. nebenstehendem Modell fiir zeit- Modell. nebenstehendem Modell fiir zcit- 
weilig gleichbleibende Belastung. weilig gleichbleibend. Verformung. 


Kinige weitere Kombinationen sind in Abb. 15 und 16 wiedergegeben. Bei dem System von 
Abb. 15 mit parallelgeschalteten Elementen (verallgemeinertes Kelvin-Modell) bleibt nach 
einer Belastung stets eine bleibende Verformung erhalten, deren Gréfe bei hinreichender Be- 
JastungsgréBe und -dauer t,/c wird (Abb. 15a). Bei dem System von Abb. 16 (verallgemeinertes 
Maxwell-Modell) geht der Relaxationsvorgang nicht bis auf t = 0 fiir t + co, sondern auf Ty 
zuriick (Abb. loa). Fir |t| < ty tritt bei y =konst. keine Spannungsdinderung auf. Ent- 
sprechende Verlaufe kénnen auch bei Uberlagerungen festgestellt werden. 

Die Unterscheidung der FlieBvorgange ist nicht nur eine Frage der verschiedenen Systeme, 
sondern auch der GréBenordnung der einzelnen Stoffbeiwerte (c, ). Es kann i. a. der all- 
gemeinste Fall als vorliegend angenommen werden, wovon lediglich fir einige extreme Bei- 
wertgréBen (~ 0 oder co) einzelne Glieder praktisch verschwinden, wodurch dann etwa die 
einzelnen Sonderfalle (Maxwell- oder Kelvin-Kiérper) zum Vorschein kommen kénnen. 

- Eine mathematische Untersuchung von nichtlinearen Ans&tzen fiir die Kigenschaften der 


Elemente steht noch aus. 


(Eingegangen am 10. August 1952.) 
Dr.-Ing. E. Mewes, Braunschweig-Vélkenrode, Institut fiir Landtechnische Grundlagenforschung. 


1 M. Reiner, Deformation and Flow. Lewis, London 1949, 
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Zur Theorie und Anwendung der schiefen Platte*. 


Von J. Krettner. 


1, Einleitung. In der Bautechnik ist neuerdings in mehreren Arbeiten! das beim Briickenbau 
auftauchende Problem der Durchbiegung einer schiefen Platte behandelt worden, bei der 
jedenfalls zwei Achsen eine schiefe Lage aufweisen, wahrend die dritte Achse senkrecht zur 
Ebene jener beiden Achsen angenommen wird. Es erscheint daher zweckmaBig, die Formeln 
der Elastizitatstheorie ganz allgemein auf schiefwinklige Koordinaten zu reduzieren, um allen 
auftretenden Sonderfragen gerecht zu werden. Die direkte Methode der Koordinatentrans- 
formation erweist sich als ziemlich umistandlich und oft auch schwierig. Ein in jeder Weise 
befriedigendes Hilfsmittel, bei dem die iibliche Gleichgewichtsbetrachtung am Elementar- 
kérper vermieden wird, mu die Methode der Tensorrechnung gelten. Im Folgenden sollen 
nun die Grundformeln fiir schiefwinklige Achsen zusammengestellt und auf die schiefe Platte 
angewandt werden, in der Absicht, die Sonderbehandlung der Durchbiegung einer schiefen 
Platte vorzubereiten. 


2. Vektor im schiefwinkligen Koordinaten, Sind x1, x? x3 die Koordinaten des Raum- 
punktes P in einem kartesischen Koordinatensystem mit den Einheitsvektoren ¢;(i=1, 2, 3), 
wobei ¢; + ¢, = 1 bzw. 0 fiir i=k bzw. i-k ist, so ist der Ortsvektor r des Punktes Pin bezug 
auf dieses System in der Form ry = 2’ x‘e; darstellbar. In diesem kartesischen System werden 
nun drei beliebige Vektoren a, (k=1, 2, 3) mittels der Gleichungen a, =) aj e; (k=1, 2, 3) 
als Grundvektoren eines schiefwinkligen Koordinatensystems (u!, u2, u3) mit den Einheits- 
vektoren 1, = ,/|0,| in Richtung der schiefwinkligen Achsen gewahlt. Dabei ist voraus- 
gesetzt, dafs die drei Vektoren q, von cinander linear unabhangig sind, ihre Determinante 
also verschieden von Null ist: A =(a,, Q),0;) +0. Der zum Raumpunkt P fiihrende Orts- 
vektor t ist in diesem schiefwinkligen System durch 


uk 


3 3 
CS hy Oy, mit uk — —— 
k=l Lat |ay,| 


(1) 
gegeben. Fiihrt man noch das zum schiefwinkligen System gehdrende reziproke (u,, U,, Us)- 
System mit den Grundvektoren q* (k =1, 2, 3) ein, so ist definitionsgemaB: 


1 Xs qa? = XM 


(Ay My M5)” (Ay My Ag)” 


a, Xa, 
(4 My A5)” 
wobei a;:a* = 1 bzw. 0 fir i=k baw. i+k ist.-Die drei Einheitsvektoren im reziproken 


System sind dabei durch n* =q"/|a*| gegeben. Im reziproken System hat der Ortsvektor 1 
die Darstellung 


— 


3 3 
r=) yn ne = Sua mit Seer (2) 
b= 


Bezeichnet man das skalare Produkt a; - a, mit Siz bzw. a’+ a* mit gi*, so erhalt man fiir die 
Projektion des Vektors t auf die schiefwinklige Achse a, die Beziehung 1 - a, = SY) u'g;, = up, 
fir die Projektion auf die reziproke Achse dagegen t+ a* = > u; gi* = uk, Damit ist zugleich 
der Zusammenhang zwischen den sog. kovarianten Komponenten uw, und den kontravarianten 
Komponenten u* des Vektors t gegeben. 

Fiir das Differential des Ortsvektors t = Sulla; bzw. r= Su; a' erhalt man dr = >i du‘ a; bzw. 
dy = S'du; a‘, woraus durch skalare Multiplikation mit den entsprechenden Achsenvektoren 
fiir die Differentiale der Komponenten dui —ai- dy bzw. du; =a; + dr folgt. 


* Die Anregung zu dieser Arbeit verdanke ich den Herren W. Miiller, Miinchen und K. Kammiiller 
Karlsruhe (Institut f. Beton und Stahlbeton),. 
1 P. Lardy, Schweiz. Bauztg. 67 (1949) S. 207 und 419; H. Favre, Schweiz. Bauztg. (1942), S$. 35, 51, 60; 
H. Favre, Bulletin Techn. de la Suisse Romane (1946), S. 321; I. Nielsen, Beitrag zur Berechnung schief- 
winkliger Platten, Akademie der Techn. Wissenschaften; Vogt, Beitrag zur Berechnung schiefwinkliger 


Platten, Diss, Hannover 1939; Rongveld, Uber das Randwertproblem schiefwinkliger Platten, Diss. 
Berlin 1945, 
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_ 3, Gradiententensor im schiefwinkligen System. Ist nun ein beliebiger Vektor p =» (ul, u*, u’) 


| 3 
durch p= 5’ v' a; gegeben, so ist sein Differential, wenn wir noch fiir die partielle Ableitung 


a1 
die Abkinzung 0p/du' =v; einfithren 
3 


Goce ete = 


i=1 i 


3 3 3 
Dv; (dt: a) = Sdr-ot'y;= dr: Yoiy,=are- B, (3) 
= 


i=1 i=l 


wobei mit % der Tensor zweiter Stufe, der sog. Gradiententensor des Vektors ) bezeichnet 
wird. Setzen wir die Komponenten des Vektors b ein, so ergibt sich: 


3 3 : 3 dub 
B= dai“ eG >" oF a, 
4 Oui \. ; Out 
Pos at eal 

Wenn wir den Vektor » in der Form ) = (Uy, Uy, Us) = Sv; a" nehmen, so wird der Gradienten- 
ensor von Ci! 


3 3 


3 3 
B= Dat Sud(> nat) =D’ aah as () 


_ 4, Deformationskomponenten im schiefwinkligen System. In einem elastischen Medium 
sollen nun den einzelnen Punkten durch irgendeine Deformation die Verriickungen ) zu- 
yeordnet werden. Exfahrt der Punkt P, die Verriickung 9, so erfahrt der von P, um dt ent- 
ernte Punkt P die Veiriickung p=, -+dr- &. 
_ Bezeichnet man mit ¥ den zum Tensor % konjugierten Tensor 

ron dv 


ee Oui a, (9) 


fo kann man den Tensor & als Summe der beiden Tensoren D = ~ (B SF ¥) und Qf = $(8— ¥) 


iarstellen, der symmetrische Teil ® bewirkt die reine Deformation des Bereichs, wahrend 
fer antisymmetrische Teil 2 eine infinitesimale starre Drehung des Bereiches um P bestimmt. 
Js ist also stets 8 —D+Y%. Wenn man also vom Gradiententensor ¥ den antisymmetrischen 


Teil 2{ abspaltet, erhalt man den Tensor der reinen Deformation 9, und es ist 


3 3 
1 a 1 , : 1 BI se 
a ay (5) Ze felt Es At \ eae AY i 
ee as et) = OF ai (Ai Oe + OKO’) « (6) 
jkalare Multiplikation von D mit dem Vektor a, ergibt den Vektor a,.°®, der durch noch- 
halige skalare Multiplikation mit a, die skalare GroBe 
3 


1 dvk dvk 
ee Davie O40,“ O) = (Fi 0 Gu + Ben) (7) 
kee 
irgibt, die fiir 0 =g in 
3 
dvk 
Doo = > Sabet (8) 
k—1 


ubergeht. Damit ist dann die Dehnung in Richtung des Vektors @, durch 


1] , Do 
| Egg = a (9) 
ind die Winkelausweitung der beiden Vektoren Q, und a, durch 


aS 1 
\f00 800 — 30 


Doo Do 
ZaD: oe 4 
| Ps beo( Fy : 5a) 


5. Der Spannungstensor in schiefwinkligen Koordinaten. Betrachtet man ein Parallelflach 
pit den Kanten parallel zu den Achsen des schiefwinkligen Systems, so sind die flachen- 
ormalen Einheitsvektoren die reziproken Einheitsvektoren n* = a /|a"| des schiefwinkligen 
! ystems. Es sei ); die reduzierte Spannung, die auf das Flachenelement wirkt, welches auf 


(10) 


mC 


. 
| 
legeben. 


Q 
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dem Vektor 1‘ senkrecht steht. Sind p,;, die Projektionen dieser reduzierten SP HaTate auf 
die Achsenrichtung n*, so gilt n*+ p;= pix, und der Vektor Pe, kann in ae Form Pi= 2 Pir Ny, 


angesetzt werden. Der SPannungeteneos & selbst ist © = s p n= Sy Yin 1; Me. 


i=1 ie 
Dasit— \a*| Q, ist und sre e;, also @)* Q,= aj, so erhalt man fiir die Spannungs- 
i=l 
komponenten in bezug auf die Achsenrichtungen folgende Werte: 


a= te (Nea SG) = 1; Sipe: At;). 10, == 10 2 Per We = > Pok (1% * Mk) = Poo» 


A tet | 
avt= ae- (av G) =a" al Mi) = & Pik Jae] Ja] (a2 + a,) (a7 + a) = Pos [0°] [071 » 
Ap g= Ap * (Ag° S) = @° Dy Pik (Me my) Ma =D pan (Qe tts) les th) = & Pix || [O*] Boi Be 
Spice aA Ft i) 1 = & Pir (Co * Mi) (ee th) = & Pix [ol [a"| a ay, - 


Fiir die Projektion der reduzierten aoe p; auf die Achsenrichtung a? des reziproken 
Systems erhalt man den Wert 


3 3 
ae pi = S pix (02+ te) = Y pin |a*| (02+ Oe) = Pio [0°] = afo- (11) 
k=1 k=1 
Bezieht man die einzelnen Spannungskomponenten auf diese GréBen, so wird 
a * 3 s 
avo Se ? Cpaae |a2| Apo » ao= > || Goi Bok az, » = 5 [ai al aga aj,. (12) 
Vgee i,k=1 ijk — 2 


Der Spannungstensor © lautet dann fiir die verschiedenen Systeme 
S =) a** 1; 1, =a .o 8 dy pease nena as = ey) Gi x CuiGrr. 
i, k i, k i, k i,k 
6. Die Gleichgewichtsbedingungen im schiefwinkligen System. Die Bedingung des Gleich- 
gewichts deformierbarer Kérper zwischen den dufBeren Kraften und den inneren Spannungen 
in irgendeinem Punkt des Spannungsfeldes lautet bekanntlich div © + # = 0. 
Fir den Gradienten des Spannungstensors © erhaélt man, wenn dieser in der Form 


3 
G= J aa; 
i,k=1 
angesetzt wird, 


3 
pad = 45 = SaegS— DS) Hreae. 


=1 
Durch Verjiingung entsteht daraus die ‘Oonret des ert ear 
3 3 3 
; daib daik 0,ik 
d = ° = ‘a: = — ee 
w6ay-G= SY Mee aa= Tem = > (So) 
a Oy (= us = aN 


Wird der Kraftvektor § in der Form § = 5) K* a, angesetzt, so lauten die Gleichgewichts- | 
bedingungen in Richtung der scab tits 4 Achsen: | 


* Qaik ke 
ga tKh=0 (k= 12,3). (13) 
Dale P= Gna poy a!) @ a*,, ist, kénnen wir dafiir auch schreiben 
3 9 * 
. ay is 
> Jai 2 4 KE 0 (b= 1,2, 3). (14) 


i=l 
7. Das Hookesche Gesetz in schiefwinkligen Koordinaten. Der Zusammenhang zwischen 


dem Spannungstensor © und dem Deformationstensor D ist durch das Hookesche Gesetz in 
seiner allgemeinen Fassung gegeben: 


G6 =2uD +ADS. (15) 
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3 3 
Darin bedeutet § den Einheitstensor J = Y gi* a; a, baw. J = » gi,aia* und die GréBe D 
ies oY 
den Skalar des Deformationstensors D. Wird die Konstante von Poisson mit x bezeichnet, 
so sind die beiden Konstanten von Lamé durch 


Bee Be k 
Dey veces und A= 25s 


gegeben. Durch zweimalige skalare Multiplikation mit den Grundvektoren a, bzw. a* erhalt 
man aus dem Hookeschen Gesetz die allgemeine Spannungs-Deformations-Bezichung. Fiir 
S Seas %, D = Dd e'* a; 0, mit D= >» e* 2, und § = > g'* 0; wird 

i, i,k 


i,k i, k 
| GF (0"> ©) = 2u 0" (a-~D) +AD ak: (a?- 8) 
oder 
| 3 
aik — 2 ygih ta{ 3s e* gus) 8. (16) 
Fir joes 


C= > 4,,.0°0", D => sa 0° und Se a a8 
i,k ok om: 
erhalt man in entsprechender Weise 
A, ° (A;° ©) = 2 ay: (a;°D) + AD a, - (a; - 9) 


oder 

Aik =2uén+aD gir mit De PSie ne: (13) 

i, k 
Durch Auflésen der allgemeinen Hookeschen Bezichung (15) nach dem Deformationstensor D 
ergibt sich 
le: ou 
D= oa [6—AaDS]. 

Man kann nun noch den Skalar D durch den Skalar S des Spannungstensors ausdriicken, 


wenn man (15) auf beiden Seiten verjiingt. Man erhalt S = 2D + 3AD=(2u+2A)D und 
damit die beiden Gleichungen, wenn man noch den Wert A|(2u +34) durch x ausdriickt, 


it LEE gh = 5 git Init. ge = a pie 


i,k 
bzw. (18) 
Gik at ai — ES oi mit See ps Ot Os 
ik 


8. Der Differentialausdruck von Laplace in schiefwinkligen Koordinaten. Setzt man fiir 
é é Ow O?w 5 

‘ ‘ a taf off Ars) ey 
a partielle Ableitung der Funktion w=w(u', uv’, u*) aatat e und ui dub Wik» 80 ist 


der Gradiententensor der Ortsfunktion 


3 
grad w= jy =>) 0 w;- (19) 
i=1 
Fiir den Operator Aw =div (grad w) = div tp folgt damit in schiefwinkligen Koordinaten 
3 3 3 3 
| divi = Sater, = Sab-( Sak win) = 3 gh win. (20) 
) ea k=1 i=1 i,k=l 


Nochmalige Anwendung des 4-Operaters ergibt schlieBlich den Differentialausdruck von 
| Laplace: 

| 3 

| AAw = div (grad Aw) = divb = 3) 0°: ba 


| o=1 


= » ar (Har (A1)e) = > ge7(Aw)o, (21) 


@ eo=l 


3 


3 
= Ds se( 3 Wines) — By LOSE Gey age 


| i,k, Q,0=1 
| 04 w 


wobei Wego = 
. 


- ist. 
Oui du’ duedus 
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9. Die schiefe Platte. Mit Hilfe dieser Formeln sollen nun die fiir die Berechnung der 
schiefen Platte notwendigen Gleichungen ermittelt werden. Die eine Kante der schiefen | 
Platte sei dabei senkrecht zur Ebene der beiden anderen Kanten angenommen. Die Achsen | 
des Koordinatensystems sollen in den Kanten des Parallelflachs fallen und auf ein kartesisches | 
Koordinatensystem mit demselben Nullpunkt bezogen werden. Der Ortsvektor Y zu einem 
Plattenpunkt ist dann in den beiden Formen darzustellen: 


3 3 
Li) SAG; und cee ey (22) 
iat 
Es ist (Abb. 1) 


1,=e,, OG, =singe, +cos@me,, My = C3, (A, AgM3) =cosy, 


J 
vee ; S722 
a =e; — tg Per, Oe rege es ae = e3, 


811 = Ugg = 833 = 1, Biz = 821 = SIN PY» Biz = 831 = S23 = Bs2 0, 


fil sense ee Ce ee 13 23 31 $3 =.) 
3 8 cos? yp” 8 eee. g cos?” & & 8 8 > 
a= law 0550) -- 0). of = sing, “6h cos ts a a—0, @=0, a= 
Fiir den Zusammenhang der schiefwinkligen und rechtwinkligen Koordinaten folgt aus (22) 
3 3 
RPE ey, == a (Gey) bzw. ub = ya = » a (e > a*) 
i=in = 
also 
ai—ultsingu?, 2 =cosgu’, ere 
23 
ui = a21_tgo x, ue = : x, te (23) 
cos p 


Setzen wir den Vektor der Verriickung in der Form 
3 3 3 
b= St: = DvaG= yw 


yh pu Hiss 


an, so erhalten wie fiir seine kovarianten bzw. kontravarianten Kom- 


Abb. 1. Koordinaten- ponenten 
system fiir die schiefe 


3 3 
Platte. = See 
is = 0° 0; =>) 2, (GG bzw. of ==) 0 ==) > (0 ee 
i=1 =a 
also aa v,=singvitcospv?, v,=0v%, | 
ide 24 
vi=vl—tgpe, v2 =—~—®, 3 = BS | ( )q 
cos Y 


Aus der Formel (8) folgt, da in dem Falle der schiefen Platte alle g;; = 1 sind, fiir die Deh- — 
nungen in Richtung der Achsen des schiefwinkligen Systems : 


Gul ay oe Oe Me us ave d 
Eu = aa 1 SID Pa? Ex. = SIN P 3 Twas £33 — 40 (25) 4 
Da g13 = 81 =0 und go. = B39 =0 ist, ergeben die Forme!n (7) die Winkelausweitungen 
Ov Ov : Ov? 
= 63 = — & =2 Dy = s+ 5, +singsG, | 
dv? GUE, sane. dv! (28) 
E93 = — 3. = 2 Das = Auk Oni ae 
Fir die Winkelausweitung — ¢,, folgt mit 
dv) dv? : dvi dv") 
2 Da = Fat gai +0 0 (Fa + 5a] 
aus (10) 
dul de® 
== OD SSS ET = 008 9 (Sia + 5a) (27) 


Um noch den Zusammenhang zwischen den Dehnungen und Winkelausweitungen im recht- 
winkligen und schiefwinkligen System zu erhalten, miissen mittels den Formeln (24) die 
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Komponenten v' durch die Komponenten 1! ausgedriickt und die Ableitungen nach den ui 
nach der allgemeinen Forme] 

| dvi + ave Oxo 

due — ue ie duo 


j 


berechnet werden. Werden die Deformationskomponenten in bezug auf die kartesischen 
Achsen durch 

= dvi = dvi duh 

Eii een und Eik =o 


ausgedriickt, so ergeben die Formeln (25) die folgenden Beziehungen: 

Er = E> Eng = Sin* M ey) + sin Ycosp eq + COs? Mey, Egg = E55, 

fia en — C08 P[ (8,4 — €n2) 81D Y + €1,¢08 P]; fig = Eg) = Cig = En (28) 
£93 = &32 = E13 SIN P ++ &y3 COSY. 


[Fir die Projektionen der Spannungen auf die Achsen des schiefwinkligen Systems er geben 
die Formeln (12) die Werte 


1 ; 
Mee cos [af +2 sin p af;+ sin! y a3z], ay. = Ep sin? p ay, + 2 sin ~ af, + ado], 
i a3; , 
1 ; (29) 
Fie = 4%, = —— [(af,+ @},) sing + (1 + sin? q) af], 


~ cosy 
—— d= gtd [lags eis 2716" G5, — 0,57 ail [sin py a}3-+ as]. 
‘EM 31 cos 5 he cos ~ 


‘Die Projektionen der Spannungen auf die Achsen des rechtwinkligen Systems erhalt man 
ebenfalls aus den Formeln (12): 


_— 1 . . * = * _— — 
a, = —— [aj + 2 sin ¢ af, + sin? p ay], Azo = COS P 492; G33) — 033 9 
cos @ (30) 
~ 1 : as = 
= ‘ * — * * ett, a gt See 
Ua eae ay) sin P43, yg = cos p (aj3-++ sin P 433), 4g, = Ago = 093 = 432. 


urch Auflésen dieser Gleichungen nach den a}, folgt der Zusammenhang zwischen den 
jProjektionen auf die rechtwinkligen und reziproken Achsen: 


+ a ae = ogi Ga at, = Gy. — tg pa aj. = —— 4, 
Gj] = COs P a, + —— Ayo P U2; 12 12 229 22 (31) 


cos p cos Pp 

at; = cos @ G3 — SIN P Gog, 933 = Ang, 033 = Gg. 
ie Formeln (14) ergeben die Gleichgewichtsbedingungen in Richtung der Achsen des schief- 
jwinkligen Systems: 


1 2 da*,, 


+ Kk =0 (ke "1, 25 8). 


COS @ 4 dui 
cl 


E k ‘ ' 
Site l+k [es Gq —2k) cos? p (€11 + €22 + &33 —— 2 sin ~ a) | ’ 
eieetiict Ie le — ba Pp (E11 + &22 F &33 — 2 sin ex) | ’ : he 
(eas le Ves 
E 


UTE) a ey ae ? 


E 
Ao3 = 7 7, 23 ¢ 
1+k 


‘fir den Gradienten, Aw-Operator und den Differentialausdruck von Laplace ergeben die 
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Formeln (19), (20) und (21) die Werte 


3 
to = gradw= >a", 
a | 


1 Cha) Oo? w : 0? w 0?w 
Aw = | aaa + a 2 sin ~ aa a | Ou? 9 


cos? p 


(33) 


1 Ge of ; o*w : 04w o*w 
=| au | aa aaetleg 2 sin®@) sagan — 4809 (aves | du? dul? }|’ 


‘a a4 w Ow a] 


| du? du? A du2? gus® | corey, du4 


cos? 


Dabei fallen fiir diinne Platten die Ableitungen nach der dritten Koordinate noch weg. 
Die Formeln (23) liefern fiir die Ableitungen der Koordinaten die Werte 


du Out Ou? Ou? 1 


Get gat, PE Get eee oes 
und damit erhalten wir fiir die Ableitung der Funktion w = w(x}, x?) =w(u!, u?) selbst 
Ow ow Ow Ow 1. dw Ow , Pw 
dat du” Ox? 8 Pau | cos p du?’ axt2 ' Qui?’ 
0?w . 0*w 1 la Dele 


axe? te? du ete cos p dui du? | cosy du?” 
Setzt man diese Werte in die bekannten Formeln fiir rechtwinklige Systeme 
Ex | ew Ow - Ex' | dw Ow x Ex dw 
1— 2? leat Me 22 }— leat Lb ayae >» eS 1-By oxox 
ein, so erhilt man aus den Gleichungen (31) unter Verwendung der Gleichungen (33) schlieB- 
lich fiir die Spannungskomponenten in Richtung der reziproken Achsen die Werte 


One 


Eu aw Shia Ew | yy Bw 
alee (1—p*) cos [av ey) a ‘ O22 (1p), con. Aiea ieaes dul ’ BY. 
* Eu3 a A 1 0? w ( ) 
W112 = (1—p?) cos Al Sees os (u ) du! du? | 
Durch Integration tiber die Plattendicke kénnen nun daraus die Momente berechnet werden: 
+ Yh 
Eh’ Ow : 0? w : 0?w Pw 
ae ee 12 (1—?) cos® = lau yee au? etic? du! du* i du? ooo | : 
—t%h 
- as ‘ Eh’ : Ow , dw ; 0? w Cu 
M3.= | a3, uu? due = dap) coargh™ ® sai8 ais el 2 sin y aul due Se ayia COS | A (3a) 
Wp 
aie Eh? a aw 
ee & 8 que Sy 
Mf,= | ay, du = 12 (1—p,?) mel Pe ae es ee PET 


—h/2 


: Ow 0?w 
+ (1 + sin? ¢) Salou Fier cos? | e 


Eh? 


Setzen wir 12 (1 — ») cos i 


° 


so erhalten wir fiir Summenausdruck 
M* = M#, + Mf, +2 sing Mi, = K Aw (36) 
und damit 
AM* = KAAw. (37) 
Damit sind die wesentlichen Formeln, die zur Berechnung der schiefen Platte notwendig 
sind, hergeleitet. In einer spiteren Arbeit soll die Durchbiegung der schiefen Platte fitr 
besondere Belastungsfalle naher untersucht werden. 


(Eingegangen am 25. August 1952.) 
Anschrift des Verfassers: Dr. Josef Krettner, Miinchen, Gentzstr. 5/III/1. 
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Anwendung der Matrizenrechnung auf betriebswirtschaftliche Aufgaben 
| Von O, Pichler. 


1, Einleitung. Fiir viele gréBere Industriewerke ist es kennzeichnend, daB sie aus vielen 
einzelnen Teilbetrieben bestehen, die unter sich stark gekoppelt sind. Der Produktionsablauf 
eines solchen Betriebssystems ist nun keineswegs immer stationar; denn die zur Verarbeitung 
gelangenden Rohstoffe andern sich, das Produktionsprogramm muf den jeweiligen wirtschaft- 
lichen Anforderungen angepaft werden, und es kénnen Betriebsstérungen in den verschiedenen 
Teilen auftreten, die sich auf das Gesamtwerk auswirken. In allen diesen Fallen muB die 
Betriebsleitung in der Lage sein, die Folgen dieser Einfliisse auf den Produktionsablauf 
zahlenmafiig zu iibersehen, um unverziiglich eingreifen und eine unter den vorgegebenen 
Bedingungen jeweils wirtschaftlich optimale Lage herbeifithren zu kénnen. Aber auch bei 
der Projektierung neuer Werke und bei wirtschaftlichen Planungen auf weitere Sicht ist die 
zahlenmaBige Kenntnis des Betriebsablaufes eine notwendige Voraussetzung. Da sich nun 
die Ablaufe der einzelnen Teilbetriebe im allgemeinen einfach iibersehen lassen, macht es 
‘keine grundsatzlichen Schwierigkeiten, die Auswirkung irgendeines Einflusses auf die ver- 
ischiedenen Teilbetriebe rechnerisch weiter zu verfolgen; aber solche Rechnungen sind meist 
sehr mithsam und liefern das Ergebnis nicht in der gewiinschten allgemeinen Form. Und hier 
sbietet sich nun, wie in folgendem gezeigt werden soll, die Matrizenrechnung Koks — H10 

als ein fiir soleche Probleme geradezu pradestinierter Kalki] an. Bevor 
die Methode und ihre Anwendungsméglichkeiten aber in méglichst all- 
igemeiner Form dargelegt werden, erscheint es geboten, sie zunachst an Wassergaserzeugung 
icinem konkreten, technischen Beispiel zu erlautern. Als solches Beispiel 
jwurde die groftechnische Herstellung von Wasserstoff gewahlt, die u. a. 
ibei der katalytischen Druckhydrierung von Kohlen, Teeren und Mineral- Entschwerelung 
jlen eine entscheidende Rolle spielt. Was die im folgenden benutzten 
ormeln der Matrizenrechnung anlangt, so sei auf das Buch von Zur- 
imiihl+ verwiesen, dessen Bezeichnungsweise weitgehend iibernommen Konvertierung 
jwurde. 


} 2, Darstellung von Betriebsvorgangen durch Matrizen bei Einzelbetrieben. 
In vielen Fallen wird bei der Wasserstoffherstellung vom sogenannten 
|Wassergas ausgegangen, das in ,,Gasgeneratoren“ etwa aus Steinkohlen- 
ikoks oder Braunkohlenschwelkoks mit Wasserdampf erzeugt wird und 
jaus Kohlensdéure, Kohlenoxyd, Wasserstoff, Stickstoff, Methan und Druckwasserwasche 
\Schwefelwasserstoff besteht. Dieses Wassergas muf nun in die fiir den 
}Verwendungszweck erforderliche Zusammensetzung iibergefiihrt und auf 
iden dazu notwendigen Druck gebracht werden”. Zu diesem Zweck wird Kompression auf 300 at 
xs zundchst entschwefelt (Abb. 1), sodann wird das Kohlenoxyd des 
s;Wassergases zum gréBten Teil in der ,,Konvertierung’: zu Wasserstoff 
fund Kohlenséure umgesetzt. Das so erhaltene Gas wird dann auf etwa Kupterwiische 
30 at komprimiert und hierauf zur Entfernung der Kohlensdure mit 
:Wasser gewaschen. Es folgt eine weitere Kompression auf etwa 300 at, 


und dann wird in einer Endstufe der Rest an Kohlenoxyd entfernt und — (0-Aéekgas eee, 
Pischg 


Kompression auf X0ot 


| ‘ 5 ‘ : é : 
damit ein fiir die Hydrierung geeigneter hochprozentiger Wasserstoff, 

‘ : Abb. 1. FlieBbild der 
jdas sogenannte ,,Frischgas‘‘, erhalten. Das entfernte Kohlenoxyd mit — wasserstofferzeugung. 
jseinen Begleitkomponenten wird gelegentlich als ,,CO-Riickgas“ be- 

Peichnet, weil es, wenn es nicht anderweitig Verwendung findet, in die Konvertierung zu- 
dctickgefiihrt wird. 

Es soll nun zunachst an dem Beispiel der Konvertierung die matrizenmafige Darstellung 


fines Betriebsvorganges gezeigt werden. In der Konvertierung wird Kohlenoxyd unter Zugabe 
| 


1 R. Zurmiihl, Matrizen. Berlin/Géttingen/Heidelberg 1950. 
2 Naheres s. etwa Krénig, Die katalytische Druckhydrierung von Kohlen, Teeren und Mineraldlen, 
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von Wasserdampf zu Kohlendioxyd und Wasserstoff umgesetzt gemaB der Reaktionsgleichung: 
CO --H,0 = CO, +- H,- 

Aus 1 Nm? CO im Eingangsgas, dem ,,Durchsatz‘ (in diesem Fall entschwefeltes Wassergas), 
entstehen also unter Zugabe von Wasserdampf 1 Nm? CO, und 1 Nm? H, im Ausgangsgas, 
der ,,Ausbeute‘*. Vollstandigen Umsatz vorausgesetzt, verschwindet also CO, wahrend die 
Mengen an CO, und H, sich je um die umgesetzte CO-Menge erhéhen. Stickstoff und Methan 
bleiben bei dem Proze& praktisch unverdndert. Bezeichnet man nun die einzelnen Kom- 
ponenten (etwa in Nm’) im Eingangs- bzw. Ausgangsgas gemab 

a ie a ee ns 


Eingangsgas Ausgangsgas 
(Durchsatz) (Ausbeute) 
Nm$ Nm? 
SOR VN ed ee ee 
cO, d, a, 
CO d, ay 
H, d, a, 
N, d, a, 
CH, d; Geshe 
so lat sich der UmwandlungsprozeB, wie folgt darstellen: 
a,=4,+4,, 
a, = 0 
es d, + ds, 
Cae 4? 
as = d;, . 


Es liegt nun nahe, den Durchsatz und die Ausbeute durch Vektoren gemafs 
D.= {d, |, und a= (a, 


d, ae 

d, a, 

d, a, 

ds | as 

darzustellen und ferner eine ,, Umsatzmatrix“ 11 gemaB 

Wes (us) ly doe, BOO 
Or. O26050% 0 
Ovceky; ple DO 
040500 ater 
0.70, 5 0420.981 


einzufiithren; dann bekommt man fiir die Darstellung des Umwandlungsprozesses die einfache 
Form 


a=lNd 


(Ausbeutevektor = Umsatzmatrix mal Durchsatzvektor). 


Man erkennt ohne weiteres, daB® bei der Konvertierung die Umsatzmatrix singular ist, da der 
zweite Zeilenvektor ein Nullvektor ist. Bei der Darstellung der tatsachlichen Betriebsver- 
haltnisse ist die Matrix noch etwas zu modifizieren, da das Kohlenoxyd nicht vollstandig 
umgesetzt, sondern ein festes, wirtschaftlich optimales Verhaltnis von CO:H, eingestellt wird 
und auBerdem noch gewisse Verluste zu beriicksichtigen sind. Diese Verhaltnisse kénnen 
durch den Ansatz 


a, = u,(d, + ds), 

a, = us(d, + ds) ; 
dargestellt werden, wo u,/(u, ++ us) den CO-Anteil in der CO -++ H,-Menge bedeutet, der durch 
betriebstechnische Ma®nahmen praktisch konstant gehalten wird und (u, +-us) die ,,spezifische 
Ausbeute“ an CO -+-H, darstellt, die nur wenig kleiner als 1 ist. Die umgesetzte CO-Menge 
betragt somit 


d, — dy = (1 — u,) dy — und, , (1) 
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um welche sich die CO,-Menge erhéht. Die Verluste an den tibrigen Komponenten kann man 
ebenfalls durch ,,spezifische Ausbeuteziffern“ u,, u4, us beriicksichtigen. Die Umsatzmatrix, 
die jetzt zum Unterschied von noch zu erérternden anderen Umsatzmatrizen mit U, bezeichnet 
werden soll, wobei fiir die spezifischen Ausbeuteziffern aus demselben Grunde noch obere 


Ue ES RSE yk: 0 
a ui Usa 2020 
0 ul ee OO (2) 
0 DROP eae or 
| 0 0 0 Oe te 


‘Auch diese Matrix ist wieder singular, da sich die dritte Zeile aus der zweiten durch Multi- 
plikation mit ul/u} ergibt. Dies bedeutet mathematisch, daB sie sich nicht umkehren, also sich 
aicht D aus a, gem4B D=11j"0 errechnen 1a8t; technologisch bedeutet dies, da® man aus der 
jAnalyse des Ausgangsgases, der Ausbeute a, nicht die des Eingangsgases, des Durchsatzes 0, 
erschlieBen kann, da ja der Umsatz so geleitet wurde, daB sich im Ausgang ein festes Ver- 
1altnis von CO:H, einstellt. Durch die Umsatzmatrix 1 wird der Betriebsvorgang mathe- 
matisch eingefangen. Dabei wird der Umstand, daB der Betriebsvorgang von verschiedenen 
iHinfluBgréBen, in diesem Fall von den verschiedenen Komponenten des Durchsatzes abhangt, 
Hadurch beriicksichtigt, da& der mehrdimensionale Vektor D in die Rechnung eingeht. Auf 
pines mu} hier aber noch ausdriicklich hingewiesen werden: die einzelnen Komponenten der 
Vektoren D0 und a diirfen nicht negativ werden. Dies bedeutet im vorliegenden Fall, da die 
Matrix 11, nur anwendbar ist, solange 

(1— ul) d,— uk d;=>0, oder 2 

3 


up 


1—u} 
bleibt, da sonst, wenn etwa d, gleich Null ware, sich eine negative CO,-Menge a, ergeben wiirde. 


Aber auch bei positivem d, gilt diese Grenze fiir & da bei ihrer Unterschreitung sich gemaB 


d. 9 

ssleichung (1) die CO-Menge erhéhen wiirde, wad dies ist nicht der Zweck der Konver- 
Lierung ; Gase, die einen derart niedrigenCO-Gehalt aufweisen, wird man im allgemeinen in 
vine andere Betriebsstufe leiten. Bei praktischen Rechnungen stért diese Giiltigkeitsgrenze 
Her Umsatzmatrix nicht, doch ist immerhin bei formalen Rechnungen eine gewisse Vorsicht 
im Platze. Wenn oben gesagt wurde, da® durch die Umsatzmatrix | der Betriebsvorgang 
mmathematisch dargestellt wird, so ist dies nicht gleichbedeutend mit der zahlenmaBigen Dar- 
Mtellung eines selbstandigen Reaktionsablaufes. Die Reaktionsbedingungen sind zwar 
‘mplizit in den Matrizenelementen enthalten; daneben aber auch die durch die Apparatur 
nd deren Bedienung bewirkte Beeinflussung des Reaktionsablaufes. Zuweilen ist es not- 
vendig, die Ausbeutegleichung 


| = UD 
1och durch einen von der Durchsatzhéhe unabhiangigen Vektor Q) zu erweitern und 
a=Ud—at 


a schreiben, wo t die Zeit und a bzw. D die Ausbeute bzw. den Durchsatz fiir diese Zeit be- 
nedeuten, wahrend der Vektor a, die vom Durchsatz unabhangigen Ausbeuteverluste in der 
Zeiteinheit zum Ausdruck bringt. In solchen Fallen ist also ein gewisser Mindesthbetrag fiir 
len Durchsatz erforderlich, damit die Apparatur iitberhaupt zum Fahren kommen kann, d. h. 
a8 die Komponenten von q positiv werden. Was die zahlenmaSige Ermittlung der Matrizen- 
slemente anlangt, so sei hier nur kurz angedeutet, da® man dabei neben theoretischen Er- 
wagungen ausgiebig von der Auswertung der Betriebsergebnisse auf statistischer Grundlage 
Gebrauch macht. (Naheres hieriiber siehe Ziff. 5.) 

Was nun die Umsatzmatrix der beiden Kompressionsstufen anlangt, so ist zu bedenken, 
JaB sich bei der Kompression die Zusammensetzung des Gases nicht andert. Man kann 
jaher, wenn man die Verluste der Kompressionsstufe auf 30 at durch spezifische Ausbeute- 


“iffexn gemaB 


0<VW= w= wv=—uj=u=u <1 


‘um Ausdruck bringt, die Umsatzmatrix der Kompression durch 


| aC = Us, & (3) 


122 Pichler: Anwendung der Matrizenrechnung auf betriebswirtschaftliche Aufgaben. 


als Skalarmatrix (€ bedeute die Einheitsmatrix) anschreiben. Entsprechend gilt fiir die Kom 
pressionsstufe von 30 bis 300 at | 

YU, =u. (4 
Fir die Kohlensaureentfernung aus dem Konvertgas ist das Auswaschen mit Wasser die ein}! 
fachste und groBtechnisch am haufigsten angewandte Form. Der Kohlensiuregehalt dep 
gewaschenen Gases betragt 1—2%; bei dieser Wasserwische werden auch von den ibngay 
Komponenten gewisse, von ihrer Léslichkeit abhangige Anteile mit entfernt. Bezeichney) 
man den CO,-Gehalt im gewaschenen Gas a mit ¢ bezogen auf das gewaschene Gas ohne COf/ 
Tae bezogen auf das gesamte Gas einschlieBlich CO, }, fihrt man ferner fiir “y 
Komponenten CO, H,, N,, CH, die spezifischen Ausbeuteziffern u8, uf, uz, uZ ein, so 14Bt sick 
der Betriebsvorgang bei der Druckwasserwasche mit praktisch hinreichender Genauigkei 
durch die Umsatzmatrix 


[also mit 


RD cue cus cul cus 
0 tery) 0 0 
0 0 uw 0 0 (5 
Palas liens eae’ 
0 0 0 0 uz 


darstellen. 

Bei der Entspannung des Waschwassers werden die geléste Kohlensaure und ihre Begleit 
komponenten wieder frei. Das so als ,,Nebenprodukt“ anfallende Entspannungsgas kan 
unter Umstinden fiir andere Zwecke (beispielsweise als Schutzgas) Verwendung finden. Ki 
soll hier als ,,Fertiggas‘ 

fs = (fi 


bezeichnet werden, um zum Ausdruck zu bringen, da8 es in das betrachtete System nicht mei 
eingesetzt wird, von ihm aus also als Fertigprodukt anzusehen ist. Der Zusammenhang) 
zwischen f, und D0, kann durch 


fs = [ 03 —cug —cug —cuj —cuz Ds = Bz Ds 
Goud slinds yO Hao 
We ee Se (5s 


0 0 0 v3 0 
| 0 0 0 0 v3 
dargestellt werden, wo 0 < uj + v3 <1 ist (mit uf? = 0). 

Nach der Kompression des Gases auf 300 at (Ul, = u, ©) erfolgt in der Kupferwasche die 
Entfernung des Kohlenoxyds durch Waschen mit ammoniakalischer Kupfersalzlésung, di¢ 
durch Entspannen und Entgasen wieder regeneriert wird, wobei die hierbei frei gewordene# 
Gase (,,CO-Riickgas“‘) wieder in die Konvertierung zuriickgefiihrt werden, soweit sie nick 
anderweitig Verwendung finden. Auch in diesem Fall werden zur Darstellung des Betriebs 
vorganges zwei Matrizen benétigt, namlich eine, welche den Ausbeutevektor f; des Frisch 
gases gemaB 


fs aes &; Ds 
und eine, welche den Ausbeutevektor a; des CO-Riickgases gemaB 
ds = Us De 


in Abhangigkeit vom Durchsatz D; zum Ausdruck bringt. (Die Wahl der Bezeichnungen ist mit 
Riicksicht darauf erfolgt, daB das CO-Riickgas u. U. in das System wieder eingesetzt werder 
kann, ein Fall, der Ziff. 3b behandelt ist). Beide Matrizen lassen sich mit praktisch hin 
reichender Genauigkeit als Diagonalmatrizen gemaB 


5 te Oe Omen | und UU; = {us 0) 10 20580 
O78 OF 0, 0 OF 02 20. 0 70 
Or 0 ee OO 10. 0% 21800 (6} 
00 208270 | OS Ur, 20 
0 0 0) 0) | OTRO ON ie 
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darstellen. Man kann, was gelegentlich zweckmaBig ist, auch beide Matrizen zusammenfassen 
jand schreiben: 

(F) — (3°) Ds 

fs) \&s 


= 3. Kopplung von Betrieben. Im vorhergehenden Abschnitt ist gezeigt worden, wie sich 
‘Betriebsvorginge in der chemischen Industrie durch Matrizen darstellen lassen. Diese Dar- 
(tellungsméglichkeit allein wiirde jedoch die Anwendung des Matrizenkalkiils wohl kaum 
vechtfertigen. Seine Brauchbarkeit zeigt sich in vollem Umfang erst bei der Kopplung von 


vinzelnen Betrieben, die in folgendem fiir verschiedene Falle betrachtet werden soll. 
i 


3a. Die unverzweigte Betriebskette. Betrachtet man zunichst eine F olge von 


4 3 , 3 d,=t 
jnintereinandergeschalteten Teilbetrieben, eine »unverzweigte Betriebskette“‘, so ; 
jsann man aus den Umsatzmatrizen der Teilbetriebe sehr leicht auf die Umsatz- 
matrix der ganzen Betriebskette schlieBen. Aus Abb. 2 fiir eine aus drei Teil- a 
etrieben bestehende Kette liest man ohne weiteres ab: 
| a1=b2 
| a, =U, d,, | 
a, = U,d, = Uy a; = Uy Udy, (7 
as = Uy d, = Uy a, = U, U, Uy Oy. | ee 
4 2=03 
sBezeichnet man den Durchsatzvektor 0, des ersten Betriebes mit 1 (,,Roh- 
oroduktvektor“), den Ausbeutevektor a, des letzten Betriebs mit f (,,Fertig- Ny 
oe) und fiithrt die Abkirzung 
{ U,, = U; U, Wy (8) a5-f 
bin, so erhalt man Saver 
4 Betriebskette. 
, j=U,yt. (9) 


“Man erkennt aus den Gleichungen (7) bzw. (8), daf§ der Hintereinanderschaltung der Teil- 
vetriebe mathematisch die Multiplikation ihrer Umsatzmatrizen entspricht. Nun ist bei der 
Matrizenmultiplikation die Reihenfolge der Faktoren im allgemeinen nicht vertauschbar. 
‘Dieser dem Nichtmathematiker zuweilen schwer vorstellbare Satz bedeutet aber technolo- 
zisch, daB die Reihenfolge der in den einzelnen Teilbetrieben stattfindenden technologischen 
Prozesse im allgemeinen nicht vertauscht werden darf. Dies ist aber ohne weiteres klar. Um- 
izekehrt 148t sich aus der Nichtvertauschbarkeit technologischer Prozesse schlieBen, da8 zur 
-ationellen mathematischen Darstellung solcher verkoppelten Prozesse ein Kalkiil mit nicht- 
vertauschbaren Verkniipfungen Anwendung finden mu. Der Matrizenkalkiil erscheint also 
teradezu pradestiniert fiir solche Aufgaben, zumal durch die Verwendung mehrdimensionaler 
‘Vektoren beliebig viele EinfluBgréBen beriicksichtigt werden kénnen. 

Die Verallgemeinerung der Gleichungen (8) und (9) liegt auf der Hand; dabei mége auch 
sleich der Fall mit einbezogen werden, daB in jedem Betrieb ein Fertigprodukt f; anfalle, 
Jessen Durchsatzabhangigkeit durch 


dargestellt sei, wahrend die an die nachfolgenden Betriebe weitergeleiteten Ausbeuten Q; 
wieder durch 


fi = Bid; (@@ =1,2,--+ m) 


i 
i 
! 
t 
l 
i 


a; = U;0; (i= 1,2, +m) 
gegeben seien. Unter Beriicksichtigung, daB 

| Dp (a (Gia) 

‘st, folgt dann mit 0, = 1 

) a;=U;Ue1--: Uh = Yat, (Geel, 2, mn) 

De Oy Ui ty a (i = 2, 3,+-+ m) (10) 
| eS B0,= Bitar a= Viti: J 

| hh = Bi Dy = eon vy 

or 
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oder in anderer Schreibweise: 


d=(d,/= | i a = (or (et ae Bi |t- 
Dy Un A, Ue fe a1 | 
Das Un—1,1} Om Ua Ve and 


Fs kann nun vorkommen, da8 nicht nur dem Betrieb 1 ein Rohprodukt 1, zugefiihrt wird, 
sondern da® auch den Betrieben k (k > 1) neben der Ausbeute ,_, ein von auBerhalb des 
Systems kommendes Rohprodukt 1, zur Weiterverarbeitung zugeleitet wird. Es gelten dann 


‘ folgende Bedingungen: 


a;= U;d,; (eer eI : 

ieee Gea ea (Umsatzgleichungen) | ai 
- er ae ti (i = 2,3, -* «m) } (Kopphungegleichungen)| 

i= ty 


Die Lésung 1af%t sich unmittelbar anschreiben, wenn man beriicksichtigt, daB wegen 
10;0; = Usa +t) = Wont, bzw. B,d; = Biaat Bil 


die aus t,(k>1) stammenden Ausbeuten und Fertigprodukte einfach superponiert |, 
werden kénnen. Fiir den Einsatz von 1; allein folgt aber analog wie fiir 1): 


De tee 
Di lee te (k <4 Sm), 
Op Uiaete (ht em 
(2 aon Ue (cies): 
wo 
Ue Uae, (ham) Ge A 
Bi, = Bi tla. (k <i Sm)y'-und Bi, = Bz 
gesetzt ist. Setzt man noch fest, daB® U;, fir i< k gleich Null (d. h. gleich einer n-reihigen, 
quadratischen Nullmatrix) sein soll, so 1a8t sich schreiben: 
d= D, ad & 0 * 0 Y | = (6;, ©+U,4,,) tv mit t= ( ty (On= i 
. . . - 0 . “i Oik = (I) 
[ Din ees, WELD A : & Un ae fiir tz k) 9 
a= | = a 0 oc. rises 
a U 0 t 
2 vay Maa Ff 2 (12) 
Om Die Urs z je Vm 
{= in = Bau 0 co Ty =(%B;,) YT. 
fe Bo Be * ‘ Vp 


ie Bint ins 7 Se is 135 


Man kann diese Gleichungen noch in eine einzige zusammenfassen, wenn man 


a ei = Oi. © + Uy St (Us, = Baz 0 fire <Ak)} 
a in; Opie ae 
j Bip =0, i+hk) (yo aa 

(i, k =1,2,+++m) | 


setzt, und hat dann den ,,Betriebsablauf*! als Folgeerscheinung des Rohprodukteinsatzes 
in einer einzigen Matrix zum Ausdruck gebracht. Die praktische Bedeutung der Darstel- 
lungen (12) bzw. (12a) liegt darin, daB sie es ohne Umweg itber samtliche Zwischenstufen 
erméglichen, sofort die Wirkung der Anderung einer Rohproduktkomponente auf eine be- 


| XXI. Band 1953. Pichler: Anwendung der Matrizenrechnung auf betriebswirtschaftliche Aufgaben. 125 


jliebige Komponente von Dd, a oder { angeben zu kénnen, eine Forderung, die sehr haufig an 
die Betriebsleitung herantritt. 

Wendet man die Gleichungen (12) auf die unverzweigte Betriebskette der Wasserstoff- 
erzeugung von der Konvertierung bis zur Kupferwasche an und zwar unter der Annahme, 
) daB das CO-Riickgas nicht in die Konvertierung zuriickgefiihrt wird, so erhalt man unter Be- 
jnutzung der Gleichungen (2) bis (6) fiir an die in Tabelle 1 zusammengestellte 


ik 
| Matrix. Wenn sich auch in diesem Sonderfall diese Matrix auf einfache Weise in allgemeiner 


| Form darstellen 1aBt, so empfiehlt es sich doch bei praktischen Anwendungen, die zahlen- 
imafige Darstellung von Q unmittelbar aus der zahlenmifigen Darstellung von I; und %; 
jzu berechnen. Ein Beispiel fiir eine solche Darstellung gibt Tabelle 2, bei der je- 
i} doch die Kompression auf 30 at mit der Druckwasserwische und die Kompression von 30 auf 
) 300 at mit der Kupferwasche zusammengefaBt sind. Da die Umsatzmatrizen fiir die Kom- 
ipressionsstufen mit sehr guter Annaherung durch Einheitsmatrizen darstellbar sind, bleibt 
; das typische Bild bei dieser Darstellung erhalten. 


7 


a ¢ 3b. Die Betriebskette mit Riickfiihrung. Es wurde bereits (Ziff.2) dar- 
i auf hingewiesen, daB das in der Kupferwasche anfallende CO-Riickgas, 
soweit es nicht anderweitig Verwendung findet, wieder in die Kon- 
d= tras vertierung zuriickgefiihrt wird. Solche Falle der ,,Riickfiihrung“ sind 
besonders fiir die chemische Industrie typisch. 
ie In der Abb. 3 ist eine solche ,, Betriebs- t; 
kette mit Riickfiihrung“ fiir den Fall é 8; oe 
ai=bz gezeichnet, da® die Ausbeute a, = UD, By, i 
des dritten Teilbetriebes in den ersten of 
zuriickgefithrt wird; das Fertigprodukt f; [iat h pa[us 
sei durch 


on 2 


' 


i =d 

F we fs a Bs Ds 

| ts & gegeben. Man erhalt zunachst wieder wie i 
bei der geraden Betriebskette oe fi A 

tm, < wit 

{ re i= U; Uy Uy, Dy a Us; dy : Nay <i 

am fs Ui, B 


Nun gilt weiter: 
Abb. 3. Betriebskette AbboA, Zallische Auff 4 
Be eet ae 2 ts .4. Zyklische Auffassung der 
mit Rickfibrung. D, = 1 = as est Ae Usy Dy Betriebskette mit Riickfiihrung. 


oder mit der Hinheitsmatrix ©: 

(© — U1) Dd, =v. 

| Daraus folgt unter der Voraussetzung, dafs (© — U,,) nicht singular ist, 
b, = (E— Uy). 

Fir das Fertigprodukt f, folgt nun 

| fo = Bo By = By dy = By Way By = By My (E— Uy) 
| oder schlieBlich mit B,, = ¥, Uy; 

fs = Bai (E — Uy)" t. 

| Anch bei der Betriebskette mit Riickfiihrung liegt die Verallgemeinerung auf m Teilbetriebe 


auf der Hand. An Stelle der Bedingungsgleichungen (11) (Ziff.3a) treten jetzt unter zy- 
| klischer Auffassung der Folge 1,2---m (Abb. 4) folgende Gleichungen: 


0; = U;D;. | 
, f= %B,d,, J 
| D;=a;, +1t;, (Kopplungsgleichungen) | 
(Gene (oe Lm): 


. (Umsatzgleichungen) | 
(lla) 


, Die Auflésung nach },, a; und f; erfolgt genau so, wie oben fiir m = 3 gezeigt wurde, und man 
-erhalt wieder unter Benutzung der Abkiirzungen ;, und &;;, (Ziff. 3a), zunachst fiir den 
| Fall, daB nur an der Stelle 1 ein Rohprodukt 1, eingesetzt wird (also 1; = 0 fiir i > 1), und 
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Tabelle 1. Matrix (= Hh fiir eine unverzweigte} 
“wu 
yy Pee erg 
Se ore 1 
(ui 1—u} —ul 0 0 
Ry ut ui 0 : 
1 ul 0 
Uy, 1 p Re 2 a 0. 
0 0 0 0 ul 
i= —us 0 0 
i. e a! i 0 0 
: 0 0 
Us, 2 0 u} Us A 
0 0 0 ut 0 
1 
[0 0 0 0 u} 
Uy (0 c(us us + u3 ul) cus ra Oe u}) cu} - a ui 
aa Fa 
si 3 0 uh ut us ul 0 0 
: Ail 
0 0 0 us ut 0 
0 0 0 0 u3 ut 
Uy Uy (0 c(us ui + uf u}) c(uz ul + u3 u}) OU Wt cui ui 
0 u3 ua ee a ; 4 
Uap 4 0 wes DET 
0 0 0 u3 u} 0 
0 0 0 0 Us te 
Wyn Wis (0 cu?(ug ud + u uf) CUE OS = es) GIR Us cu’ u ud 
0 DIC Us us us uh 0 
Wy 5 0 OS ai us us us 0 0 
0 he Op 0 
[0 0 0 0 ue u? u! 
Uy [v3 ut vi (1—u4)—c(u3 ub+-u3u}) —v? uy —e(ugus + u345) —_ ae ul —c us ul | 
0 v3 u v3U 
Bs 3 | 0 o ui , 08 ui 0 0 | 
0 0 0 Os 0 
Leu 0 0 0 ou 
By Uh (0 cv' (uz ud+u3 w}) cv} (us uh + v3 u}) cv? u3 ul cv? us u}} 
0 v3 us ud ORUR Ti = 0 0 
2b a 5 0 2 we us ve us ul 0 0 
; 0 oes ages vue ul 0 
0 0 0 0 Oe went 
unter der Voraussetzung, da®B (© — U,,;) nicht singular ist: 
Dy cs. ( aa Vas > 
D: = Usa (© — Uni) 1h (i = 25,3, >=" a) 


= (6 — Ue Gs | 
f= Ba (E— Uy) y. J 
Durch Einsetzen dieser Lésung in (lla) wird leicht ihre Richtigkeit bestatigt; dabei ist, um 
Dy = (E— Wa) 1 = Ge + = [Un (E — Un) 34+ OC] y 
nachzuweisen, die Identitat 
(E—W7I=uUu(E— WN7+E€=(EC—W73U+C (13) 


zu beachten, die aus 


(9 2) eee 


GC=u+C—w 
folgt, wenn man beide Seiten, von rechts bzw. von links mit (© — M1) multipliziert. 


Wird an der Stelle k das Rohprodukt 1, eingesetzt, so erhalt man bei dem zyklischen 
Charakter der Betriebskette mit Rickfithrung die Lésungen ganz einfach dadurch, da® man 


a 
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9. 
U;, baw. B;. 


ou, E 0720) 20, 40 | 
Ope ee 0 On 20 
ORee 0 = aA 0.65.0 
Z re ieee oan | 
ONO. 0) 0 4 ed 
—U, | 0 cuz cuz cuz cus OM cures tcus cue 
0 uz : 0 Oise Oe Or 0 
OO se sr i OF 0 S80 0 
OR Ogee Ose 0 U Os Waa 
{ Ga00 8 Oat 07-050 20” ut | 
Bu,u, | 0 cus cus cus cui U4, (0 BU RUM CU Ore 
O GR Oe Yo On Oo e 
i OOF GR A, OMe OM mee a0 e a0) 
1 | OmeOT Oe 3) 0 DO Or GR 
Peas 0. 0 un! O10) 0) ot | 


eqooe o| ooorFo 
ooeoo;ooreoeo 
==) a) 


— i — i — 


iy Uy (0 cuz us cujuz cugu? cu?u? Ocufus cufus cupuz cuzue au 07050 0| 
| 0 uu 0 0 0 Oujus 09 0 0 5 0 uso 0 Oj 
0 0 usu 0 0 0 0 usue 0 0 5 0 0 u0 0 
be One uint. 0 Th Er eae a 5 lo oo uwol 
tO 0) 0 ul ue Pees 0 ula us} |[0 0 0 0 us| 
vi —cu3 —cu —cu} —cu $ —cu3 —cu3 —cu} —cu3 
0 Ue aU 0 0 ga 0 0 
0 OND ve 0 0 0 v3 0 0 
org "3 0 Oo oe 6 
SAC Ok ee ae Oo ee OIE 92 
0 es cvus cvjui cvjue|| u, [0 cvjus cvjus cvius cvju3 0 0 oF Oy wo OW 
0 v 0 0 0 Ores 0) 0 0 0 0 eh De ie) 
0 Th Oe OD Xb Os TOF fetus <0 0 v3, 0 Y O oO 0 
Geol 0. ius” 0 0 0 0 vfu2 0 Oe lo 0 0 v30 | 
0 0 0 0 v8u3 Oa 0-0 0 ub? 0 0 [0 0 0 0 o3| 
iden Zeiger 1 mit dem Zeiger k vertauscht. Die Matrizen U;;, und %;,, die gemaB 
eal ees Uy, Gale 12, 
Bn = Vi i (i + fy), (14) 
Vin = Ve 


‘Gegensatz zur geraden Betriebskette (s. 
definiert. Fiihrt man weiter noch die Abkiirzung 


Uke = (© — Uy_1,5)7* 

ein, so erhalt man fiir den Rohprodukteinsatz an der Stelle k die Lésung 
Veet (© + Una Var) Ve ml 
oe tal ara Mar Ue (i +k 
p= al, Man Te 

i= Vin Ver Tk: 

Hier ist wieder die Identitat (13) zu beachten, aus der 


® = (© — U4) 7 = €4 Uae (E— Vana) 1 = C+ Uns. Un = E+ Wn Ue+,x (18a) 


‘definiert sind, sind nun ebenfalls zyklisch aufzufassen, d.h. es ist Wy. = = I,, zu setzen. Im 
Ziff. 3a) sind nunmebhr U;;, und %,, auch fir i <k 


(14a) 


pate natin. 
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folgt. Durch Superposition erhalt man schlieBlich fiir den allgemeinen Fall, bei dem an jeder 
Stelle Rohprodukte eingesetzt werden, 


d= D ae Un Une ie ‘ US Ves Y aa (On c + Ue k Was) t, 
De Wee ie ie Lien ees r,| (© = n-reihige Einheitsmatrix) 
: . . . . : (0,;, =1 fir t—& 
Os Meio ee ie ie : ae | Uy, =0 , t+hk) 
C= | 07) == | Ya Un Ure Tye : Ung. Nes Ly || = (Win Unx) Y, 
M5 User Uy Wen Use + Usm Unm | | Ye (15) 
: a ee z 5 | 
a | tea Na us Wher r ie mm } Tk | 
j= hh =a Bi Wn Bio Us, U One Ue TY, ra (Biz ex) i 
le Bor Uyy Bae Ugg + Bom Umm | | Te 
: soe fad a : | 
ies Bint Ya Ome Ue 2 Soe, tie Tr | 
Auch hier 1aBt sich, wie oben, zusammenfassend schreiben 
(= Ate bin € + Usa,e Wen Leen te 
a | Nis Une (15a) 
ies Bie Tee | 


Tabelle 3 zeigt ein Zahlenbeispiel fiir eine Betriebskette mit Riickfiihrung, dem dieselben 
Umsatzmatrizen der Teilbetricbe wie in Tabelle 2 zugrunde gelegt sind. 

Der mathematischen Voraussetzung, daB (© — U,,;) nicht singular ist, einer Voraus- 
setzung, zu der noch die Bedingung hinzukommt, daB® die Komponenten von (€ — UL,,;)-!1 
nicht negativ werden diirfen, solange es die von t nicht sind, entsprechen technologische Be- 
dingungen fiir die Durchfithrung des Verfahrens mit Riickfiihrung. Es gibt Faille, wo diese 
Voraussetzung durch besondere Ma®nahmen erzwungen werden mu. Man kann die in solchen 
Fallen auftretenden Verhiltnisse leicht ttherblicken, wenn die Matrix UU, eine Diagonal- 
matrix ist. Bezeichnet man die i-ten Komponenten von t, }, und a,, mit r;, d und a‘™, 
so gelte also etwa 

al™ =u, d™, 


Ist nun u; = 1, so hei®t dies, daB die zuriickkommende Menge a$”) gréBer als die eingesetzte 
Menge d{) oder mindestens ihr gleich ist. Es ist also keine Abnahme von d erfolgt und 
- es leuchtet ein, daB unter solchen Um- 

standen ein stationarer Kreislauf un- 


moglich ist. Die Gleichung 
d= (© — Uy) 2 


r 
u aoe eee 
oy 4 
2 lautet in diesem Fall fiir die i-te Kom- 
ah ae 
do nr 3 
: ; fm a 7 fm 


vA 
™ l—u, 
a 


d\) wird also fiir u; 21 und r;>0 

Abb. 5. MeGaaloien Be einer singuldiren gaits (© —U,, 1): entweder unendlich und (G _— ia 

a) Einschalten einer zusatzlichen Betriebsstufe. : x . p 
b) Zwischenproduktentnahme. singular, oder es wird negativ und es 
entspricht (€ — U,,,) ebenfalls nicht 
den Bedingungen, die oben als notwendig bezeichnet wurden. 

In solchen Fallen geht man z. B. in der chemischen Technik so vor, da® man entweder 
eine zusatzliche Betriebsstufe mit der Umsatzmatrix U, einschaltet(s. Abb. 5a, wo die Betriebs- 
kette durch ein einziges Kastchen dargestellt ist) und dadurch an Stelle der Matrix (© — Ua 
die Matrix (© — U, U,,) tritt, oder man fihrt einen Teil der Ausbeute A, aus dem System 
ab (Zwischenproduktentnahme, s. Abb. Sb). In diesem Fall tritt an die Stelle der Matrix 
(© — U1) die Matrix (© —AUn), wo 0 <A <1 ist. 
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4, Wirtschaftliche Probleme. Durch die Gleichung (12a) bzw. (15a), also durch 


1=QMr mit 1 =/d 
( 


j 


ist der durch die Durchsatze 0, die Ausbeuten a und die Fertigprodukte | gekennzeichnete 
Betriebsablauf { als Funktion der eingesetzten Rohprodukte 1 dargestellt. Man bekommt nun 
einen bequemen Zugang zur Behandlung wirtschaftlicher Probleme, wenn man den oben ein- 
gefithrten Begriff des Betriebsablaufes etwas erweitert. Bei wirtschaftlichen Problemen spielt 
der Betriebsaufwand eine wesentliche Rolle. Es ist zweckmiaBig, die verschiedenen Aufwand- 
arten, z. B. Rohstoffe, Hilfsstoffe, Energien (soweit diese von aufen bezogen und nicht in 
dem betrachteten Betriebssystem selbst erzeugt werden), Lohnstunden, kurz alle von auBen 
an das Betriebssystem gelieferten Leistungen zu einem erweiterten Rohproduktvektor t* zu- 
sammenzufassen. Auch der oben eingefihrte Begriff des Durchsatzes soll etwas abgewandelt 
werden. Dabei werde angenommen, daf sich die Leistung eines Teilbetriebes durch einen 
oder mehrere skalare Parameter D,, die Durchsatze der einzelnen Teilbetriebe, eindeutig 
kennzeichnen laBt, und es mégen diese Durchsitze D, wieder zu einem, den Ablauf des ganzen 
) Betriebssystems kennzeichnenden ,,Durchsatzvektor“ D) zusammengefaBt werden. SchlieBlich 
;mégen die Mengen der verschiedenen Fertigprodukte durch die skalaren Gréf8en F, gekenn- 
izeichnet und zu einem Vektor zusammengefaft werden. FaBt man die so eingefiihrten Vek- 


(toren gema} 
(* = /r* 
() 
j 


zu dem ,,Betriebsvektor“ [* zusammen, so enthalten dessen Komponenten die fiir die betriebs- 
iwittschaftliche Beurteilung des vorliegenden Betriebssystemes notwendigen GréBen. Wahrend 
jmun im vorhergehenden Abschnitt der Betriebsablauf | als Funktion { = Ot des Rohprodukt- 
vektors dargestellt wurde, ist es fiir wirtschaftliche Betrachtungen zweckmaBiger, den Be- 
itriebsvektor in Abhangigkeit von den Fertigprodukten darzustellen. Dabei ist jedoch zu be- 
achten, da unter Umstianden fertige Verkaufsprodukte zwangslaufig als Nebenprodukte 
\(,.Kuppelprodukte‘‘) anfallen und daher nicht als unabhangige Veranderliche in der gesuchten 
(Darstellung erscheinen kénnen. Die Fertigprodukte, die nach Ausscheidung dieser Kuppel- 


jprodukte iibrig bleiben, mégen zu einem Vektor f zusammengefaBt werden. In vielen Fallen 
(wird es jedoch auch dann nicht méglich sein, den Betriebsvektor [* eindeutig als Funktion 


ivon f darzustellen, man wird vielmehr zur Erfassung von betrieblichen Nebenbedingungen, 


‘die unabhangig von der durch f ausgedriickten Produktionshohe sind, weitere Betriebsparameter 
theranziehen miissen, die zu einem vektoriellen Parameter tt zusammengefaft seien. Ist nun 
reine lineare Darstellung méglich, so 1aBt sich diese in der Form 


(*¥= a ea B* 


nN (16) 


schreiben. « 

‘Der aus den das Produktionsprogramm bestimmenden Fertigprodukten J und den betrieb- 
lichen Nebenbedingungen i zusammengesetzte Vektor ) mige als ,,Produktionsvektor“ be- 
‘eichnet werden. Dann 1aBt sich Gleichung (16), wie folgt, in Worte fassen: 

| 4 ; : 
| ,,Betriebsvektor = Strukturmatrix - Produktionsvektor“ , 

wo die Matrix %* als ,,Strukturmatrix“ bezeichnet werde, da sie einen weitgehenden EKin- 
iblick in das betrachtete Betriebssystem vermittelt. Ein Beispiel fiir eine solche Struktur- 
‘matrix zeigt Tabelle4, wobei gemaB Abb.6 ein aus acht Teilbetrieben bestehendes Betriebs- 


system zugrunde gelegt ist. Dabei wurde f = j angenommen. Unter den Nebenbedingungen n 
ifindet sich dabei auch die Zeit t, deren Einfiihrung praktisch wohl immer notwendig sein wird, 
m etwa die in der Betriebswirtschaft haufig als ,,fixe Kosten“ (z. B. Amortisation) bezeich- 
neten Teile des Aufwandes oder die von der Produktionshéhe unabhangigen Verluste zum 
‘Ausdruck zu bringen. (Man vergleiche damit auch das Glied a) tin der Gleichung a= UL D—(, t, 
Ziff. 2. Es sei an dieser Stelle lediglich darauf hingewiesen, daB man bei der Aufstellung 
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Tabelle 4. Strukturmatrix $* fiir ein Betriebssystem mit 8 Betrieben. 


\ 


Dy« | fv n’ 
DM/Einheit Fertigprodukte nate 
Einheit °C-h 
Rohstoff ] t —50,— 1,26 0,6 
Rohstoff 2 t —30,— 0,72 0,8 1,6 
c* | Strom 10°kWh | —25, — 1,6 i 1,3 
Dampf t — 4,— 2,0 QBS 1,5 —0,4 
Lohne h — 1,30 7,0 Oy, 6,5 
Sonstiger Aufwand | DM —1 38 25 24 
Durchsatz D, 103Nm? Pees ] 
Ds 103Nm? 0,9 I 2 
| D, 10°Nm4 3 
D, 10?Nm? 
Dd D; 10°Nm$ 
D, t 
De t 
D, t 
Fertigprodukt F, DM 245,— 
j F, DM 170,— ; 
IPG, DM 160,— 


bie=e.(—50, 30, —25, —4, —1,30, 1,0, 04. 0,.00.20, 0, 05,0, Ao mide eo) 
SS ee 

— bi, = bi 

g’ = b/, B* = (65,30, 42,94, 41,05, 1,60, —998,50). 


einer Strukturmatrix im allgemeinen von entsprechenden Matrizen der Teilbetriebe auszu- 
gehen hat und dabei von dem Matrizenkalkiil in ahnlicher Weise Gebrauch machen kann, 
wie es im vorhergehenden Abschnitt gezeigt wurde. Nahere Ausfiithrungen dariiber sollen 
einer besonderen Mitteilung vorbehalten bleiben; ein 


by fo einfaches Beispiel wird jedoch weiter unten noch ge- 
bracht. 
Betrieb2Z 
0;=Lp Bei dieser Darstellung empfiehlt es sich, die ver- 
L schiedenen Aufwandarten, soweit irgend méglich, 
zunachst mengenmaBig, also etwa in t, kWh, Lohn- 
stunden usw. zu bringen, damit, wie jetzt gleich ge- 


zeigt werden soll, der Einflu® der Preise auf die 
Aufwandarten besonders herausgestellt werden kann. 
Dazu fiithrt man fiir den Aufwand 1* eine geldmaBige 
Bewertung durch den Vektor b{s ein. Dann ist der 


Betrieb ¥ 
Dy=Ly 
Yy 


Ze Z. 


Rebrieh 6 Reriee7. Bekriebé gesamte Betriebsaufwand in Geld ausgedriickt durch 
Os-F; Dy-k = 
ot 712 Dy=ls Doers 
lA b fa gegeben, und genau so ]aBt sich der durch den Ver- 
Abb. 6. FlicBbild fir das Beispiel in Tabelle 4. kauf der Fertigprodukte f erzielte Erlés durch 
Es bedeuten: he bf 
R, und R,: Rohprodukte me if. 
D, bis Dg: Durchsatze 17 + “ “ Be og 1 
Psst ay ates ausdriicken. Das wirtschaftliche Ergebnis ist damit 
F, bis F,: Fertigprodukte. durch 
G = F— R* 


gegeben. Fihrt man nun cinen Bewertungsvektor 

bjs = (— bye, 0--- 0, bi), 
den ,,Preisvektor“ ein, der zwischen —}{« und 6, genau so viele Nullen enthalt wie der Vek- 
tor D Komponenten hat, so wird schlieBlich durch 


C= bye * = be B* p = Q'p (17) 
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das wirtschaftliche Ergebnis einmal als skalares Produkt aus Preisvektor bjs und Betriebs- 
vektor [*, und auBerdem auch als skalares Produkt aus einem Vektor g’ und dem Produktions- 
vektor ) dargestellt. Der Preisvektor b/. vermittelt ein Bild tber die von ,,aufen“ gegebenen 
wirtschaftlichen Bedingungen (Preise, Lohnsatze usw.); die Bedeutung des Vektors g’ wird 
oaxrkannt, wenn man 


qG’ = bie B* = (81, Bo, ° °° Bm) 


ae é 
| wl Ae: 
Pl 


setzt und damit 
G=g'p = 2 Be Pi 
=) 
schreibt. Andert sich die Produktion um die »,Produktionsspitze*+ AP,, so andert sich das 
wirtschaftliche Ergebnis um den ,,Spitzengewinn“: AG und es gilt: 
AG. a Sk A P,, 5 


Jaraus ergibt sich die Bedeutung der Komponenten g;,, des Vektors g’ fiir viele wirtschaftliche 
Jntersuchungen (s. auch Tabelle 4). 

Zur Erlauterung soll ein einfach gehaltenes Beispiel der Kostenermittlung fiir Dampf und 
Strom dienen. Dabei werde von einer aus einer Kesselanlage, Gegendruckturbinen und 


Mondensationsturbinen bestehen- 

len Energiezentrale ausgegangen, ee D, 
Iche ein Industriewerk mit Hoch- ihe 

ac Ry, RoR =Autwond 

uckdampf, Gegendruckdampf und Dy, Bh, 0, =Durchsiize An 


i 


Strom zu versorgen hat. Die ge- = 
ahlten Bezeichnungen fir die ae Arg 
Jurchsatze, Ausbeuten und,, Fertig- Gegendruck- eae EAT ASE 
ue Dp ag 2D; 
Samtliche Aufwandarten sollen (Fa) (fs) 
Azz Ase A33 


rodukte*‘ gehen aus Abb. 7 hervor. 
yereits geldmaBig zusammengefabt ; 
ein und mit R bezeichnet werden. : | 1 


ue : i : 5 fj 
Die Fertigprodukte F, und dex Auf Abgaten: — Hochdruck- Gegendrutk- Strom 
vand R seien zunadchst, wie dies (ferkgprodukt) dampt damp 

ei solchen Problemen meist der Abb. 7. Schema einer Energiezentrale. 


Mall ist (vgl. auch die Bemerkung 
jur Aufstellung der Matrix 5$*), in Abhangigkeit von den Durchsiatzen D;, gegeben und zwar 
vurch folgende, stark vereinfachte Matrix: 


Bl 4 Aqg— Ag = Ay | =| 1 W3 (D, 

F, An, 0 Ws 0 D, (18) 
F; Ay + Ags 0 1 1 D; 

R R,+ Rk, +R, Ty i, Ts 


Es werde nun zunachst auf die im kaufmannischen Rechnungswesen iibliche Art die Kosten- 
srmittlung fiir Hochdruckdampf, Gegendruckdampf und Strom gezeigt. Man geht dabei 
Rep oe verse, wie folgt, vor: 
| Kesselanlage: Aufwand: R, =7r, Dj = Ay, 
omit Hochdruckdampfkosten (d.h. Aufwand je t Dampf): r,. 
Kondensationsturbinen: 

direkter Aufwand: R, = r; D; = rz As; ; 

indirekter Aufwand! im Dampf A,,: r, Ay = 1, Wz Agg 
Gesamtaufwand: (rz + 1, W5) Ags 


omit Stromkosten: rz +1, U3. 


1 (Es ist A,, = w3D, = w3Ag3). 
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Gegendruckturbinen: 


direkter Aufwand: Re =", Dy = 1 Aa 
indirekter Aufwand!im Hochdruckdampf: 7, Ay, = 7 2 Age» 
Gutschrift fiir Stromabgabe: — (rg + 1 W3) Age » 
Gesamtaufwand: [r,(w, — ws) + (T. — rs) ] Age 


oder, da A, = w,D, = W,Azo, 


Gesamtaufwand = In (1 3 | ee Apacs 
« 2 


somit Kosten fiir Gegendruckdampf: r, ( = i) + Sap ae 
2 
Man itberzeugt sich unter Benutzung von (18) leicht, daB 


R=1,D,-+r,D,+1%D, =n + é ( ea = epee (r,w5 + 15) Fy (19) 
2 . 


gilt, womit der Aufwand auf die , Fertigprodukte* F,, F, und F, aufgeteilt und somit deren 
Kosten ermittelt sind. Wesentlich bei dieser Kalkulation ist, da® der von den Gegendruck- 
turbinen gelieferte Strom zu den sich in der Kalkulation der Kondensationsturbinen er- 
gebenden Stromkosten (r3 +-T Ws) gutgeschrieben wird. Dieser Weg ist jedoch keineswegs der 


allgemein itibliche; es wird vielfach auch so verfahren, daB man den Gegendampf mit 2 k, DM 


je t gutschreibt, wo k, =r, die Hochdruckdampfkosten (in DM/t) und i, bzw. 1, die Enthalpien 
von Gegendruck- bzw. Hochdruckdampf bedeuten. Man spricht dann von einer kalorischen 
Bewertung des Dampfes (im Gegensatz zur energetischen Bewertung, wie sie oben angewandt 
wurde) und bekommt fiir den Gegendruckstrom im allgemeinen wesentlich geringere Kosten 
als fiir den Kondensationsstrom. Diese beiden Kalkulationsméglichkeiten wurden schon viel 
erértert, ohne da bis jetzt wohl eine endgiiltige Einigung im Streit der Meinungen erzielt 
wurde. 

Zum Vergleich soll nun der Aufwand in Abhangigkeit von den gelieferten Mengen an 
Hochdruckdampf, Gegendruckdampf und Strom (den ,,Fertigprodukten“ F,, F, und F’,) unter 
Benutzung des Matrizenkalkiils ermittelt werden. Aus der Gleichung 


Fi) =|1 —w, —ws| |D, 

Be 0 We 0 D,; 

is 0 il 1 D, 

folgt 

We 

Dy) == \l i ae WA AEG 
1 

D, 0 ag 0 F, 
] 

D, 0 — - 1 ie 


Damit ergibt sich als zwangslaufige Folge der Gleichung (18) 


R = (1, Tp, Ts) ID, == (rj, 7, 73) Jt ze = ws] |Fy) =F, + In (2 Z| ae I dbs 
D, 0 oe ay 1 + (7, w3+ 73) F's. 
. We 
|D,| Oa. Zaseee state 


ae | "a 


in Ubereinstimmung mit (19); unter den Voraussetzungen (18) erscheint somit die energetische 
Bewertung als die gegebene. Die in diesem Beispiel erérterte Bewertungsfrage ist typisch 
fiir viele Kalkulationsprobleme und der hier aufgezeigte Weg zu ihrer Lésung fihrt in vielen 
Fallen zum Ziel, vor allem auch dann, wenn die Verhaltnisse nicht so einfach liegen, wie in 
dem behandelten Beispiel. 

Auf eine Anwendung der Gleichungen (16) und (17) sei hier noch eingegangen. Bei vielen 
wirtschaftlichen Untersuchungen ist die maximale Erzeugungsméglichkeit eines Werkes an 


1 (Esist A,, = w, D, = w, Ag). 
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‘Fertigprodukten, seine ,,Kapazitat‘, von entscheidender Bedeutung. Da nun bei einem Werk, 
idessen Einzelbetriebe stark miteinander gekoppelt sind, diese unter Umstanden fiir die Er- 
zeugung verschiedener Fertigprodukte in Anspruch genommen werden miissen, ist es in solchen 
Fallen nicht méglich, die Kapazitat fir ein Fertigprodukt unabhangig von den daneben noch 
zu erzeugenden Fertigprodukten anzugeben. Die dabei auftretenden Verhiltnisse lassen sich 
jnun mit Hilfe der Strukturmatrix auf einfache Weise iibersehen. Die Kapazitaten der Einzel- 
tbetriebe, d. h. ihre maximal méglichen Durchsatze, kénnen namlich als bekannt vorausgesetzt 
jwerden. Sie mégen zu einem ,,Kapazitatsvektor“ 
i= K, a dD, max | 
K, 2 max 
Les Dy max 


zusammengefaft sein. Andererseits sind die Durchsatze der Einzelbetriebe durch 


b= Be (7) = Bey 

1 
#gegeben, wo {8** aus den zur Bestimmung von D dienenden Zeilen der Strukturmatrix ‘S* 
tbesteht. Die Durchfithrbarkeit der durch p = (1) dargestellten Produktion ist dann an 
|\folgende ,,Kapazitatshedingung* gekniipft: 
| dD = B** p(s) f, (20) 
}wo D(<) f soviel bedeuten soll, wie 
D;< K; fiir alle z. 
$Gilt nun fiir mindestens ein i das Gleichheitszeichen, so ist fiir die Produktion ) ein Grenzfall 


jdurch den Betrieb i als Engpa® gegeben. Man kann aber von keiner Maximalproduktion 
‘sprechen, solange fiir ) keine MaBbestimmung vorgenommen wird. Sind die Kapazitaten der 


“Behandlung dieser Frage ist zu beachten, daB sich fiir die Durchsatze D; keine negativen Werte 
sergeben diirfen. Die Kapazitatsbedingung (20) ist daher durch Einfiihrung zweier Kapazitats- 


£(<)d = $** p(s) f. 


‘Mit den Zahlenwerten des Beispiels in Tabelle 4 erhalt man, wenn man £ = 0 setzt und weiter 
‘beriicksichtigt, da die Betriebe 3, 4 und 5 gleichartig und aufeinander umschaltbar sein sollen 


| G = 65,30 F, + 42,94 F, + 41,05 F, — 982,50 soll absolutes Maximum werden 
sunter den Nebenbedingungen 

! Ore) O10 04 A (<) (30 

| 0 O:0n deeta0 ll 2) 31 

| 0 yen tae Re 60 

0 [eo 0 0) Ol 180 10 

| 0 ee WE Ae Na 5 

| 0 Vas Re al aa | 10 


+Um in diesem Falle das absolute Maximum von G zu bestimmen, ist zu beachten, daB aus 
«den letzten drei Nebenbedingungen bereits 


(Dj= i < [10 


Del ihe 
| Dilinsy| Balm chh0 
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folgt. Mit f a LOlie j wird aber | D, =—=27 
5 Daa 35 
10 D,-+ D,+ Ds} 60] , 


d.h. der dafiir erforderliche Durchsatz D, liegt um 4-10? Nm*/h itber der Kapazitat von 
31-103 Nm3/h, wihrend die itbrigen Kapazitaten ausreichen. Will man die notwendige Sen- 
kung AD, =—4- 103 Nm?/h lediglich durch Anderung von F,;=10 um AF, exzwingen, so 
wird AF, = AD,/0,9 und AG = 65,30 AF, = (65,30/0,9) AD, =72,56 AD,, wahrend fiir F, baw. 
F, die entsprechenden Werte AG =(42,94/1) AD, = 42,94 AD,, bzw. = (41,05/2) AD, = 20,52 AD; 
lauten. Die geringste Abnahme von G (gegeniitber dem fiir f erreichbaren Betrag) wird daher 
durch AF, =0, AF, =0, AF, = —4/2 = — 2 erreicht. Man erhiilt also fiir [ = [10] ein ab- 
solutes Maximum G,,,,,—= 213,60 DM/h. : 

Man kann das Optimalproblem noch etwas allgemeiner fassen, wenn man beriicksichtigen 
will, daB z.B. auch die Grenzen fiir die Rohstofflieferungen 1 und die Absatzméglichkeiten 
fiir die Fertigprodukte f einen Engpa8 darstellen kénnen. Man braucht dazu nur fiir den 
Betriebsvektor [* Grenzvektoren [* und [* einzufiithren und hat damit folgende, wohl fiir die 
meisten praktisch auftretenden Falle hinreichende Fassung des wirtschaftlichen Optimal- 
problems: 

»,Es wird der Vektor p gesucht, fiir den 

C= ¢?P 


ein absolutes Maximum mit den Nebenbedingungen: 


I*(s) = B* p(s) i 
darstellt.‘¢ Als Beispiel (nach Tabelle 4) sei wieder angenommen, da die Rohstofflieferung R, 
auf 12 t/h beschrankt sei. Dann ergibt sich in ahnlicher Weise, wie oben, ein Maximum 


Gmax—= 27,84 DM/h fiir F = {6,03 
5 
9,79 


Zu beachten ist bei dieser Formulierung, daB in dem Vektor p = i weitere Neben- 
nN 

bedingungen 0 enthalten sein kénnen, die im allgemeinen nicht variiert werden diirfen. Da G 
eine in ) lineare Form darstellt, handelt es sich bei der Lésung des Optimalproblems um die 
Ermittlung eines Randextremums. Mit der Aufstellung der Strukturmatrix $8*, aus der sich 
gemaB (17) leicht der Vektor g’ =bj« $* ergibt, ist ein wichtiger Schritt zur Lisung dieser 
Aufgabe getan und weitgehend vorbereitende Rechenarbeit geleistet fiir solche Falle, wo aus 
allgemein wirtschaftlichen oder innerbetrieblichen Griinden sich die Grenzvektoren [* bzw. [* 
und damit auch der optimale Vektor ) andern. 


5. Ermittlung der Matrizen fiir die einzelnen Teilbetriebe. Bei allen Untersuchungen wurde 
bis jetzt davon ausgegangen, da sich die einzelnen Teilbetriebe eines Betriebssystems durch 
Matrizen, also durch eine Anzahl linearer Funktionen kennzeichnen lassen, und der Schwer- 
punkt der Untersuchungen lag darin, wie sich aus diesen einzelnen Teilmatrizen eine das ganze 
Betriebssystem hinreichend kennzeichnende Strukturmatrix aufbauen 14Bt. Die Anwend- 
barkeit der geschilderten Methoden hangt also davon ab, wie weit solche lineare Darstellungen 
iiberhaupt méglich sind und wie sie sich praktisch ermitteln lassen. Darauf soll jetzt noch 
etwas eingegangen werden. 

Zunachst sei bemerkt, daB die Theorie vieler technologischer Prozesse noch lange nicht 
so weit entwickelt ist, wie etwa die Theorie der elektrischen Maschinen, der Wasserkraft- 
maschinen oder Dampfturbinen, bei denen sich eine Charakteristik mit gro®er Genauigkeit 
vorausberechnen aft. Es ist daher im allgemeinen auch nicht zu erwarten, daB man allein 
auf Grund theoretischer Uberlegungen die Charakteristik eines Teilbetriebes, d.h. seine Aus- 
beute und seinen Energiebedarf in Abhangigkeit von seinem Durchsatz und sonstigen vor- 
gegebenen EinfluBgréBen mit hinreichender Genauigkeit berechnen kann. Man ist daher, 
wenn man zu einer Charakteristik kommen will, auf eine Auswertung der Betriebsergebnisse 
angewiesen. Es ist also notwendig, dafi die Betriebsergebnisse laufend gemessen, und die 
MeSergebnisse zusammengestellt und ausgewertet werden. Dies ist aber eine Voraussetzung, — 
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jdie fiir einen nach wirtschaftlichen Gesichtspunkten geleiteten Betrieb schon seit langem als 
Iselbstverstandlich gilt. Liegen nun fiir eine gréBere Anzahl Zeitabschnitte mit verschiedenen 
Durchsatzen MeBergebnisse vor, so ist es naheliegend, nach den bewahrten Methoden der 
jmathematischen Statistik nach den Zusammenhingen zwischen Ausbeuten, Energieverbrauch 
jjusw. zu suchen. Bei der praktischen Durchfiihrung hat es sich dabei als zweckmaBig erwiesen, 
ibei den Zeitabschnitten, wenn nicht besondere Vorkehrungen getroffen werden, nicht unter 
yeinen Monat herunterzugehen; liegt geniigend Zahlenmaterial vor, so ist die Auswertung nach 
‘},, Vierteljahrspunkten ‘ vorzuziehen. Man darf jedoch bei diesen Auswertungen nicht rein 
}formal vorgehen, sondern mu immer die speziellen betrieblichen Verhaltnisse im Auge be- 
Hhalten und soweit, als irgend méglich, auch theoretische Erwagungen heranziehen. Dadurch 
qwerden den vorzunehmenden Ausgleichsrechnungen vielfach Nebenbedingungen auferlegt. 
fEs kann auch zuweilen zweckmabig sein, besonders wenn mehrere unabhangige Veranderliche 
jzur Beschreibung herangezogen werden miissen, eine theoretisch faBbare Gréfe mit ihren im 
fBetrieb erreichten Werten zu vergleichen und so die besonderen Betriebskonstanten heraus- 
qzuschalen. Wenn nun aber auch die Statistik einen breiten Raum bei der Ermittlung der 
/Betriebscharakteristik einnimmt, so gibt es doch auch Falle, wo eine unmittelbare Betrachtung 
wrascher zum Ziele fiihrt. So wird es vielfach méglich sein, den Arbeitskraftebedarf durch eine 
igenaue Untersuchung der bei den verschiedenen Betriebsbedingungen notwendigen Arbeits- 
iplatze unmittelbar zu bestimmen. Gerade dieses Beispiel zeigt iibrigens deutlich, daB bei einer 
tgenauen Betriebscharakteristik auch mit Unstetigkeiten gerechnet werden mu, wenn z.B. 


egestellt werden kénnen. Wenn sich dabei zu Beginn der Untersuchungen gréiere Abweichungen 
ivon einer strammen linearen Bindung ergaben, so lag — abgesehen von MeBfehlern — die 
+ Ursache meist darin, daf zunachst irgendein Betriebsparameter nicht bericksichtigt wurde, 
.h. mathematisch gesprochen, da8 zu wenig unabhangige Verdnderliche zur Beschreibung 
des Betriebsablaufes herangezogen wurden. Diese weitgehende Méglichkeit der linearen Be- 
ischreibung liegt in der Struktur der chemischen Betriebe begriindet, derzufolge kleinere 
sDurchsatzanderungen durch ,,Tiefer‘‘- oder ,,Héherfahren“ der einzelnen Maschinen und 
<Apparate, gréBere Anderungen dagegen durch Ab- oder Zu- 
Bei slten von Maschinen und Apparaten bewirkt werden. Da- 
udurch entsteht fiir die Betriebscharakteristik eine Art Treppen- 
ukurve, die, wenn die Stufen nicht zu hoch sind, mit hinreichen- 
oder Genauigkeit durch eine Gerade ersetzt werden kann. Wenn 
sman nun bedenkt, da& die einzelnen ,,Treppenstufen“ nicht in 
Jallen Betrieben eines Betriebssystems gleichzeitig erreicht werden, 
‘so ist einleuchtend, daB fir das gesamte Betriebssystem die 
jlineare Darstellung mit wesentlich gréBerer Genauigkeit als 
‘fiir die einzelnen Teilbetriebe méglich ist, was sich praktisch te ee 
simmer wieder bestatigt hat. Es ist jedoch nicht immer méglich, 
Nauftretende Unstetigkeiten einfach zu glatten. Wenn die Beibehaltung von Unstetigkeiten 
ifiir die Weiterrechnung notwendig ist, kann man im Interesse einer glatten Rechnung haufig, 
iwie folgt, vorgehen. Es sei etwa eine Aufwandart A, z.B. der Dampfverbrauch, gemaf nach- 
ostehender Beziehung vom Durchsatz und der AuSentemperatur abhangig: 


| 


LDurchsatz 


t 
| A=aD+b,fodx fir 9<%, 
; 0 


i 
! 


=aD+b,fOdx fir 99> %, 

0 

swo A bzw. D Aufwand bzw. Durchsatz fiir den Zeitabschnitt t und @ die AuSentemperatur 

t t 

‘bedeuten. Wenn die Glieder b, [o dx baw. b, [0 dx beispielsweise den fiir die Raumheizung 
0 0 


10 
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bendtigten Dampf zum Ausdruck bringen, so wiirde b,=0 heiBen, dah beim Uberschreiten 
der Temperatur 9, keine Heizung mehr erforderlich ist. Schreibt man nun 


A=aD + }byt,+ byt, 
mit = [ ode, W=9 fr 05%, 
O == 0 fur? 8S" 05; 
und t= [ v'dx, 9” =0 fir 0 = By, 
w=) fir 0>H, 


so hat man formal die Unstetigkeit (bzw. den Knick) bei @) dadurch beseitigt, daB man an 
t 
Stelle des einen Parameters f 8 dx zwei Parameter t, und t, eingefiihrt hat. t, bedeutet dabei 


das Temperaturzeitintegral fiir den Teil des Zeitabschnittes t mit einer Temperatur < # und | 
das Entsprechende gilt fiir t,. Auf diese Weise wird erreicht, daB man mit der, allerdings um 
eine Spalte erweiterten Matrix glatt weiterrechnen kann. Soll aber aus der Strukturmatrix B* 
fiir irgend einen Produktionsfall p (p enthalt als Komponenten die Parameter t, und #,) der 
Betriebsablauf [* berechnet werden, so ist es nur notwendig, vor Beginn der eigentlichen 
Rechnung die Werte fiir die Komponenten #, und t, zu ermitteln, und die u. U. vollautomatisch 
durchgefithrte Rechnung (s. nachster Abschnitt, Ziff. 6) lauft ebenfalls glatt. ErfahrungsgemaB — 
stellt die Aufstellung der Matrizen fiir die einzelnen Teilbetriebe weitaus den gréBten Teil | 
der fiir den Aufbau der Strukturmatrix des Gesamtwerkes zu leistenden Rechenarbeit dar. — 
Es ist aber zu erwahnen, daB bereits die Teilmatrizen ein brauchbares Hilfsmittel fiir die wirt- — 
schaftliche Uberwachung der einzelnen Betriebe darstellen, wenn man sie als Mafistab bei der 
Beurteilung der laufenden Betriebsergebnisse verwendet. Auf eine Méglichkeit, auch nicht- 
lineare Zusammenhange mit dem Kalkiil noch zu behandeln, sei abschlieBend kurz hinge- 
wiesen. Es sei etwa die Vektorfunktion 


Yi =e =%% 4% = %U[x 

Bh) Xo 

ee Xn 
gegeben, ferner etwa die Funktion: 

Ym+1= 0 * Uy+ Ym Xn 5 
dann ]aBt sich schreiben 
)=/% = | 0} 1% = Ur. 
Vm+1 a a Ym Xn 


Die Zahlenrechnung geht dann so vor sich, daf§ zunichst der Vektor {), berechnet wird; damit 
ist dann die m-te Komponente des Vektors tj, und damit auch die (n +1)-te Komponente 
Ym%n des Vektors ¢ gegeben, so daB nach der Zwischenrechnung fiir die Ermittlung von yp %,_ 
auch y,,;1 errechnet werden kann. Man bleibt so, abgesehen von einer kurzen Unterbrechung, 
im Kalki, und dies ist fiir die praktische Zahlenrechnung nicht unbedeutend. 


6. Praktische Anwendungen. In den vorausgehenden Abschnitten wurde gezeigt, in welcher — 
Weise sich der Matrizenkalkiil bei der Darstellung von Betriebsablaufen bei Betriebssystemen — 
verwenden und die meist sehr umfangreichen Zahlenrechnungen auf die Grundoperationen — 
der Matrizenrechnung zuriickfiihren lassen. Das Ergebnis dieser Rechnungen ist eine Struktur- — 
matrix, in der in iibersichtlicher Weise die Struktur des Werkes beschrieben wird und die 
eine zahlenmaBige Kenntnis der Betriebszusammenhange vermittelt, die sonst, besonders bei — 
umfangreichen Werken, sich kaum in dieser Ausfihrlichkeit gewinnen 148t. Voraussetzung 
ist allerdings, daB eine lineare Darstellung ttberhaupt méglich ist. Diese Voraussetzung ist — 
jedoch in vielen Fallen praktisch mit hinreichender Genauigkeit erfillt; auch 148t sich haufig — 
Linearisierung durch einfache Mittel erreichen. Der Hauptvorteil fiir die praktischen An- — 
wendungen der Matrizen auf Aufgaben wirtschaftlicher und technologischer Art besteht nun 
darin, da®B mit der Aufstellung einer Matrix fiir viele solche Aufgaben der gréBte Teil der 
Rechenarbeit bereits vorbereitend erledigt und das Ergebnis dieser vorbereitenden Arbeit in 
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tier Matrix gespeichert ist. Die Durchrechnung eines speziellen Falles besteht dann vielfach 
aur noch in der Multiplikation einer Matrix mit einem Vektor, und diese kann in kiirzester 
(Zeit und ohne besondere mathematische Kenntnisse durchgefiihrt werden. Beides ist in der 
‘Praxis von groBer Bedeutung. Man bringt die Matrix zweckmafig in folgende Form: 


4 n 
Xy Xn 
Pia Bre Grew Pin bia 
PX PGi, 
ee eee | 
Pma ’ Pmn oles 
Pra Xq tag OS be Be 


Be schiedenen he aoe ark Nebenbedingungen, die an Stelle der Ziffern 1 bis nin dee 
werksiiblichen Nomenklatur eingetragen werden. 

Die Zeilen 1 bis m dienen zur Berechnung der gesuchten Gréfen, z. B. der benétigten Roh- 

toffe, der Durchsatze der Teilbetriebe usw.; auch sie werden namentlich an Stelle der Ziffern 1 
ibis m eingetragen. 
Die Rechenvorschrift lautet nun einfach so: ,,Man trage. in die oberste freie Zeile die vor- 
azegebenen Werte fiir den Produktionszustand, die ,,Ausgangszahlen“ (in groBen Buchstaben 
argestellt) ein. Um den Wert fiir irgendeine Zeilengréfe zu berechnen werden die in dieser 
WZeile stehenden Koeffizienten mit den senkrecht dariiber stehenden Zahlen der obersten Zeile 
tmultipliziert und in die dafiir vorgesehene Freizeile eingetragen. Die Summe dieser einzelnen 
Produkte gibt den gesuchten Wert der jeweiligen ZeilengréBe und wird in die letzte freie 
( palte eingetragen.* 


Ht 


Mit dieser einfachen Rechenregel kann einem Nichtmathematiker in kiirzester Zeit die 
“Multiplikation einer Matrix mit einem Vektor erlautert werden, so da8 er praktisch damit 
‘rechnen kann. Wenn mit einer Rechenmaschine gerechnet wird, welche die einzelnen Produkte 
‘sleich aufzusummieren gestattet, wird man in den Matrizenformularen die Zwischenfreizeilen 
Wweglassen, dafiir aber oben und rechts mehrere Freizeilen und Freispalten vorsehen, damit 
fein Formblatt fiir mehrere Rechnungen benutzt werden kann. Gibt man als ,,Ausgangszahlen“ 
adie Spaltenvektoren einer zweiten Matrix ©, vor, so kann auf diese Weise auch die Multi- 
iplikation zweier Matrizen durchgefiihrt werden. 


Wenn man eine Matrix §$ in dieser Weise anordnet, d.h. an die einzelnen Zeilen und 
‘Spalten gleich ihren ,,Namen“ schreibt, so eriibrigt sich eine besondere Erklarung, was in der 
|Gleichung =x unter den Komponenten der Vektoren t) und y zu verstehen ist. Dieser 
Vorteil ist, so trivial er erscheint, bei umfangreichen Matrizen mit etwa 50 Zeilen und Spalten, 
pres i in der Praxis keine Seltenheit ist, nicht zu unterschatzen. Ein weiterer Vorteil ist fiir die 
‘Praxis darin zu erblicken, da®& sich die Berechnung des Produktes ‘$y mit Hilfe von Loch- 
‘artenmaschinen durchfiihren 14Bt. Zu diesem Zweck wird die Matrix ‘3 in einen Stapel 
(Lochkarten abgewandelt, die nur einmal, solange die Matrix Giiltigkeit hat, von Hand gelocht 
‘werden miissen. Soll nun $¢ fiir einen bestimmten Vektor x berechnet werden, so brauchen 
ur dessen Komponenten von Hand gelocht zu werden; alle anderen Rechnungen werden 
lautomatisch durchgefihrt. Locht man in die Lochkarten der Matrix auch noch die Kompo- 
oenten des Bewertungsvektors bys (Ziff. 4) ein, so hat man die Méglichkeit, die Komponenten 
des Betriebsvektors {* (Ziff. 4) nicht nur mengenmafig, sondern auch geldmafig ausrechnen 
zu lassen und bekommt damit auch automatisch das Ergebnis G. In dieser Méglichkeit, die 
‘Matrizenrechnung durch Lochkartenmaschinen durchzufihren, liegt vielleicht ein besonderer 
Anreiz fiir gréRere Werke, dieses Verfahren anzuwenden. 


i 7, Zusammenfassung. Der Betriebsablauf eines aus vielen miteinander verkoppelten Be- 
trieben bestehenden Industriewerkes lat sich in einer der Struktur des Werkes angepaBten 
Weise durch eine Matrix beschreiben. Dadurch wird es méglich, fiir viele wirtschaftliche Auf- 


10* 
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gabenstellungen den gréften Teil der erforderlichen Rechnungen vorbereitend zu erledigen 
und deren Ergebnis in einer ,,Strukturmatrix“ zu speichern. Um einen vorgegebenen Betriebs 
fall durchzurechnen, bedarf es dann nur noch einer Multiplikation dieser Matrix mit einem; 
den betreffenden Einzelfall kennzeichnenden Vektor, eine Aufgabe, die auch von nicht be- 
sonders mathematisch vorgebildeten Kraften, ja sogar rein maschinell mit Lochkarten- 
maschinen durchgefiihrt werden kann. Die besondere Anpassungsfahigkeit des Kalkils zeigt 
sich in vollem Umfang bei der Aufgabe, aus den Matrizen der Einzelbetriebe die des ganzent 
Betriebss ystems aufzubauen. Die Nichtvertauschbarkeit der Matrizenmultiplikation und diet 
Singularitat von Matrizen finden dabei ihre sinnvolle Deutung in den betrieblichen Ver-} 
haltnissen. 

Eine Anwendung der Matrizenrechnung auf betriebswirtschaftliche Aufgaben erscheint vor 
allem da aussichtsreich, wo ein aus vielen Einzelbetrieben stammendes Zahlenmaterial unter} 
besonderer Herausstellung der kausalen Zusammenhange zusammengefafBt und fiir praktische 
Anwendungen nutzbar gemacht werden soll. Der dazu nétige mathematische Apparat ist i 
meist durchaus einfach. So zeigt es sich, daB ein von den Mathematikern schon langst ents} 
wickelter Kalkiil mit Vorteil auf Probleme der modernen Wirtschaft angewandt werden kann, 


(Eingegangen am 16. September 1952) 
Anschrift des Verfassers: Dr.-Ing. O. Pichler, Leuna (Krs. Merseburg) Wickenweg 9. 
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Die diinne Kreisplatte mit pulsierender Randbelastung 
in ihrer Mittelebene als Stabilititsproblem. 


Von K, A. Reckling. 


1. Einleitung und Problemstellung. Der Einflu8 einer konstanten Randbelastung auf die 
| Biegeschwingungen von diinnen Kreisplatten wurde zuerst von Federhofer! untersucht. Kommt 
t noch, wie es in der vorliegenden Untersuchung angenommen wird, eine periodisch verinder- 
liche Randbelastung mit vorgegebener Frequenz hinzu, so tritt die Frage in den Vordergrund, 
| bei welchen zusammengehérigen Werten von Randbelastung und Frequenz die erzwungenen 
1 Radialschwingungen der Platte instabil werden, bzw. wann die durch eine kleine Stérung 
* eingeleiteten Biegeschwingungen mit der Zeit immer starker anwachsen. Dieses ist die iibliche 
Fragestellung bei allen Stabilitatsproblemen der Elastokinetik, von denen als einfachstes 
| Beispiel der durch pulsierende Lingskrafte belastete gerade Stab mehrfach behandelt wurde ®. 
| Das in einer fritheren Arbeit des Verfassers? zur Lisung eines dynamischen Stabilitatsproblems 
benutzte Verfahren der Stérungsrechnung fiihrt auch 
_ im vorliegenden Falle zum Ziel. 
Die Platte sei nach Abb. 1 am Rande durch eine in 
_ ihrer Mittelebene wirkende, allseitig gleichmaBige Radial- 
kraft P, + P, cos wt belastet, wobei der statische, als 
Druck positiv gerechnete Anteil P, und die Amplitude 
P, des mit der vorgegebenen Frequenz w schwingenden, 
zeitlich veranderlichen Anteils je Langeneinheit des 
| Umfangs gerechnet sind. Die Platte soll am Rande in 
gleicher Hihe gestiitzt und dort in beliebiger Weise ein- 
/ gespannt sein. Thre Dicke A sei klein im Vergleich zum 
Radius R. Der Einflu8B der Dampfung wird bei der 
Untersuchung nicht beriicksichtigt. 

2. Die Differentialgleichung der Biegeschwingungen 
und der Ansatz zu ihrer Lésung. Wenn r und 9 als Polar- 


P,+P,cos@t 


koordinaten eines Plattenpunktes eingefiihrt werden, 2 hot ih mare 
. oe . . . 4, 
mit m die Masse der Platte je Flacheneinheit und mit ZN : ae = pase, 
= 3 = ° “Ge . . . n= 
N= Gh /6(1 ”) ihre Steifigkeit bezeichnet wird, a Abb. 1. Kreisplatte mit pulsierender Radiallast. 


lautet die Differentialgleichung fiir die Auslenkung 
w(r, y, t) aus der urspriinglichen Plattenmittelebene', wenn die konstante Belastung durch 
die zeitlich veranderliche ersetzt wird, 


a wi Aw += (Py + P, cos wt) dw=0, (1) 
| -wobei zur Abkiirzung der Operator 
) i era y Gas 1 2 
US or® 1 r Or | r? dy? = ira dy? c 
 eingesetzt wurde. Nach Einfithrung von 
, Sn Re t 
wo, eet = Oh, i=: tO. 
(3) 
NB, ps Nelad coe 
far ip fin ip,? Ler K 


und Px als erster statischer Knicklast der Platte geht die Differentialgleichung (1) der Biege- 
schwingungen itber in die dimensionslose Form 


4 + AAn + pw? (e& +e, cosat) Ay = 0, (4) 


1 K. Federhofer, Ing.-Arch. 6 (1935) S. 68. 
2 Vel. Schrifttumsangaben der in Fuf$note 3 genannten Arbeit. 
3K. A, Reckling, Ing.-Arch. 20 (1952) S. 137. 
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worin Ableitungen nach t mit Punkten bezeichnet sind und der Operator A(0, p) entsprechend 
(2) definiert ist. Diese Differentialgleichung ist unter Anpassung an die jeweiligen Rand- und 
Anfangsbedingungen zu integrieren, und schlieBlich ist zu untersuchen, bei welchen Frequenzen 
o« und Belastungen ¢ die Auslenkungen 7 mit der Zeit tT anwachsen. 

Aus der Entwicklung der Lésung von (4) nach wachsenden Potenzen des Parameters ¢, 


1 (QP t) = 2 ef Mn (0, 9 7) (5) 

folgt das System ai 
Ngo + AAny + w? & AN = 0, (6) 
Hn + AAnn + & Ayn = — wcosat An, (nm =1,2,..., m) (6a) 


von m +1 partiellen Differentialgleichungen, dessen Kinzelgleichungen sukzessive zu lésen sind. 


3. Freie Biegeschwingungen. Durch Gleichung (6) werden die freien Biegeschwingungen der 
Platte beschrieben. Fiir den Spezialfall 7, = 0 und ¢, = 1 stellt (6) die Differentialgleichung 
der Plattenknickung dar; ihre Integration ergibt bei Anpassung an die jeweiligen Rand- 
bedingungen die Kigenwerte 4 und damit die Knicklasten, von denen fiir die vorliegende 
Untersuchung nur die erste interessiert und fiir das folgende als bekannt angenommen wird. 
Auch die Lésung der vollstandigen Gleichung (6) ist bekannt!, soll hier aber wegen ihrer 
Bedeutung fiir die anschlieBende Rechnung kurz wiedergegeben werden: Der Ansatz 


foe} 


1o (0: 7) = 3S Vaeel@- #) Paro) (1) 


p= 
liefert aus (6) fiir die Zeitfunktionen die Differentialgleichungen 


Do kp ie kp Do, kp =0 
mit den Integralen 
Do, ip = Akp COS Wupt + By sin Wx, T , (8) 
welche beliebigen Anfangsbedingungen angepaf®t werden kénnen, sowie fiir die Ortsfunktionen 
die Differentialgleichungen 
ANU ip a i ee AU gt "Onn Ue (9) 
die durch 
Uo, kp (Q> P) = Uo, kp (2) sin (pp + Hp) (10) 
und die gewohnlichen Differentialgleichungen 
wt ’ 2, 2 4 1 , a fe a 4 c 
(2 Uo,np)”” 4 ( P + bP 00) M6. 4 08 £ 


Q 


erfiillt werden. Letztere lassen sich mit Hilfe der symbolischen Schreibweise nach (2) in die | 
Form 


De 2 
[Pe ~ 00bp) ap =O (1) 


2 2 2 | 

| D(0) — (ftp +%5)| [| B10) + (ttp —2)} mao] = 0 (1a) 

mit | 

i 1 ee 249 | 

i} <2 (GPa 40h; F whey), hp = Mbp, (1b) 
bringen, also in zwei Gruppen von Besselschen Differentialgleichungen aufspalten, von deren 
Integralen nur J,(A,,@) und J,(iA,po) als Lésungen in Frage kommen, da die Durchbiegung 
fiir @ =0 endlich bleiben mu. Wenn man die durch die Beziehung T, (Axp ©) =P Jp (i Any 0) 
definierten modifizierten Besselschen Funktionen einfiihrt und beachtet, daB die Platte am 


Rande in gleicher Hohe gestiitzt, also U,kp (1) =O sein soll, erhalt man die Kigenfunktionen 
der freien Biegeschwingungen 


Tp (Akp) PA sn) ( ) 


mit den beliebigen Konstanten C,,. Der Index k bezeichnet die Anzahl der Knotenkreise, 
wahrend p die Zahl der Knotendurchmesser angibt. 


Uo, kp (0) oo Crp De (Akp 0) ae 


1 Vel. z. B. FuBnote 1 auf Seite 141. 
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4. Frequenzgleichungen; Orthogonalitat der Eigenfunktionen. Die Kigenwerte A, und Anp 
— und damit nach (lla) die Eigenfrequenzen Wp — erhalt man aus Frequenzgleichungen, 
deren Gestalt sich nach den Einspannbedingungen am Plattenrande richtet. Ist die Platte 
am Rande elastisch eingespannt, so ist dort die Bedingung 1 


Uo, wp(L) + (v + €) up, xp(1) = 0 (13) 


zu erfiillen, wobei c den Einspannungsgrad angibt und die Striche Ableitungen nach @ be- 
deuten. Fithrt man (12) hierin ein, so erhalt man als allgemeine Frequenzgleichung 


aS (2 c) | Dep Jp (Axp) Iya (ep) eae dnp Jp-1 (Axp) Ip (xp) ] = 
= (Akp + ip) J (Arp) I, (Arp) =0, 


i 

vaus der sich leicht die Gleichungen fiir spezielle Einspannbedingungen herleiten lassen; z. B. 
fiir freie Auflagerung mit c=0 und fiir vollkommene Einspannung mit 1/e =0. Wenn man 
reine elastische Einspannung ¢=1—vy vorschreibt, also mit der. Einspannungsbedingung 
Uo, kp (1) + Uo, ep (1) = 0 das Verschwinden der mittleren Kriimmung am Plattenrande ver- 
langt, so erhalt man die einfache Frequenzgleichung J,(A,») = 0, die als Naherungsgleichung 
fiir die freie Auflagerung verwendet werden kann (vgl. Ziff. 7). Die beiden Eigenwerte sind 
nach (11b) durch die Beziehung 2%, —Ajip = we) miteinander verkniipft und werden fir 
verschwindende statische Vorbelastung gleich. 

Fiir die spatere Rechnung kommt es auf die Orthogonalitat der Eigenfunktionen an. Das 
Kigenwertproblem (11) 148t sich in der abgekiirzten Form L[u,,] = @ wip Upz,) mit den Rand- 
bedingungen (13) schreiben, wobei der Differentialausdruck L[u,;,,] selbstadjungiert ist. Aus 
dem fiir zwei Eigenfunktionen up, und up,n, gebildeten Integral 


(14) 


: : 
f (Uo, kp L [Upan | cae tg np Lug, kp l) do 


e=0 
erhalt man durch Teilintegration 
; 


: 


‘Da im dritten Randteil auf der rechten Seite 


1 
2 2 pe: ” " 
(kp ca np) iy QU, kp Uo, np do = [uo, kp Up, np Uo, np Uo, te lo=6 ae 
e=0 


Vf / 
Uo, kp “o,np — “0, np “o, kp 
Y o=0 


" / ” Hf 
— [Uo, kp Uo, np — Uo, np U0, Lalomt ia 


[uo, kp Uo, Rp? om Uo, np Up, poo =0 
sein mu, hebt er sich gegen den ersten fort, wahrend der zweite wegen der Randbedingung 
(13) verschwindet. Das Eigenwertproblem ist daher fiir alle Randbedingungen und alle 
Schwingungsformen selbstadjungiert, d.h. die Eigenfunktionen sind mit der Belegungs- 
funktion g orthogonal und kénnen auferdem durch passende Wahl der hbeliebigen Kon- 
eenten C,p in (12) normiert werden. Somit gilt die Orthogonalitatsbezichung 
) 4 ( 0 fir kAn 

ang 15 
Lieto (0) Uonp (0) de od far Pen: (15) 


5. Erzwungene Schwingungen. Nach diesen Vorbereitungen kénnen die Gleichungen (6a) 
sukzessive integriert werden. Die erste Gleichung fiir n =1 lautet mit (7), (10) und (12) 


nt 4dr + weg Am = >) Cop bp Jp Aap) + | “ 
p 1 
egy sell) he 2) sin (pp +9) (D_1, xp (t) + Oy, xp(7)) » | 
Ip (Akp) 


wobei zur Abkiirzung fiir die Zeitfunktionen 
Dz kp (T) = App 60s (p+ 90) T+ Byp sin (xp + 74) T (q=0, +1, +2,...) (17) 
gesetzt wurde. Mit dem Ansatz 
m1 (0; P> tT) = = > [u_a.up (0) P-1, 4p (T) + U4, ep (Q) P41, xp (7) ] sin (PP + %) (18) 
7, 


1 Vgl. z. B. Biezeno-Grammel, Techn, Dynamik, S. 440. 
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kommen aus (16) die gewéhnlichen Differentialgleichungen 
2 2 
D0) oy Gz +%| [Po cr (ctu pe | ust, uh = | 


/\\ 
_ Crp; 72 Jp (Akp) a. 
= OF ly Jo (tap) + Fp ee Ty (iup0)| 
J 


Tap = > (Vu ea)? F 4 (my F 0)? + He) » | 


(19) 


(19a) 


mit 


Orkp = 5 VW Eo)” + 4 (Oep F 0)? — pw? £6) ; | 


Die homogenen Gleichungen entsprechen denGleichungen (lla) und haben daher als Integrale 
die Besselschen bzw. die modifizierten Besselschen Funktionen mit den Argumenten 0+;)0 
bzw. G34p@. Eine Partikularlésung der inhomogenen Gleichungen kann leicht durch einen 
der rechten Seite entsprechenden Ansatz gefunden werden, so daf§ man schlieSlich als voll+ 
stindige Lésung von (19) 


ux, kp (0) = Crk p IE (Cxxp 0) 4 Cxkp is (Oxkp 0) a Orkp a (Arp 0) — Oxkp I (Arp 0) (20) 


mit den Konstanten 


242 © are 
ae LU kp ~“kp = Le Akp Crp (20a) 


= = 3 a= = = = 
"PT 2 [Ay — # ey Aky — (Ory F 2)? "PO [dtp + pt ag dip — (cup € 0)F 
und den durch Anpassung der Liésungen an die Randbedingungen zu ermittelnden Konstanten | 


Crkp» Crp erhalt. Damit ist die erste Naherung 7, bekannt und kann auf der rechten Seite 


der Gleichungen (6a) fiir n =2 eingesetzt werden: 

- 2 = _— — 

UP ae AAn, ae We An, = — pe > {[e_np Cos ak (OE. 0) wn C—kp Oke is (x5 0) ae 
kp 


== Cpp Ot bpd p (O.%p0) — C4.kp O+ kp T, (O+1p@) <5 . 
(a _ap + O14 4p) Akp Jp (Ap 0) — (& ap + &4.kp) Akp Tp (Ap 0) 1p, up (0) + @y) 
+ [e+] Do ap(t) +L ++.) By op Wf sin (PP +%) = 


a oe {Frep (Q) Po, xp (t) + ***"} sin (PP + Go) - 


Hier sind nur die bei Zeitfunktionen mit q=0 stehenden Ortsfunktionen hingeschrieben, da 
sie allein fiir eine Naherung bis ej einschlieSlich in Betracht kommen. Die F,,(0) sind bekannte 
stetige Funktionen von 0, die nach den Eigenfunktionen der freien Schwingungen entwickelt 
werden. Da Letztere nach Ziffer 4 ein orthogonales Funktionensystem bilden, kann die Ent- 
wicklung in der Form 


Fry (Q) = & af) Uo, np (0) (22) 
mit den Entwicklungskoeffizienten 
1 
a(n) = J ¢ Fev (0) to, mp (@) de (22a) 
vorgenommen werden, womit die Gleichungen (21) itibergehen in 
ee 2 
2+ AAn,+ we Ans = - > ee oh al”) Uo, np (0, P) Po,kp(t) + - (23) 
k p n 
Hieraus liefert der Lésungsansatz 
Me = 3 XS Of) Us, np (0, 9) Of) (0) (24) 
P n 


zunachst 


FETA am DB") + (AA Ug np + p? & AU, np) BOY] = | 


“ESS Satie ee 9) 
a2 


» was fiir n ~k durch 
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Dy) = Dy, up (26) 
erfiillt werden kann, wodurch (25) mit Beachtung von (9) tibergeht in 


fi 2 / 
n 2 2 Nad a 
= 2 Ch (Onp — kp) Ug, np = Gy > >, al”) Up, np « 
Pie Pp n 


} Die Striche an den Summenzeichen deuten an, daf® die Glieder mit n—k fortzulassen sind. 


Der Vergleich der linear unabhangigen Funktionen Uy np (@. ~) liefert die Entwicklungs- 


i koeffizienten 
G2 er 
Ce (ot, = of) al” (n Ak). (27) 
| Fir n=Ef versagt der Ansatz (26). Fir die Zeitfunktionen gelten jetzt die Differential- 
_ gleichungen 
= - 2 (*) 
k A 
OY) + a}, Of) = _ o: One (28) 
Ch 


Da die Lésungen der homogenen Gleichungen dieselben Frequenzen zp Wie die Zeitfunktionen 


®),xp haben und daher der .,Resonanzfall*+ vorliegt, 14Bt sich (28) erfiillen durch 


OY (t) = t(Agp sin Opp T — By cos Wpp T) (29) 


) und 


(eA) 
Os, og ay (30) 


| Damit sind simtliche Anteile der Lésung 7, gewonnen und die Einzellésungen 7, und 7, werden 
) nach (5) zur Gesamtlésung 


= 2 D> {pe tel, ef) Up np) (Ap. COS Wy, tT + B,, sin @, Tt) + | 
5 ; (31) 
+ ¢ cf) Up, kp T (App Sin Opp, T — By, cos Mp T)} sin (PY + Pp) | 


» zusammengesetzt, wobei der fiir die Naherung bis ¢{ unwesentliche Anteil 7), fortgelassen wurde. 


6, Stabilitatsuntersuchung. Zundachst wird die Lésung (31) umgeordnet zu 


4 = 2X [Avy BL) (QP 1) + Bey 8h) (0, P 7) ] = 2 Me » 
P P 


' indem alle bei den Konstanten Akp bzw. By, stehenden Ausdriicke zusammengefafBt werden. 


Fiir jede Kinzellésung 7, ist sodann zu untersuchen, wann sie mit der Zeit anwachst. Die 
Beantwortung dieser Frage liefert ein bei solchen Problemen zuerst von W. Kucharski an- 
gewendetes, aus der Theorie der gewéhnlichen Differentialgleichungen mit periodischen 
Koeffizienten auf partielle Differentialgleichungen iibertragenes Verfahren’ *: Man bildet die 


neuen Funktionen a) _ st) (c+ P) mit P= und entsprechend ai) und beriicksichtigt nur 
ihre Zeitabhangigkeit. Da diese Funktionen auch Lésungen der Differentialgleichung (4) sind, 
lassen sie sich durch lineare Substitution 

zl) = al) af) + aly gf) = 1,2) (32) 
aus den g() aufbauen, die noch bestimmten Anfangsbedingungen unterworfen werden, so da 
die Floquetsche Theorie anwendbar wird. Die Bedingung fiir die Stabilitat der Lésungen 7, 
lautet dann ‘ 

| ott?) | zak. (33) 

Die Substitutionskoeffizienten bestimmt man aus (32) in gleicher Weise wie in der friiheren 
Arbeit des Verfassers2 durch Vergleich der bei gleichen Zeit- und Ortsfunktionen stehenden 
Ausdriicke. Fiihrt man noch das Frequenzverhaltnis 


Sup =—- = : (34) 
kp rp Arp /23, == i a, 


1 W. Kucharski, Ing.-Arch. 18 (1950), S. 385. 
2 Vel. FuBnote 3 auf Seite 141. 
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ein und nimmt in (33) das Gleichheitszeichen, so liefert die Ausrechnung fiir jede Schwingungs- 
form eine Schar von Grenzkurven zwischen stabilen und instabilen Lésungsgebieten in den| 
(e(*P), d,,)-Ebenen aus den Gleichungen 


ae ie ONG | 
a aie ely [2 bmp (4 1 — 008 : 
o(kp) — dip 72> He eq Okp : (35) 
x - 


5 Das 
Tt al®) Sie 
P. Okp 


Die Instabilititsgebiete nehmen ihren Ausgang von den durch e(*?)=0 gegebenen Punkten | 
der 6,,-Achsen, also von 04,=2/r bzw. von den kritischen Frequenzen 
} 

2 On 


Xkrit = = 


Obwohl (35) bei Einsetzen dieser Werte einen unbestimmten Ausdruck von der Form 0/0_ 
liefert, ist doch 
lim e(*P) = 0, 

okp > 2/7 


Weitere kritische Frequenzen treten bei einer Naherung bis ef nicht hinzu, da die Ent- | 


wicklungskoeffizienten a\") wegen (22a) und (21) nur mit den Koeffizienten (20a), d. h. nur_ 


fiir «a =—20,, gegen unendlich gehen. 


Das Problem ist damit in erster Naherung fiir alle Randbedingungen gelést, wenn die | 


Entwicklungskoeffizienten af") aus (22a) berechnet sind. 


7. Ein Beispiel. Der Kinfachheit wegen wird das Beispiel einer Platte behandelt, deren 


mittlere Kiiimmung am Rande veischwinden soll. Nach Ziffer 4 erhalt man den Eigenwert | 


App aus der k-ten Nullstelle der Besselschen Funktion p-ter Ordnung. Berechnet man das 
allgemein zu u’’(1)-++ vu’(1) proportionale gréBte Randmoment M, fiir diesen Einspannungsfall 
und vergleicht es mit dem bei vollkommener Einspannung (1/c =0) und gleicher Auslenkung 
u(0) in Plattenmitte auftretenden Randmoment M,, so findet man fiir die Grundschwingungs- 
form k—0; p=O der nicht statisch vorbelasteten Platte (¢)=0) M,=0,14 M,. Fiir die 
statisch vorbelastete Platte sowie fiir alle anderen Schwingungsformen wird M, sogar noch 
kleiner im Verhaltnis zu M,. Die gewahlte Randbedingung kann also als brauchbare Naherung 
fiir die frei aufgelagerte Platte angesehen werden, da sowieso in der technischen Praxis keine 
vollkommen freie Auflagerung vorkommt. 

Die Rechnung gestaltet sich nun besonders einfach, weil die Funktionen J, (xp @) in samt- 
lichen Gleichungen als einzige Ortsfunktionen auftreten. So erhalt man 


U4, kp (2) = Oxkp Jf (Akp 0) mit O=+kp 
nach (20a) und 


li (2) = hip (G2 ais X + kp) Jp (Axp 0) = al) Uo, kp (0) a 


In der Entwicklung (22) ist also a) —0 fiir k An, und die fiir die Stabilitatskriterien 


mafgebenden Koeffizienten ai) lassen sich leicht berechnen: 


(ie Mep (=p FO hp) wan, 
kp Crp (ip — &4) (4 — ip) 


Wenn man noch den ersten Eigenwert js = Ag) des statischen Knickproblems einfiihrt, erhalt 
man damit aus (35) die Beziehung 


22 
jie De), (O% =A) (+ il + cos") 
aa 7) 0400 kp \ kp ; 0 
el 2) ee ees | ee (36) 


a 


- oo SAGE 

7 sin —— 

Okp 
fiir die Grenzkurven zwischen stabilen und instabilen Lésungsgebieten. Erreicht die statische 
Vorbelastung die Knicklast, d.h. wird ¢)=1, so ergibt die erste Gleichung ¢\°° = 0, also 


gleichfalls Instabilitaét in allen Fallen. Die Engebnisse der fiir verschiedene statische Vor- 
belastungen und die drei ersten Schwingungsformen durchgefiithrten Auswertung mit 
oy = 2,4048, Ay) = 3,8317 und Ay) = 5,5201 sind in Abb. 2 dargestellt. Die zugehérigen di- 


| 
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bei der ersten Schwingungsform am 
ausgepragtesten. Die sich links von 
Orp = 2/3 haufenden Gebiete sind 
nicht eingetragen, da sie praktisch 
bedeutungslos sind. Fiir die erste 
Schwingungsform ohne _ statische 
Vorbelastung stimmt (36) iibrigens 
genau mit dem fiir den langspulsie- 
rend belasteten, frei aufgelagerten 
(Stab ermittelten Stabilitatskrite- 
rium! iiberein. 

Fir die héheren Schwingungs- 
formen interessiert praktisch héch- 
stens noch der Frequenzbereich in 
der Nahe von «=2«a,,, der aber 
schon fiir die Schwingung mit einem 
Knotenkreis sehr schmal wird. Die 
‘dargestellten Bereiche werden daher 
fiir die praktische Beurteilung des 
)Stabilitatsverhaltens im allgemei- 
nen ausreichen, besonders wenn der 
)stets vorhandene stabilitaétsverbes- 
isernde EinfluB der Dampfung mit 
beriicksichtigt wird. 
| Es ist von besonderem Interesse, 
}dafs Formel (36) auch Instabilitats- 
gebiete fiir eine statische Vorbela- 
istung auf Zug (¢)<0) liefert: Fiir 
fden Fall k —0, p=0 und ¢, =— 0,5 
‘ist das von by) =2 ausgehende noch 
‘stark ausgepragte Instabilitats- 
igebiet in Abb. 2 eingetragen. Die 
‘Plattenschwingung ist also selbst 
dann noch instabil, wenn die kon- 
istante Zugkraft P, gréBer als die 
Amplitude P, der pulsierenden Last 
iwird. An diesem Fall kann man 
besonders deutlich erkennen, wie 
ein auf den ersten Blick durchaus 
stabil erscheinender Schwingungs- 
zustand doch bei gewissen Fre- 
“quenzen instabil werden kann. 


Weel 


1 Vel. FuBnote 3 von S. 141. 


2 2 
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Mit wachsender Vorbelastung ¢) werden die Instabilitatsgebiete breiter und sind ferner 


mensionslosen Kigenfrequenzen folgen fiir verschiedene Vorbelastungen ¢) aus 
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Abb. 2. Instabilitiitsgebiete fiir die drei ersten a aplan a Sake e 
und verschiedene statische Vorbelastungen P, = e, P K: 
[ fiir e, = 0 (schraffiert) ; fiir €, 1/4 -—: fiir e,=1/2; 
O—O-— fir e, 2 (Zug- MWeckclasenueilt 
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Gewasserkunde. Von Dr. Walter Wundt, Honorarprofessor an der Universitit Frei- 


burg i. Br., Mitglied der Deutschen Akademie der Naturforscher in Halle a. d. S. Mit 
185 Abbildungen. VI, 320 Seiten. 1953. Ganzleinen DM 34.50 


Inhaltsibersicht: Einleitung. — A. Der Kreislaut des Wassers. — B. Die bautechnische Problemstellung: 
Wasserschutz. ‘Wassernutzung. — ©. Gewasserkunde und Klima: Die meteorologischen Elemente, Klimatypen. 
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und Gerdlle. — E. Grundwasser, Quellen und Boden: Grundwasser und Quellen. Bezichungen zum Boden. — 
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Pa Mensch. Naturschutz. — G. Messung und Statistik in der Gewasserkunde: Direkte Messung. Indirekte Messung, 

tatistische und analytische Methoden. — H. Die unregelmaBigen Schwankuhgen der Wassertiihrung: Allgemeines. 
Hochwasser. Niedrigwasser. Klimatische Schwankungen. — J. Die regelmaSigen Schwankungen und die mittlere 
Wasserfiihrung: Die regelmaBigen Schwankungen. Die mittlere Wasserfiihrung. — K, Vorratsbildung und Wasser- 
wirtschaft. — Anhang: Tabellen. — Literaturverzeichnis. — Orts- und Sachverzeichnis. 
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_ Geschaftsfiihrer der Deutschen Forschungsgesellschaft fiir Bodenmechanik (Degebo), 

Technische Universitat Berlin. Mit 498 Abbildungen. XI, 464 Seiten. 1950, 
Ganzleinen DM 43.50 


Inhaltsiibersicht: I. Teil: Voruntersuchungen vor Baubeginn: Arbeitsprogramm und Gang der Unter- 
suchungen. — Untersuchungen im Geliande: Flichenuntersuchungen des Untergrundes. Punktuntersuchungen 
des Untergrundes. Untersuchungen des Grundwassers. — Untersuchungen im Laboratorium: Normalversuche. 
Besondere Untersuchungen. — IJ. Teil: Untersuchungen am fertigen Bauwerk: Beobachtungen und Messungen 
an Bauten ags Erde: Nachpriifung der Verdichtung. Setzungen. Erschiitterungen. — Beobachtungen und 
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und Spannungen. Risse. — Beobachtungen und Messungen des Grundwassers: Grundwasserstand, Auftrieb. 
Grundwasserstrémung. — SchluBbemerkung. — Literatur- und Sachverzeichnis. 
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Handbuch des Wasserbaues. Von Dipl.-Ing. Dr. techn. Dr.-Ing. he. A. Schok- 


litsch, Professor der Universidad Nacional de Tucuman, Argentinien, In zwei Banden. 
Zweite, neubearbeitete Auflage. ' 


Erster Band: Mit 722 Textabbildungen und 87Zahlentafeln. X, 478 Seiten. 
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Inhaltstiibersicht: Meteorologie. — Gewisserkunde und Hydraulik. — Die Wasserversorgung. — Die 
Ortsentwasserung. 


Zweiter Band: Mit 1327 Textabbildungen und 26 Zahlentafeln. VIII, 594 Seiten. 
1952. (W) 


Inhaltsibersicht: Stauwerke und Entnahmeanlagen. — Wasserkraftanlagen. — Meliorationen, — Der 
Flu8bau. — Verkehrswasserbau. Sachverzeichnis. : 


Das Werk wird nur noch geschlossen abgegeben, 
Preis fiir das Gesamtwerk, zwei Bande in Ganzleinen: DM 155.40 


Der Grundbau. fin Handbuch fiir Studium und Praxis. Von Dipl.-Ing., Dr. teehn., 


Dr.-Ing. h. c. Armin Schoklitsch, Professor der Universidad Nacional de Tucuman, 
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Inhaltstibersicht: f. Allgemeines zur Elektrizititsversorgung. — Il. Krattwerke: Wirmekraftanlagen. 
Wasserkraftanlagen. Der Einflu8 des zeitlich veriinderlichen Verbrauchs auf die Kraftwerke. Verhalten der ; 
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X. Freileitungen: |}Aligemeines. Die Seilschwingungen. Isolatoren fiir Freileitungen. Maste und Leitungsanord- 
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nnd der elektrische Schaltfehlerschutz. — XV. Fernsprechen, Fernmessen und Fernsteuern: ernsprecher. 
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veristelten Leitungen. Leitungsberechnung und Leitungsbemessung nach der Schnittmethode. Die Berechnung 

des Spannungsabfalls bei Dreiphasenstrom unter naherungsweiser Berticksichtigung der Induktivitét und der 
Phasenverschiebung. Niederspannungsmaschennetze. Die Berechnung der Induktivitat und Kapazitét von 
Netzen. Die Berecbnung von Wechselstromnetzen unter Beriicksichtigung der Induktivitat. Berechnung des 
wirtschaftlichen Leitungsquerschnitts. Berechnung der Ubertragungsspannung einer Fernleitung, Ringleitungen. 
ZusammenschluB von verschiedenen GroSversorgungen zur Verbundwirtschaft. Die Hochspannungsleitung ohne 
Spannungsabfall. Leitungen groBer Linge. Die Stabilitat von Leitungen. Uber die Méglichkeiten der Gleichstrom- 
Hochspannungsiibertragung. — XX, Berechnung der KurzschluBstréme:. Die Berechnung des Dauerkurzschlu6- 
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Elektrische Schaltvorgiinge in geschlossenen Stromkreisen von Stark- 


stromanlagen. Von Reinhold Ridenberg, Dr.-Ing.; Dr.-Ing. e. h.; A. M. hon. 


Gordon McKay, Professor fiir Elektrotechnik an der Harvard University, Cambridge, 
Mass., USA. Vormals: Chef-Elektriker der Siemens-Schuckertwerke, A.-G., Berlin. — 
Honorarprofessor an der Technischen Hochschule Berlin. Beratender Ingenieur der — 
General Electric Co., Ltd., London. Vierte, vermehrte Auflage. Mit 735 Abbildungen 
und einer Tafel. XII, 599 Seiten, 1953. Ganzleinen DM 67.50 


Inhaltstibersicht: A. Lineare Ausgleichsvorgange in geschlossenen Stromkreisen. Einfache Strom- 
kreise: Einschalten und Abschalten von Stromkreisen mit Selbstinduktion, Ladung und Entladung von 
Kapazititskreisen. Allgemeine Schaltgesetze. Resonanzerscheinungen. Ausgleichsstréme in Schwingungskreisen. 
Kinschalten von Schwingungskreisen. Schaltschwingungen beim Abschalten. -Magnetisch verkettete 
Stromkreise: Ruhende Stromkreise mit Wechselinduktion. Schalten yon Transformatoren. Wirbelstréme 
in massiven Magnetkernen. Kommutierungsflu$8 in Kollektormaschinen. Freie Drehfelder in Mehrphasen- 
maschinen. Plétzlicher Kurzschlu8 von Drehstrommaschinen. StoSkurzschluBstréme in der Praxis, Mechanische 
und thermische KurzschluSwirkungen. LEigenschwingungen in Kollektormaschinen. — Wirkung von 
Schwungmassen: Anlauf von Motoren. Einschalten von Gleichstromankern. Anlauf yon Asynchron- 
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